facility form coa 


THE SEE SITIVITY OF BUCKLING OF AXIALLY COMPRESSED 
FIBER-REINFORCED CYLINDRICAL SHELLS TO 
.SMALL GEOMETRIC IMPERFECTIONS 



June 1969 




THE SENSITIVITY OF BUCKLING OF AXIALLY COMPRESSED 


FIBER-REINFORCED CYLINDRICAL SHELLS TO 
SMALL 'GEOMETRIC IMPERFECTIONS 

Michael F. Card 


ABSTRACT 


An index of imperfection sensitivity is obtained by investigating 
the character of the initial postbuckling region for a geometrically 
perfect cylindrical shell. Using techniques of perturbation theory, 
equations governing postbuckling. behavior of a multilayered, orthotropic 
cylinder are developed. The theory considers the effects of nonlinear 
deformations induced by loading prior to buckling. Solutions for post- 
buckling displacements are obtained by converting the governing equations 
and boundary conditions to a system of second-order differential equa- 
tions, which are then written in terms of finite differences at stations 
along the length of the cylinder. The resulting algebraic equations are 
solved by matrix algebra exploiting Gaussian elimination. A computer 
program was developed to solve the resulting large systems of simultan- 
eous equations and to perform the numerical differentiation and integra- 
tion necessary to calculate the imperfection sensitivity index. 

Relative values of the imperfection sensitivity index are investi- 
gated for three types of fiber -reinforced cylinders: helically wrapped 

glass-epoxy and boron-epoxy cylinders, and aluminum cylinders overwrapped 
with boron-epoxy. The helical wrap angle was varied in cylinders of 
fixed geometry to identify wrapping configurations of minimum imperfection 



sensitivity. Sensitivity indices were obtained from the present theory 
as well as from an extension of existing theory. Agreement between 
sensitivity indices predicted by the two theories was reasonably good. 
Computed results for buckling of the three types of cylinders suggest 
that substantial differences (up to 27 percent) can exist between classi- 
cal and consistent theory predictions for buckling., depending on the 
wrapping configuration. The differences are a result of load-induced 
prebuckling deformations which are retained in consistent buckling 
theories but are omitted in classical theories. 

Sensitivity indices computed for the three types of cylinders 
suggest that 4-5° waps in helically wound glass-epoxy and boron-epoxy 
shells are desirable. For wraps possessing small imperfection sensi- 
tivity, boron-epoxy cylinders are less imperfection- sensitive than glass- 
epoxy cylinders. The most attractive configuration identified in the 
study was a circumferentially overwound boron-aluminum shell which was 
found to have a positive postbuckling slope. 

The behavior of fiber-reinforced cylinders was investigated experi- 
mentally by conducting compressive buckling tests on twelve 50-inch- 
diameter glass-epoxy shells with variable wrap configurations. Elastic 
constants computed from test data were in good agreement with analytical 
predictions. Depending on wrap configuration, compressive buckling loads 
varied from 50 to 90 percent of consistent theory predictions. Compari- 
sons of imperfection sensitivity indices with test data show that the 
indices are in fair agreement with experimental buckling data trends. 

The theoretical and experimental studies reported herein suggest that 



certain fiber- orientations can enhance the compressive buckling strength 
of geometrically similar cylinders due to the fiber configuration's 
reduced sensitivity to geometric imperfections in the cylinder wall. 
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IV. INTRODUCTION 


The purpose of the dissertation is to investigate the sensitivity 
of fiber -reinforced cylinders to small imperfections in the geometry of 
the wall of the cylinder. The development of boron and high-modulus 
glass fibers as reinforcement suggests that fiber-wound or fiber- 
reinforced shells possess great technological potential. Recent experi- 
ence at NASA (refs. 1 to and elsewhere has indicated that the state 
of the art of these new materials is such that their elastic behavior 
can be predicted with reasonable engineering accuracy. Furthermore, 
efficiency studies of cylindrical shells loaded in axial compression 
(ref. 5) have shown that they are competitive with the best metallic 
shells as lightweight load-carrying members. Thus, studies of filamen- 
tary shells are of current technological interest- 

If fibrous shells are ever to be accepted as reliable load-bearing 
elements, their structural behavior must be thoroughly understood. One 
of the most important factor's in the buckling behavior of a shell is its 
sensitivity to small deviations from perfect geometry. The sensitivity 
of axially loaded cylindrical shells to initial imperfections has long 
been suspected to be a major cause of lack of agreement between theoret- 
ical predictions and experimental buckling data for shells loaded in 
axial compression. Recent tests on carefully manufactured cylindrical 
specimens (refs. 6 to 8) indicate that shells of near-perfect geometry 
can develop loads which approach predictions for idealized shells. The 
problem, however, still remains to find some method of making practical 
shell structures more predictable and, hence, more reliable. 
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Hy the nature of their construction, fiber-reinforced shells usually 
possess orthotropic material characteristics. In studies of the behavior 
of orthotropic shells, there is some evidence that, given the proper 
orthotropy and geometry, orthotropic shells can- behave with relatively 
low sensitivity to small initial imperfections. For example, buckling 
data for axially loaded longitudinally stiffened cylindrical shells 
(which can be considered to be a class of orthotropic shells) have been 
shown by the author to be in good agreement with theoretical predictions 
for ideal cylinders (ref. 9). Furthermore, postbuckling studies (refs. 10 
to 14) have indicated that minimum postbuckling loads for orthotropic 
cylinders can be much higher than isotropic cylinders with similar radius- 
thickness ratios. These observations have been used in attempts to 
predict "knockdown factors" (reductions in theoretical predictions) to 
correlate buckling theory and experiment (ref. 15) * The insensitivity 
of orthotropic shells has not been established in general and, in fact, 
imperfection sensitivity of stiffened shells in certain geometric ranges 
has been noted in references 1 6 and 17 • 

An interesting peculiarity of fiber-reinforced shells is the fact 
that by suitable layering, two-dimensional isotropic relations for the 
elastic constants of the shell wall can be obtained. Thus, by varying 
the orientation of fibers in a shell of fixed geometry, one can encounter 
for a shell with specified length, radius, and total, thickness, an iso- 
tropic configuration among configurations of varying degrees of ortho- 
tropy. Hence, the designer of fiber -reinforced shells, as opposed to 
metallic shells, has at his disposal an additional degree-of -freedom — 
orthotropy. 
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Based on the limited evidence of the reduced sensitivity of certain 
orthotropic shells to initial imperfections, we are confronted with the 
following interesting possibility in design. It is conceivable that for 
certain geometries we can find two fiber -reinforced shells, one possess- 
ing orthotropy, the other isotropy, such that the isotropic shell will 
buckle at a higher load theoretically. However, because of the greater 
sensitivity of the isotropic shell to initial imperfections, the ortho- 
tropic shell practically will carry a higher load and, even more important, 
may be more predictable in its behavior. 

t , 

The goal of the present research is to study the relative imperfec- 
tion sensitivity of fiber-reinforced cylinders in order to identify which 
fiber configurations are least / sensitive to initial geometric imperfec- 
tions. The investigation is accomplished by studying the behavior of a 
perfect shell in its initial postbuckling region. An index of imperfec- 
tion sensitivity as suggested by Koiter (refs. 1 6 to 2k) is determined 
by perturbing the shell displacements in the neighborhood of the bifur- 
cation point. The effects of nonlinear displacements induced by loading 
prior to buckling are considered. The theory is derived by employing 
principles of perturbation theory and represents an extension to cylin- 
drical shells of the plate postbuckling theory developed by Stein 
(ref. 25). Relative values of the sensitivity index are found for 
helically wrapped cylinders of glass-epoxy, boron-epoxy, and boron- 
aluminum loaded in axial compression. Wrapping configurations were 
varied to discover configurations of minimum imperfection sensitivity. 
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A series of compression tests was conducted on twelve glass-epoxy, 
filament -wound cylinders to investigate experimental 'buckling 'behavior. 
Theoretical imperfection sensitivity trends are compared with the 
experimental results. 



V. LIST OF SYMBOLS 


A 

^21^31 




c l> c 2 

-4 

D 


D- . 
10 


D 21^ D 31 ,E 31 


V*y 


E 


x 


W G 1 
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xy 




M x^^V' J -^xy 


cylinder length 

square matrices associated with left-hand side of 
equations ( 31 ) and (40)., respectively 
coefficient in load-displacement expansion (eq. (9)) 
imperfection sensitivity index 

extensions! stiffness of cylinder wall (see eq. ( 3 )) 
constants in homogeneous solutions (see eqs. ( 27 ) and ( 3 6 )) 
postbuckling displacement vector defined in equation (45) 
bending stiffness of cylinder wall (see eq. ( 5 )) 
column vectors associated with right-hand sides of 

secular equations (see eqs. ( 31 ), (40)., and ( 71 ) to (74)) 
equivalent Young's moduli for helically wrapped layer 
of cylinder wall 

effective Young’s modulus for contraction of composite 
cylinder wall 

constants associated with axisyrametric prebuckling 
solution for w (eq. (l8)) 

equivalent shear modulus for helically wrapped layer 
of cylinder wall 

stiffness of cylinder wall associated with coupling 
between extension and bending (see eq. ( 5 )) 
moment resultants 
stress resultants 
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number of full waves in cylinder buckle pattern in 
circumferential direction 
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a 
A 
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B *'W 

A ij 


A 


radius of cylinder measured to arbitrary reference 
surface 

total wall thickness of cylindrical shell 
displacement functions of x associated with first 
perturbation (eigenvector components) 
normalized eigenvector components 

displacements of shell wall in x- } y-, and z-direetions, 
respectively 

column vector formed from u, v, w, and evaluated 

at each finite difference station 
surface and normal coordinates for shell 
helical wrap angle 

cylinder prestress parameters (see eq. (l8)) 
shear strain in shell at reference surface 
distance between adjacent finite difference stations 
virtual work for an axially loaded cylinder 
normal strain 
curvatures 

structural constants defined in equation ( 51 ) 
magnitude of applied compressive load 1 at ends of 


cylinder 
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u ■ 




y 


R 


x 


I 

a 

T 

Superscripts: 


Poisson’s ratio for isotropic material 
equivalent Poisson's ratios for helically wrapped., 
orthotropic layer 

effective Poisson’s ratio for composite cylinder wall 
perturbation parameter, measure of amplitude of initial 
postbuckling displacements 
amplitude of initial imperfection 
normal stress 
shear stress 


i 

* 

T 


ith layer of cylinder wall 
evaluated at A = A cr 
transpose of a column vector 


Subscripts: 


0,1,2, 3 


s 

x,y 

cr 


prebuckling, first, second, or third perturbation, 

respectively (when used on load, strain, or displacement 
variables ) 

imperfect cylinder buckling 

longitudinal and circumferential directions, respectively 
at buckling 


iso isotropic ' ■ . 

exp experimental 

Primes are used to denote ordinary differentiation with -respect to 
the variable x. A subscript preceded by a comma denotes partial 
differentiation with respect to the subscript. 



VX. IMPERFECTION SENSITIVITY INDEX 


There are several studies of the sensitivity of shells to initial 
imperfections available in the literature ( see., e.g., refs. 2 6 and 27). 

The majority of these studies are restricted by the adoption of an 
initial imperfection of specified geometry. Practically speaking, 
however, real initial imperfections have a somewhat random character. 

An example of this fact is shown in figure 1. Measured radial displace- 
ments from perfect geometry are shown for a 10-foot diameter ring- 
stiffened cylindrical shell (ref. 28). The initial radial displacement 
wi divided by the shell thickness t is plotted against length for 
various circumferential stations about 6 inches apart. It is evident 
that considerable judgment and ingenuity is required to represent such 
initial imperfections analytically. Thus, in the present state of the 
art, it seems impractical to hope for an analysis which will yield the 
exact reduction in buckling load for a wide range of realistic imperfect- 

ion geometries. However, a method of assessing relative sensitivities 

/ 

f 

of shells to initial imperfections does exist. ' _ 

In his doctoral dissertation, Koiter (ref. 19.) proposed a theoretical 
method of obtaining an estimate of the sensitivity of a broad class of 
shells to initial imperfections. In the Koiter theory, it is postulated 
that the sensitivity of a structure to small initial imperfections can 
be assessed by a knowledge of the perfect structure’s behavior in the 
initial postbuckling region. In his thesis, Koiter demonstrates that 
for imperfections having the shape of the perfect shell buckling mode. 
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Inward 


J 10 


the sensitivity of the shell could he investigated with a knowledge of 
either of two fundamental parameters. 

These parameters were identified hy expanding the applied load A 
as a function of a scalar perturbation parameter £. The parameter 5 
is arbitrary and physically can be considered to be any convenient, 
suitably small measure of the growth of displacements in the shell just 
after buckling. The expansion of load in terms of this scalar measure 
of postbuckling displacement can be written as 

^ ~ \jr + a ^cr^ + ^cr^ 2 + . . . 

where A cr is the bifurcation or buckling load. Koiter has demonstrated 
that when small initial imperfections are taken into account, the 
relationship between the buckling load of the shell in the presence of 
the imperfection and the buckling load for a perfect shell is related 
to the parameters a and b. For example, Koiter has shown for a shell 
governed by "cubic" behavior (a = 0 ), the relationship between the 
buckling load A s in the presence of an initial imperfection of 
magnitude £ and the buckling load for a perfect shell A cr is given by 



The above relationship is restricted to the so-called "classical" theory 
of buckling in which the effects of deformations prior to buckling are 
ignored. Recent results (refs. ,23 and 2^) , in -which these effects are 
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considered, also conclude that b is a measure number of the reduction 
in buckling load' due to imperfections. 

Thus, depending on the type of shell, either a or b is a 
fundamental parameter in assessing relative sensitivity to small imper- 
fections, that is, the imperfection sensitivity index. Hence, the goal 
of employing Koiter’s theory to assess imperfection sensitivity is to 
understand the character of the postbuckling behavior of the structure 
in the neighborhood of a singularity, namely,; the point of buckling. 

, > t 

* i» 

This type of consideration is somewhat analogous to studies in the theory 

i * 

of nonlinear vibrations where one examines the character of singularities 
in vibrating systems (e.g., refs. 29 and 30)... In fact, -in this field 
there is even Poincare’s "index of singularity." Koiter’s work relates 
the character of the singularity occurring at buckling (bifurcation) to 
the structure’s sensitivity to initial imperfections. 



VII. DERIVATION OF GOVERNING EQUATIONS 


Assumptions 

In developing the theoretical estimate of imperfection sensitivity 
of the cylinder , several basic assumptions are made. The cylinder is 
assumed to he composed of several orthotropic layers, each with one of 
its principal axes of orthotropy alined with the cylinder axis . The 
effects of transverse shear are ignored. At the onset of loading, the 
cylinder is assumed to be perfectly cylindrical so that its shape is 
changed only by deformations induced by load. The cylinder is assumed 
to buckle into a unique buckling mode, so that the behavior of the shell 
can be investigated in the neighborhood of a bifurcation point. Eor a 
cylinder, this restriction precludes consideration of cylinders that 
"buckle" in an axisymmetric mode. Axi symmetric "buckling" in a cylinder 
is simply unchecked growth of the axisymmetric ’deformations' induced by 
loading and is often described as collapse rather than buckling. Since 
no bifurcation occurs, the present theory cannot assess relative 
imperfection sensitivity for this phenomenon. 

The theory is developed by adopting the Koiter parameters to obtain 
the sensitivity index of a perfect cylinder. The development of the 
theory differs from that existing in the literature (refs. 1 6 to 24-) in 
that the index is identified by employing perturbation techniques similar 
to those of reference 25. Nonlinear equilibrium equations formulated 
from Donnell-von Karman strain displacement relations are perturbed about 
the bifurcation point to find equations governing the behavior of the 


12 



15 

shell in the initial postbuekling region. The perturbation is shown to 
produce secular terms which are constrained to vanish. The constraint 
conditions yield relationships between the sensitivity index and the 
buckling mode. The sensitivity index is then computed by using the 
constraint condition together with a suitable definition of the 
perturbation parameter. 


Equilibrium Equations 

Nonlinear equilibrium equations are formulated using the Donnell - 
von Karman strain displacement relations. Under the assumptions, the 
virtual work expression for a multilayered orthotropic cylinder loaded 
in axial compression can be written as follows: 


n A o2kR 

51 = J 0 J 0 ( V e * + V s y - W + + w + W 3 * ay 


r>2j£R 

J 0 
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dy = 0 


( 1 ) 
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In equations (2), Cj^,. D^j, and ICy are structural stiffnesses of 

the multilayered wall associated with extension, bending and coupling 
between bending and extension, respectively. If the stress-strain 
relations for the ith cylinder layer are written as- 
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' '4 


-\ 


1 - U.- 1 -!!?- 
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E . I. 
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(4) 


then, the structural stiffnesses can be defined by the following 


integrals: 
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The integrals appearing in equations ( 5 ) can be evaluated for an arbitrary reference surface by 
using equations (83) and (85) of reference 31. The elastic constants for fiber-reinforced structures 

\JT 

appearing in equations ( 4 ) can be predicted -analytically from the constituent properties of the 
fiber and matrix. The constants for a layer in ; which the fibers are unidirectional can be obtained 
from any of references 32 "to 35 . For layers with fibers oriented at angles (see fig. 2 ), 

an "equivalent" orthotropie layer is defined through compatibility considerations and the trans- 
formation equations of orthotropic elasticity (see ref s._ 35 to 37). The formulae employed herein 
are summarized in appendix B of reference 1 . 

Equation (l) can be expressed as a function of the .displacement variables u, v, and w. 
Integration by parts yields the following equilibrium equations which are valid in the postbuckling 


region : 




Figure 2.- Geometry of fiber-reinforced cylinders. 
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1L. + 1L-,, = 0 
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" **>c,xx ' m xy,xy " “y^yy + " ,1T x w , 


- 2N w •- N„w _ = 0 
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The boundary conditions re suiting. from the variations* with the load 
specified and the- cylinder supported at the ends, can be written as 


N x + A = 0.’ 
* \ 




= 0 or- y = 0 


■ 


Mx- = 0 or w. x =. 0 


w = Q 


( 7 ) 


J 


Perturbation Procedure 

Following the Koiter procedure-, the displacements are expanded as 


u = i^(£) + Uj^i + Ugi 2 + . . . 

v = v 0 (|) +■ v^l + v 2 | 2 +....) 

W = Wq ( | ) + w,_£ + wg£ 2 + .... ^ 


(8) 


where the subscript 0 -denotes- the axisymmetric state of deformation 
induced by load -and £ is the perturbation parameter. 

-The applied compressive load A is expanded as 

A = A cr + aA cr l + lA cr £^ + . . . . 


( 9 ) 
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where A cr is the applied load at "buckling. At this point , it .should * 
be observed that Koi'ter ' s work was based on -classical- buckling analysis 
so that the dependence of the axi symmetric deformatiqns (uq, v 0 , w 0 ) on 
I i-s not clearly noted. Fitch (ref. 2k) has recently generalized 
Koiter’s work to include the effects of deformations prior to buckling. 
In Fitch’s work, the dependence of the axisymmetric deformations on | 
can be found by expanding these deformations in a Taylor series about 
the buckling load A cr - Thus, 


u 0 = u* + (A - A cr )u*^ + |(A - A cr ) 2 u*^ + . . . . 

v 0 = v* + (A - A cr )v*^ + |(A - A cr ) 2 v*^ + • • • X 

w 0 = w* + (A - A cr )w*^ + §(* - A cr) 2 + • • * ■ 


( 10 ) 


where the * denotes quantities evaluated at A = A cr . For example. 


w, 


.* 


8^w 


°,xA 5x3 A 


A = A 


cr 


Using equations ( 9 ) and (10) 


.u = u* + (aA u* + u )£ + (— a 2 A 2 u* + bA u* + u j£ 2 + . . . 
0 cr 0^ l' \2 cr 0 ^ cr 0^ 2 ) 




v = v* + (aA v* + v_ ) ^ + ( - a 2 A 2 v* + bA v* + v ]g 2 + .... I 
0 crO jA l' l 2 C **°,AA cr "°,A 2 ? 


w = w* + (aA w* + W- ) ^ + ( ^ a 2 A 2 w* + bA w* + w_ ]| 2 +' , 
0 cr °,A 1 \ 2 cr °,AA cr °,A , 2 ) 


. ,( 11 > 



19 


Zeroth Perturbation 

If equations (ll) are substituted into the equilibrium equations (6) 
and the boundary conditions (7)> a series of equations corresponding to 
powers of £ are obtained. The equations corresponding to the zeroth 
power of % can be written as follows 
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where the prime denotes ordinary differentiation' with respect to x. 

The stress resultants of equations (l2) and (13) are defined hy substi- 
tuting equations (14-) into the constitutive relations of equations (?) . 

Equations (12) and (13) can. be identified as equations determining, 
the axi symmetric deformations' occuring in. the cylinder prior to buckling. 
The asterisk indicates that the deformations resulting from these 
equations must be evaluated at the inception of buckling, that is, at 
A = Acr • 

Solutions to equations (12) can be obtained by integrating 
equation (12a) and applying boundary condition ( 13 a) so that 


N* = - A. 

•“■ft 


cr 


(15) 


From the definition of , it follows that 

x o 




= u*' + ifw** 12 „ Si ' -212 
0 ^ 0 C U 0 C 1L H 


Zisr 

C 11 


(16) 


Integration of equation (12b) together with either of boundary conditions 
(ljb) imply that both v* and N^y are zero everywhere. Substitution 
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of equation (l6) into equation (12c) yields the following ordinary 
differential equation in w*. 



ikV*'”'* 

Oil/ 0 


— / 

2/ K ll C 12 

_R\ C n 





1 ^12 
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(IT) 


For simply supported shells (Mi = 0) with houndary conditions applied 

o 

at x = ± A/2 when A' is the cylinder length, the solution of (17) can 
he written as 

w* = F si'nh 0x sin yx + F cosh px cos yx + ?\ P (l8,) 

u x u • cr x 

where 
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P l = 


c 12 r 


C 11 C 22 " C 12 


F, = 


2 Pll 


G, -cosh ^ cos 2 ^ _ g sink ^ sin- ^ 
1 2 2 2 2 2 


^^P 7 ^sinh 2 ~ sin 2 ^ + cosh 2 cos 2 ^ 


F„ = 


Gp cosh — cos — + Gn sinh — sin — 

2 2 2 1 2 2 


k!i n 


2 .ID-, n — )• P 7 [ sinh 2 — sin 2 — + cosh 2 — cos 2 — • 

V C U / \ 2 2 2 2 j 


with 


Gl - A, 


cr 


K H + R C 12 / d 11 C 21 - K 12 V p2 7 2J 


"11 C]L1 ^ C ll C 22 " C l£ 


G 2 = 2 P 7 R 


C 12'[ I> Upll “ 


■12 


J ll\C| L iC 2 2 “ c | 2 j 


For clamped boundary conditions, the radial displacement w* is given 
by (l 8 ) with 


A 


. . , BA 7 A , PA . 7A 

P sinh cos — 7 cosh sin ~ 


F 1 " ^cr 2 ! 


. . 7A 7A , . .. PA , PA 

p sin — cos — + 7 sinh — cosh — 
2 2 2 ' 2 ' 
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f 2 = 


7 sinh cos ^ + p cosh sin ^ 
2 2 2 2 


„ . 7 A 7 A , . , PA , pA 

p sin — cos — + 7 sinh — cos±l — 

2 2 2 2 


Note, in general, that quantities such as w Q ^ can be found from 
equations (l8) and ( 19 ) by replacing A cr by A, performing the 
indicated partial differentiation and evaluating the result at A = A cr . 

hi 

For the quantity ,xK’ "hhe result can be written as 

w o,xA ~ cosl1 3 X sin 7x(FqP - + sin k P x cos 7 x(Fi 7 + F 2 p)^ 

+ x sinh px sin 7X P^(-FiP - FgA) - A^A^lA + ^2(3) 

— ■ 

+ x cosh px-cos 7X p -^(F^A + FgP) t 7 -^(F^P - Fg 7 ) (20' 


with ‘ 


7 ' V Win Ci, -K?,> 


, , *> $ 

\ 


, • ( - Qn 

P v = - — — 

8p(b-| -[ C-] ■] - K^,)= 

* t 

For simple support boundary conditions, if 


Ft = 


and F 2 = 



F, , . cosh B oos ^ sinlx M sin £ 

1,A G- 2 2 G 2 2 


+ A sinh M cos A^P^ - G 2 7 , A ) - A cosh A sin + M,*) 

- ^ 7 p) isinh 2 A sin 2 A + cosh 2 A cos 2 2*) 
g 2 2 2 2 2/ 


~T 73 a/p^ sinh A cosh H - 7^ sin A C os A 


~ — 7 -> sin — - 
2 >A 2 


F, - , - slnh ££ sin 2£ - ^ oosh 

G 2 2 G 22 


+ A cosh A sin ^W, A - G 1 P,a) ' gg sintl T co ^ T< G 1 7 ,A + G 2 P,a) 


pA 7A, 


+ — A (7P) A /sinh 2 A sin 2 A + cosh 2 A cos 2 A 
G^ } \ c. ^ 


2Ggi ( RA pA j 7A 7A 

+ —A 7PAf p ;A Sinh A cosh A - 7,7, sin ^ 008 “ 




with 


*U 


°l,3v = P i[p 2 - * + 2 ^r(W,X - yy.t) - - 2 i 

L V°ii ‘ *ii J 


+ W 

w 


For clamped boundary conditions, if 


*1 " 


^cr p l G 12 

a 


and F 2 = - 


G 


F- 


-> = + ^ * -2:1 $ , sinh cos — - 7 \ cosh ^ sin — ^ 

G G \ > A 2 2 > A 2 2 j 


. A ^crPl 

T — =— 




2 G 

Acr^lG lg 
2G £ 


cosh M cos 2^(pp -v - 77 -v) - sinh H sin 2£(p7) ^ 

2 2 ■ >A 2 2 jA 


<P^ sin /A + AP7^ cos 7A + 7^ sinh pA + A7P^ cosh pA) 



( 22 ) 


^2,A = ' c ° sh sin + 7p sintl dos ^ 


A ^cr^l 


2 -G 


sinii J~£ sin ~(p0^ - 77,?\) + cosh M cos ^{p7),A 


A cr p l G 21 


+ — — (p ^ sin 7A + A37 ^ cos. 7A + 7 ^ sirih 3A + A70 cosh f3A) J 

2 G . ’ ■ s' 


It should be noted that equations (l8) through ( 22) were developed by exploiting symmetry of the 
boundary conditions using a coordinate x with origin at. the middle of the cylinder. For subsequent 

_v v * . ' „ ... 

developments , a more convenient coordinate is ,x whose '-origin is at the end of the cylinder. 
Equations (l8) through (20) may be written for this coordinate by using the transformation 
x = x - A/2. - - 


ro 

o\ 


First Perturbation 1 _ 

The equations corresponding to the first power of £ can be written as follows: 


N + N - =0 

xi,x Vi,y 


N xy x ,x + N y lf y = 0 




N- 


~ *" yi i ■X" 

“ ^x^.xx “ 2 ^5cy 1 ,xy “ M yq,yy + “ “ Nx o Wl ^ “ Nx i w 9 #*x " w y 0 w l/yy 


= 0 


( 23 ) 
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with boundary conditions 


%! = 0 


V 1 = 0 or lT xyi = 0 

= o 


V = 0 or w l,x = 


= 0 




In equations (23) .and (24 ) and subsequent perturbations 


(24) 


e Xj_ ” u i,x + w o,x w i,x 


v i 

e = v* + — 

yi v i,y R 




* 


7 xy, = u i,y + v i,x + w o,x w i,y 


*X4 — w i,xx 


/ i = 1,2,3, 


(25) 


tt yi " " w i,yy 


rt^y^ = - 2^^^ 


and the corresponding stress resultants are obtained by substituting 
equations (25) into the •• constitutive relations (2) 

To obtain equations (23')., the .equilibrium equations resulting 

f -- 1 
' “ ‘ I ' 

from the zeroth .perturbation (eqs.. ( 12 )) y we're used, to eliminate terms 

* « , ^ 1 H ^ 

* ‘ f I t ‘ ^ ^ 

such as aA cr l , I x obtained from the perturbation. * Also in obtaining 
° * •’ ’ r -“ 
the boundary conditions .exploitation was made of. the fact that terms 
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having the functional form aA f* , vanish at the boundary if f* is 

cr o,A o 

zero there. Equations (23) and (24) are similar to the usual equations 
governing buckling of a pre stressed cylinder except for the absence of 
a free buckling parameter. The buckling parameter is absent since, in 
the present problem, the displacauents of the shell have been expanded 
about a single, known point of bifurcation. Thus, as opposed to the 
usual buckling problem in which a multitude of eigenvectors are possible 
solutions, only one eigenvector is an admissible solution of (23) -namely 
that corresponding to = - A cr - If both the buckling load A cr 

as well as the number of circumferential waves n in the buckle pattern 
are established by a conventional buckling analysis (e.g., ref. 38) 
equations (23) and (24) can be converted into a set of ordinary 
simultaneous differential equations and boundary conditions by the 
following substitution: 

U 1 = U L (x) 

V 1 = 7 l( x ) 

w-j_ = Wp(x) 

As was done in reference the eigenvector corresponding to = - A cr 

can be found by the method of finite differences and determinant plotting. 
The components of the resulting eigenvector can be expressed as 


cos 


sxn 


ny 

R 


s 

R 




cos 


ny 

R 


J 


(2 6) 
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U 1 =vc l u l( x ) cos ? 




E 


V 1 - c i v i( x ) sin Y f 


( 27 ) 


wi = CiWxCx) cos 


where cp is an arbitrary constant and U]_, Vi, and W]_ are suitably 
normalized components. 

Second Perturbation 

The equations corresponding to the second power of £ can be 
written as follows: 


U + N 
x 2 ,x xy 2 ,x 


= - aA. 


‘cr 


c ll( w o,xA w l,x> + c 66 v o,xA w l,yy 


L 


§( c ll 4,x + c 12 W 1 ,y)^ x “ c 66( w ljX w l ; y^y. ^ 28a ) 


K xy 2 ,x + N y 2 ,y = “ a ^cr 




c 66( w S,xA w l,y> x + c 12 w o,xA w l,xy 


L 


' |(°12 w l,x + C 22 w l,y)^ y “ C 66( w l,x w l,y) 


( 28 b) 



%p 

- ~ + "r“ " N x 0 “ Nx 2 w 0,xx " w y 0 W2 ^ 

= a ^cr w $_,xxA + N Xq, A w l,xx + ^yQ,A w l,yy + C ll w ^xA w $,xx w l,x " r w $,xA w 1,x 

+ K n( w 5,xA' w l,x) )33c - - K 12 w S,xA ’•l.JW + ^66H,xK v l,yO + ' J (K 11 w l,x 

+ K 12 W^ y )^ + |(K 12 w^ k + K 22 w^y)^ + 2Kg6( w l,x w l ; y^ X y 
- ~( c 12 n x + c 22 w l,y) +-“(Cn wf,* + C 12 w^y^xx 

cit\ * 

+ N Xl w i ; xx + 2N xy-j_ w l,xy + %•]_ ¥ l,yy 


(28c) 
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with boundary conditions. 

K x 2 + | C 11 w l,x + | C 12 w l,y = 0 

N xy 2 + C 66 w l,x w l,y “ 0 or v 2 = 0 
w g = 0 

^x 2 + | K 11 w l,x + | k 12 w l,y = 0 or w 2 ,x 

Equations (28) and (29) can be converted into a set of ordinary 
simultaneous differential equations and boundary conditions by the 
following substitutions 


(29) 


= 0 


2ny 


u 2 = U 20 (x) + Ugit 251 ) cos ~ + U 22 (x) cos 


Vp = Vpn (x) sin — + Vpp(x) sin- 
c cu. R R 


(30) 


W 2 “ W 20 (X) + W 21< x) C0S Y + W 22^ COS IT 


J 


Since the differential equations of (28) are linear, solutions 
corresponding to each set of Ug^, V g^ , and (i = 0,1,2) can be 

superposed. However, if the particular solution for Ug-^, Vp-] , and W g-j 
is sought by algebraic means from equations (28) and (29), certain 
difficulties are encountered. If equations (28) and a set of boundary 
conditions from equation (29) are cast in finite-difference form after 
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substitution of (30 ) , the resulting simultaneous algebraic equations 
can be written in matrix notation as 

= a *cr c l D 21 ^l) 

where Zg-^ a column vector whose components .are Ug-j_, Vgj , Wgp 

4 

and evaluated at each difference station and is a column 

vector containing elements corresponding to the terms multiplied by a 
on the right-hand side of equations (28) and (29). Because the 
trigonometric function associated with Ug^, ^ 21 ’ 311(1 W 21 laas 

argument ny/R, the determinant of Agj_ has a form identical to that 
used in determining the buckling load A cr . Hence, 

det Aqi = 0 (32) 

which demonstrates that the particular solution resulting from (31) is 
a secular term which must be removed from the perturbation. A condition 
for removal of the secular term is the requirement that 

A 21 Z 21 = a ^cr D 21 = 0 ^3) 

If Z21 is identified as the eigenvector Zi corresponding to the 
solution of equations (2J) and (24) , then the desired condition for 
removal of the secular terms is 


A l z l “ A 21 Z 21 = 0 " aA cr d 21 


or 


Z 1 A 21 Z 1 " a7v cr c l Z l D 21 


= 0 


(34) 
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T T 

where Zj_ is the transpose of Zj_. Since Zj_ D 2 ;j_ is finite and 7\ cr and 

cj_ must he finite for buckling to occur, it follows that 

a = 0 (35) 

Thus the only admissible solution corresponding to the functions U , 

V p-, , and i s the solution to the homogeneous system 

A 21 Z 21 = 0 

Hence, the general solution to equations (28) and ( 29 ) with a = 0 
becomes 


u 2 = U 2 q(x). + CgU-^x) cos — + Upp (x) cos 


v 2 = c^ (x) sin ^ + V 22 (x) sin ^ 


( 36 ) 


w 2 =? W 2Q (x) + c 2 W 1 (x) cos ^ + W 22 (x) cos 


J 2"l v 


R 


R 


J 


where c 2 is an arbitrary constant and Uj_, V^_, and W-j_ are the 
normalized components of the buckling eigenvector (see eq. ( 27 )). 

» 

Third Perturbation 

If a = 0, the equations corresponding to the third power of £ 


can be written as follows 





IT + N = - bA 
Xj,x xy^,y cr 


c ll( ¥ o,xA w i,x^ x + c 66 w o,xA v i,yy 


A 


C ll( w l,x v 2 ,x) x + Cl2( v l,y ¥ 2,y) x + %6^ l,x w 2,y + ¥ l,y w 2,x) 


*y 


I + u = - bA 
“xy^x "y^y cr 


c 66 ( w 0 ,xA v l,y) jX + °12 ” 0 ,xA *l,xyj 


C 66< W 1,X v 2,y + »l,y V 2,x\ t _ + C X2K,x W 2,x> )y + °22K,y w 2,y)^ 


L 


Ny* 

" 2H xyj,xy “ ®V 5 ,yy + E " W x Q %xx “ n x^ w O,xx ” W y Q w 3,yy 


= bA, 


cr 


+ K 11^ W 0^xA W 1,x\xx + 2B €6^ w O,xX W l,y^ xy + ^K^xA w l,x^yy + K 01^ w l,x W 2,x^ 
K,y "2,7^ + < v l,l w 2,x> >yy ] + K 22 ( Wl ,y Vs,y) jJy + 2K SS<''l,x »E,y 


,XX 


+ K- 


•j 

+ w l,y v 2,x) #x y + ( C ll w l,x w 2,x + C 12 w l,y w 2,y) w 0 ; xx “ *( C 12 w l,x w 2, 


x 


\J 4 

P=~ 


®Xq^A ¥ 1 >xx + \ w $,xxA + ^y^A w ljyy + ^il w $,xx w $,xA ¥ l,x ¥ $,xA ¥ l,x 1 (57) 



+ c 2a -w 1)y v 2>y ) + ^ jX + c 12 + Sf5L(o 12 vf >3[ 

+ C 22 w 1 ^y) + 2Cgg W 1;X W-L^y Wi ;X y + 1 ^ 1 w 2,xx + w l,:*x 
+ ^ w 2,yy + N y 2 w l,yy + 2N xy 1 w 2,xy + ^xyg w l,xy 

with boundary conditions 

+ C 11 w l,x w 2,x + C 12 w l,y w 2,y + ^cr 0 !! %xA w l,x = 0 
W xy 3 + c 66( w l,x w 2,y + w l,y w 2,x) + b *cr c S6 w 0,xA v l,y = 0 or v 3 

w^ — 0 

Mx 3 + &□_ w 1 ^ x w 2 ^ x + K 12 w-^y Wg^y + tA cr K lx w 0 ^ xA V lf3t * 0 or 
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Equations (37) and (38) can be converted into a set of ordinary- 
simultaneous differential equations and boundary conditions by the 
following substitutions: 




u 3 = U 30 (x) * u 3l( x ) cos ^ + U 32 (x) cos QjL + U 33 (x) cos 222 
v 3 = V 33 (x) sin 22 + V 3 2 (x) sin ^52 + V 33 sin 2SZ 


} (39) 


w 3 = ^30^ + w 3l( x ) cos ““ + W 32 (x) cos ~~ + ^3 3 (x) cos -221 


R 


J 


Paralleling the second perturbation,. U-^, V^, and can be 

shown to be secular terms which must be eliminated. The resulting 

A j 

algebraic equation in this case can be written as: 


A 31 Z 31 “ bA cr c l D 31 + C 1 E 31 


m 


where Djq is a column vector corresponding to the terms multiplied 
by b on the right-hand sides of (37) and (38) and E31 is a column 
matrix generated from the terms that are independent of b. In 
contrast to equation (31) > the additional term multiplied by 
occurs in (4o) .because, when the general solutions (eqs. (27) and (36)) 
are substituted into the right-hand sides of (37) and (38), terms such 
as cos ny/R cos 2ny/R re stilt which can be expanded as 


cos 


ny 

R 


cos 


2ny 

R 


l/ 3ny , 

— cos — — ^ + cos 
2\ R 
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thereby yielding additional secular terms that are independent of b. 
The condition for removal of the secular terms that is the analogue of 
equation (3*0 can be written as 


7,? A 7,. a bTv 


_T 


which has the solution 


(4l) 


b 


~ c l z i E 3i 

Vr 2 ! D 31 


(42) 


Equation (42)' together with consideration of the equilibrium equations 
and boundary conditions corresponding to the first, second and third 
perturbations can be shown to imply the following proportionalities: 


w l> N xy x > %i * ' t,1//2 ' 

w 2> V ^xy 2 > ^y 2 “ 13 ^ 

W 5 , ff^, % 5 « b3/2 , 


(*»3) 


It is interesting to note that the counterpart of equation (42) was found 
by Stein in a study of the postbuckling behavior of plates. (See 
equation (19) of reference 25 *) 

The proportionalities indicated in equations (43) have also been 
found in existing applications of the Koiter theory. Eor example, Fitch 
and Cohen in references 23 and 24 have obtained an explicit expression 
for b from energy considerations. In terms of the present notation, 
the results of references 23 and 24 can be expressed as; 
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If 


b = - 


A is dy + 2 / / B dx, dy 

surface surface 




cr 


// 


(¥ 4 ) 


C dx dy 


surface 


where 


A " ^Xg w l,x + N y 2 W l>y + ^xyg Wl > x + 2 Cl1 Wl ,x + C d- 2 W l,x w l,y 


+ g ^22 w l,y + 2 ^gg w^y 


B = w l,x w 2,x + % x w l,y w 2,y + ^(^x^y + w X , y w 2 , x ) 


C = N* -. w? + N* w? + 2(N w* , w, + N„, w£ - w T ,,) 
Xq,A 1>x yQ,A l*y OjxA ~,x xy^ l*y 


Using the proportionalities indicated in equation (43), it can he seen . 

i 

that equation ( 44 ) becomes independent of b and hence equation ( 44 ) 
must be the counterpart of (42). The general theory of references 23 
and 24 can be specialized to a "classical" buckling theory by taking pre- 
buckling deformations consistent with a membrane state of stress. The 
index b has been obtained under these assumptions by Budiansky in 
references 20 to 22. For the cylinder, the use of classical buckling 
theory requires derivatives of w* to vanish. Inspection of equation ( 44 ) 
with this assumption reveals again that the equation is independent of b. 
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Determination of t 

To determine b explicitly, the perturbation parameter | must 
be defined appropriately. If the parameter £ is considered to be a 
scalar measure of, the growth of the displacements in the shell just after 
buckling, then £ must be related to the magnitude of these displacements. 
The displacement vector D is defined to be the following: 

T) = (u - u*)? + v~f + (w - w*)fef (4-5) 

so that to the second approximation using (ll) with a = 0 


— > 

D = 


hiE + (bA cr Uo^ + u 2 ) g 


i + 


( T l& 


+ V 2 £ 2 )j + 


wi£ + (bA cr wJ jA 


+ w P ) i c 


£ 


(k6) 


The parameter £ in the present theory is taken to be such that the 
square of the perturbation parameter is equal to the square of the 
magnitude of the vector S’, (if • D)., averaged over the cylinder surface 
and nondimensionalized by the square of the cylinder thickness t.. Hence 


6 2 = 



D • D dx dy 


surface 


'2*RAt 2 


(4-7) 


Using equation (46) 


2jtRAt 2 £ 2 = £ 2 jT jT (u 2 + v 2 + w 2 )dx dy + 0(g3) 


m 


surface 



Using the general solution for the first perturbation (eq. (27)) and 
expressing in terms of b from equation (42 ) , b can be written as 


Thus, 


b = 


- 2jcRAt 2 Z^E 5] 


\ 7? 


(49) 


to within o(g 2 ), the expression for b becomes 


b = 


- 2At 2 Z^E 51 


r 

D 3i J + vf + wf)dx 


0 


(50) 



VIII. NUMERICAL SOLUTION OF POSTBUCKLING EQUATIONS 


The solution for the state of axi symmetric prestress at "buckling 
has already "been presented in closed form (eqs. (l8) and ( 19 ))* The 
solution to the equilibrium, equations and boundary conditions resulting 
from the first and second perturbations (eqs. (23), ( 24 ), (28), and (29)) 
were found for simply supported and clamped shells by using numerical 
methods. Following the finite -difference technique suggested by refer- 
ences 38 and 39, the system of ordinary differential equations resulting 
from substitution of the general solutions (eqs. (27) and (36)) axe cast 
into a set of simultaneous second-order differential equations and 
boundary conditions involving the fundamental displacement variables 
U, V, and W and an auxiliary variable M x - By employing central 
differences, the equations are converted into simultaneous algebraic 
equations which are- then solved on a digital computer by matrix methods. 
Because the resulting algebraic equations contain banded matrices, it is 
expedient to employ Gaussian elimination in the solution. 

General Solution 

In formulating the postbuckling solution, it is convenient to con- 
sider the left-hand sides of the equilibrium equations in general. When 
the general solutions for the first and second perturbations are substi- 
tuted into the equilibrium equations and the definition of the moment 
resultant M x is introduced, the following set of second-order differen- 
tial equations and boundary conditions are obtained: 


in 
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A 11 TJ " + A 12 U + A 13 v ' + A 14 w " + %5 W ' + A l6 W + A 1T M ’ = D 1 

AgjU’ + ^ 22 ^" + Ag^W 1 + A^jgW = 

a 51 U' + A 32 V" + A 35 V + A^W" + A 55 W’ +, A 3 gW - M" = D,. 


A 41 U ' + A 42 V + a 43 w " + A 44 w ' + a 4 5 w - + M = \ j 


where 


A 11 = c .ll ~ 


K‘ 


D 


-11 

ii 


A 12 " 




a 13 " 


L , „ K 12 K ll\mn 

l Ci2 + 066 ■ -^r)* 


A n 1, = C 


* 

W r 


14 - -ii 577ro, x 


A 1 5 - 


11 


1: 

4i 

D l: 


j) W 0,*x + (' 


KlgKjjNl 

12 " D 11 h 


- + K 12 + 2K 66 


k 11 d 12 


D 


11 


a - r> * /mn\2 

A l6 ~ ~ C 66™0 } xi R J 


K- 


A. 


17 D. 


11 

ii 


^1 “ ~( C 12 + C 66^ Y 
A 22 = c 66 



43 

A 23 ■ -°22 

- < K 12 + 2I W Y 

Ag 5 = -(c 12 + C 66 )w^ x ^ 

A _ p * am C 2^ rsn „ /mn\3 

^6 - -C 66 w 0 ,xx “a “ -^ 2 — " MtJ 

A 31 = IT + (K 12 + 2K 66)(f) 2 - c H w O,xx 

A 32 = ~ 2K 66 f 

C 22 mn /mn\5 p * mn 

A 33 " ~ ^2 + K ^2(irj - C 12 w 0,xx X 

A 34 “ \r - ~ “ (4Dg 6 + D 12 )(y) + K H w O,xx 


a _ C 12 , (v . ov- \/mn\2 * _ * * 

A 35 [Jr - + K 12 + 2K 66 J'Yj w 0,x ■ C ll w O,x v O,xx 

A 35 = + z <*66 - K 12](f) 2 - Irj^xx + l( C22 - 

-SfK 


o* 




U 20> V 20 J W 20> M 20 
U l, V x , W x , M x 

^ 22 ^ ^ 22 J ^ 22 ’ ^22 

and (i = 1 to 4) the corresponding loading terras. 

If the central difference formulae 




are employed, equations ( 51 ) can be written at the ith station as 


+ Bj_Zj_ + Cj_Z$_+2_ = D'i' (i = 1,2,-. ..,N) 


(53) 



**•5 


■with 


B 


where A 
so that 


A 11 

" 1 A 13 

A 14 " 2 ^5 • 

. - A A-, 7 
2 17 


^22 

^4 ~ § ^5 

0 

- § a 31 

A 32 

V " 1 A 35 

-1 

j |4l 

0 

a 43 “ f A 44 

0 


-2A 11 + a2a I2 0 ’ 

0 -2A 22 + a2a 23 

0 -2Aj2 + A 2 A^ 3 

0 A \ 2 


-2A 14 + a2a 16 0 

_2A 2l{- + a2a 2 6 0 

-2A 54 + A 2 A ? g 2 

-2A 43 + A\ 5 A 2 


- A 11 

rol t> 

<*l 

a 14 + f A 15 

A A 

2 17 

A 

2 <21 

A22 

^4 + § ^5 

0 

A A 

iV 

A 32 


-1 

f \i 

0 

a 43 + f A 44 

0 


is the width of a finite -difference space and i is numbered 
i - 1, W corresponds to the boundaries of the cylinder. 
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Boundary Conditions 

For the present investigation, the so-called "classical" simple 
support and clamped boundary conditions were employed, that is, for simple 
support 


W x = " 

v = 0 

> 

w = 0 

at the boundaries, while for clamped support 

N x = -A " 
T - 0 
w = 0 I 


( 54 ) 


( 55 ) 


at the boundaries. 

Using the results of the first and second perturbations (eqs. (27) 
and ( 36 )) and the finite -difference approximations, the boundary conditions 
can be written as 


A 1 Z 0 + B 1 Z 1 + ' C 1 Z 2 = D 1 i 

A K%-1 + %% + %%+l = % 


(56) 


where for simple support 
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and. Ajy, and are given by similar expressions differing only 

in that Wq x is evaluated at station N rather than station 1. 
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•For clamped support, the last row of equations (56) is modified 
so that 


) — 

0 

•0 

-1 0 

*v 

II 

•\ 

1 pT* 


0 

— 

0 .0 


) - 

tv 

0 

f °_^ 

= ) 

/ 

J 


The contribution by points exterior to the boundary {denoted by the 
subscripts 0 and N+l) introduced by the finite-difference scheme is 
eliminated by writing the equilibrium equation at the boundary (see, for 
example, ref. 38), thereby obtaining two simultaneous equations in three 
unknowns at each boundary. If the exterior points are algebraically 
eliminated, the equilibrium equations and boundary conditions take the 
following form: 
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( 58 ) 
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and. 




^(i) 

Vi(i) 


Wi(i) 


where 


Mi(i) 


i denotes the ith station 
i = (1,2 ,..., N) 


Since is a handed matrix, it is convenient to triangul arize a i:l 

"by Gaussian elimination so that equation (59) becomes 


A n z n “ 0 


(60) 


where 




R 1 C 1 


R 2 C 2 


B, 


: W-1 %-l 


N 


and 


with 


\ = B 1 


R i = B i " A i P i-l i = 2,3 , • • • ,N 


A “] A 

P 1 = B 1 C 1 




p i = R i c i i = 2 , 3 ,...,w 


(61) 
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where 


®1 - B 1 “ A-lA-j^B-l 


c-l = c 1 . - a 1 a£ 1 c 1 


% = % - a A \ 


% = % ' % c jst % 


% % - Visr % 


®N “ % “ C B% 3 % 


Solutions to equations (58) for each of the perturbations are discussed 
in the following sections. 


Solution for First Perturbation 

For the first perturbation, the right-hand side of equation (58) 
becomes a null vector so that the equation becomes 

■ 0 (59) 


where 
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B 1 C 1 




Ag B 2 C 2 
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with m = 1 
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A condition for solution of equation (60) is that 

A 

det A- | - | = 0 


so that 


det(R 1 Rg. . .Ejj) = 0 

Using Potters' method (ref. 42) applied to a general eigenvalue problem, 
it can be shown that det = 0 is a sufficient condition for solution. 
This determinant, however, contains spurious poles which can be removed 
conveniently by the method outlined in reference 40 using the condition 


det R-^ det Rp 

- . . . det % = 0 

det Rq det R2 


(62) 


For the first perturbation solution, A cr must be such that 
equation (62) is satisfied. For the present study, the buckling load was 
determined from a computer program developed in reference 41 by using the 
techniques of reference 38, and equation (62) was used to verify the , . 
accuracy of the present equations. The buckling solution of reference 41 
differs from the present solution for the zeroth and first perturbation 
only in that deformations prior to buckling are- found from a numerical 
rather than closed-form solution, and that provision is made to- investigate 
various values of n and the applied load in order to find the buckling 
load A cr by determinant plotting (see also ref. 38 ). Once the buckling 
load -has been determined, its. corresponding eigenvector can be found by 
a method suggested by reference 38. 



52 


■ If 


a ll 

a 12 

a 13 

a l4 

®21 

^22 

^3 

^2b 

a 31 

a 32 

a 33 

a 34 

a 41 

%2 

a 43 

a 44 


it can be shown that the components of the eigenvector at the boundary 


are given by 


U-^N) = det 
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(63) 
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With the components of the vector Zj]_(W ) known, the eigenvectors at 
other stations are found, from the recurrence relation 


Z i;L (i)= +1) i = 1,2,..., N-l (64) 


After obtaining Z at all stations, the resulting vector was normalized 
for the present computations so that 


with 


z ll = 


Z 11 

t 


(65) 


z ll^^ " \ 


rv 1 '' 
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»x(i) 


i) 




i = 1,2, ,N 


This normalisation was used as a convenience to compare results from 
conventional Koiter theory (e.g., eq. (44)) with the present theory. 


Solution for Second 'Perturbation 
For the second perturbation, the equilibrium equations can be 
written as 


A 2m Z 2m “ D 2m 


m = 0,2 


( 66 ) 
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with' 
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where m takes on the value 0 to evaluate the axi symmetric part of the 
general solution (eq_. (36)), for example,. W 2 q 7 and m has the value 2 
for the asymmetric part, for example, W 22 - In equation (66), the compo- 
nents of the vector D^, can he written at the ith station as 
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where the bar over the stress resultants indicates the part of the stress 
resultant that is independent of y (viz., N x ^ = N x cos ~). 

Following the development used to eliminate exterior points in the 
general solution (see section entitled : Boundary Conditions) 


W 1 ) - W 1 ’ - Vi\w" 

> 

= D 2m( M > - j 


(68) 


where A^, A^, C-^, and C-|_ are evaluated for m =• 0 or m = 2. 


For simple support boundary conditions 


r 




“an - -* a < 


0 


i = 1,N 


( 69 ) 


and for clamped boundary conditions, Dg m becomes a null vector. It 

should also be noted that rigid body displacements will occur in U 20 

if KL is prescribed at both ends of the ’ cylinder . To suppress these 
*20 

deformations, the boundary conditions of equations (68) and (69) were 
modified so that the axial displacement was fixed at one end (i.e.. 


u 20 (n) = 0). 

Thus, the solution for displacements resulting from the second 
perturbation requires solutions to two sets of simultaneous algebraic 
equations having the form 
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Bi c x 


a 2 b 2 c 2 


a 3 b 5 C 3 


Z 2m( 2 ) 

Z 2m(3) 


D 2m( 2 ) 

D ?m (3) 


A W -2 %-2 C N -2 Z 2 m ^ N_2 ) I ) 2 m ^" 2 ) 

%-l %-l %-l D 2m( W -D 

L A U ' % J W W ) ^(*0 

with m = 0,2. 

To solve these equations, Gaussian elimination (sometimes referred to as 
Potters' method, ref. 42) was employed on a digital computer - . 


Evaluation of b 

As mentioned in the derivation of the governing equations, the 
determination of b requires identification of the 'secular terms occurring 

i 

in the general solution of the third perturbation. To accomplish, this, 

the general solutions of equations (39) are substituted into the equilib- 

ny 

rium equations, and terms whose coefficients' are either cos — or 
sin “ are collected. The resulting equilibrium equations associated 
only with secular terms can be written as 


A 11 Z 31 “ bA cr c l I) 31 + g 1 E 31 


in which 
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with' 


and. 
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Following the procedure for development of boundary conditions 

d 31 (i) = S 51 (1) - a 1 a£-4> 31 (i) 
s 31 (it) - d 3 i(n) - (N) 



(72a) 
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and 


E 5 l(l) = - AiA^iCl) 


E 51 (N) - E 31 (H) - CiCi 1 E 51 (N) 


(72b) 


where A^, C-^ and C-^ are evaluated for m = 1. 

For simple support boundary conditions 
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(73b) 


For clamped boundary conditions, D^i and Ey_ become null vectors. 


Thus the vectors and Ey can be written as 
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Finally, using the equation for b (eg.. ( 50 )) previously developed 


b = 


-2t 2 AZ^E 51 




(U^ + vf + wf )dx 


( 75 ) 


Program SICK 

The numerical computations to determine the sensitivity index b 
were performed on a CDC 6600 series computer using a Fortran program 
entitled SICK (Sensitivity Index for Cylinders ex Koiter). A listing 
of this program is contained in appendix A. In the program, four key 
subroutines were employed: MATRIX, P0TTER, DIF, and SIMP. MATRIX is a 

CDC general purpose program which permits addition, subtraction, multipli- 
cation and inversion of matrices. P0TTER is a Langley Research Center 
program which solves banded matrix equations with a three-vector bandwidth 
by Gaussian elimination (Potters' method, ref. k-2 ). DIF is a Langley 
program which numerically differentiates a function using a 10 -point 
Gauss quadrature method. SIMP is an elementary program which integrates 
a function by Simpson's rule. 

The numerical results were ' found to converge with about 40 finite- 

t 

difference elements. The results presented herein were obtained with 
60 elements. In the numerical approach employed, a major consideration 
was computer storage required to solve,, the .equations resulting from the 

- * f 

second perturbation (eqs. ( 70 )') using the subroutine PATTER. The choice 
of 60 elements was largely dictated by the fact that it led to a require- 
ment of 70,000 octal storage units by program SICK. This storage was the 



maximum obtainable on the Langley remote terminal system. Running times 
for a typical problem were quite excellent with the solution requiring 
about 25 seconds’ computing time. 



IX. INVESTIGATION OF RELATIVE IMPERFECTION SENSITIVITY 


OF FIBER-REINFORCED CYLINDERS 


The imperfection sensitivity of certain fiber-reinforced cylinders 
was investigated with the Koiter index b determined from both the 
present theory (eq. (75)) as well as the theory of references 23 and 24. 

To employ the latter theory, equation (44) was employed together with 
the normalization described in equation (65). Three types of cylinders 
were considered: glass-epoxy, boron-epoxy, and aluminum reinforced with 

an overlay of boron-epoxy. The cylinders considered had a fixed radius- 
thickness ratio R/t of 100 and a length- radius ratio a/r of 0.7- 
The helical wrap angle of the cylinders was varied in an attempt to 
identify configurations of minimum imperfection sensitivity. 

The investigation of the sensitivity index was limited to a few 
specific configurations because of the tedious and expen sive nature of 
the calculations. It was found that buckling loads using the theory of 
reference 4l could be extracted by computer with only moderate difficulty. 
However, many additional computations were required in certain instances 
to obtain a buckling mode shape sufficiently accurate to pursue post- 
buckling calculations. 

The elastic constants for the single fiber-reinforced layer in each 
cylinder were estimated by the procedure summarized in an appendix of 
reference 1. In this 1 procedure elastic constants along principal axes 
of unidirectionally reinforced material' were obtained from the theoretical 

* ’t 

study by Hashin and Rosen (ref. 33). The upper bound of the random 
fiber array model described therein was employed in the present study. 
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Elastic constants for the helically wrapped layers with filaments 
oriented at ±a from the cylinder axis were estimated using transforma- 
tion equations of orthotropic elasticity together with the concept of 
effective stiffness of cross -plied laminates (refs. 35 .> 36 ). This 
technique has been shown to give predictions in reasonable agreement 
with experiments to determine elastic constants of glass-epoxy (ref. l) 
or boron-epoxy (ref. 4) cylinders. Constituent material properties and 
layer properties are summarized in table I. 

G lass-Epoxy Cylinders 

The cylinders investigated were composed of a single helically 
wrapped layer (lot,). The constituent properties and fiber volume fraction 
vp selected for the investigation are typical of well-made filament- 
wound cylinders. Buckling predictions for cylinders whose wrap angle 
varies from 0 to 90° are shown in figure 3. The predictions are based 
on boundary conditions of simple support in which 


M X1 = 0 • 



The solid curves of figure 3 labeled "consistent theory" -were obtained 
from reference 4l from a theory which considers nonlinear effects of 
change in shape of the cylinder due to load prior to buckling. The dash- 
dot curves of figure 3 labeled "classical theory" were obtained from the 
theory of reference 2 which assumes a membrane state of stress in the 
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cylinder prior to "buckling and ignores changes in shape due to loading. 

The dashed, curve labeled "axisymmetric collapse" can he obtained from 
either classical theory with n equal to zero (ref. 2) or by a limiting 
solution of equations governing axi .symmetric deformations prior to 
buckling (e.g. eq. ( 17 ))-. 

From figure 3 it 'can be seen that on the basis of classical theory, 
cylinders with helical wrap angles from 30 to 60 ° are expected to collapse 
axi symmetrically. Consistent theory predictions, however, do not suggest 
this mode, but rather gives loads based on asymmetric modes which are 
substantially lower than axisymmetric collapse predictions. The relative 
position of axisymmetric collapse predictions to consistent theory 
buckling predictions is extremely important in the present considerations 
of imperfection sensitivity, because the sensitivity index b is 
undefined when- -axisymmetric collapse is. the predicted mode. of failure.. 
Thus, because consistent theory predictions are less than axisymmetric 
collapse predictions for all values of the helical wrap ’angle., the 

. <f , 

r \ 

sensitivity index can be defined everywhere. - 

Substantial differences between consistent theory and- -classical . 

t t ' 

theory predictions are indicated in figure 3* The differences are 

1 i +1 

functions of the helical wrap angle and are greatest, for a equal to 
4-5°. At this angle, consistent theory predictions are about 80 percent’ 
of classical theory. For isotropic cylinders with the same radius- . 
thickness ratio, length-radius ratio and boundary conditions, consistent 
theory predictions are also about 80 percent of classical theory 
(ref. by ) . The results suggest that prebuckling deformations in 
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filamentary shells possessing high shear stiffness are an important 
consideration for accurate buckling predictions. 

Imperfection .sensitivity indices for the simply supported glass- 
epoxy cylinders are shown in figure 4. Calculations -of the index b 
were made using the present theory (solid curve }, as well as the theory 
of references 23 and 24 (dashed -curve) . The results for the sensitivity 
index b 'have been normalized, by the index bpso for an isotropic 
cylinder having the same R/t and a/r -as the present cylinders. The 
values of b^ so obtained from the theory of references §3 and, 2b agreed 

* s 

well with unpublished data (ref. 44) by J. W. Hutchinson. In plotting 

results on figure 4, computed data have been linked by • straight -line 

\ 

segments to indicate the points at which discrete calculations’ were made. 

T * / 

The results shown in figure 4 indicate similar • trends' in* relative * 

* ■*. i ' 

imperfection sensitivity predictions based on the two theories, considered. 

' * T 
The differences in magnitude of the imperfection sensitivity indices 

* f ’ ~ ' * 

are attributed to differences in definition of the perturbation parameter 
g. In figure 4, the normalized sensitivity- index .h/bpso takes on 'its , 
smallest values for helical wrap angles frcm 30 to 60°. Thus, both 
theories suggest that glass-epoxy cylinders with these wrap angles are 
desirable configurations since they have a* relatively low index. Because 
of the high index values for a equal to 75°, this configuration is 
presumed to he undesirable. 

The results of a recent imperfection- study by Khot (ref. l4) using 
a classical approach which ignores load-induced deformations prior to 
buckling hut considers an initial imperfection of specified shape , should 
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also be noted. For three-layer glass-epoxy cylinders possessing 0, 

-a, +a, or -a, +a, degree configurations , Khot found that maximum 

compressive buckling loads for small amplitude imperfections occur for 
shells where a is 45°. This result coupled with the present results 
suggests that the use of wrap angles near 45° in glass filament-wound 
cylinders may improve their buckling strength due to reduced sensitivity 
to initial imperfections. 


Boron-Epoxy Cylinders 

Buckling predictions for single layer, simply supported, boron- 
epoxy cylinders are presented in figure 5- In contrast to the glass- 
epoxy cylinders, consistent theory predictions for boron-epoxy cylinders 
indicate that axisymmetric collapse will occur for some wrap angles. As 
mentioned previously, the Koiter sensitivity index is undefined for 
these wrap angles. Fote from- figure 5 that consistent theory 'and 
classical theory predictions do not yield identical wrap angles for 
axisymmetric collapse. For consistent theory, collapse is predicted 
for wrap angles from 15 to 55° whereas for classical theory, collapse 
is predicted from 30 to 60 °. 

As was the case in glass-epoxy cylinders, agreement between , 
consistent and classical theory for boron-epoxy cylinders is a function 
of the helical wrap angle. Maximum differences between the 'two theories- 
occur for wrap angles of 45 and 75°. For these wrap. angles, consistent 
theory yields predictions which are roughly 85 percent of classical • 
theory. Both classical and consistent theory buckling loads are 






70 


relatively high for cylinders with a wrap angle of 15°. The magnitude 
of load is believed to be a consequence of the high axial extensions! 
and axial bending stiffness of boron filaments in this configuration, 
coupled with unusually large values of the product of Poisson's ratios 
(see table I). Since most of the structural stiffnesses of the 
shell are proportional to 1/ji - p x My), large values of |i x Py result in 
increased stiffness. 

Imperfection sensitivity indices for boron-epoxy cylinders are 
presented in figure 6. Again, agreement in trend between the present 
theory and the theory in references 23 and 24 is good. Because of 
axi symmetric collapse, the sensitivity index is undefined for helical 
wrap angles where 15 ~ a "h 35* Note from the figure that values of 
sensitivity indices adjacent to the zone of axisymmetric collapse may 
be quite different. For a equal to 10°, the sensitivity index (as well 
as the buckling load) is relatively high. On the other hand, for a 
equal to 4o°, both the sensitivity index and the buckling load are low. 
These results demonstrate that cylinders which buckle near axisymmetric 
collapse loads are not necessarily desirable from the standpoint of 
reduced imperfection sensitivity. 

A comparison of the magnitude of the sensitivity index for glass- 
epoxy cylinders (fig. k) with boron-epoxy cylinders (fig. 6) in regions 
of minimum imperfection sensitivity reveals that boron-epoxy shells -are 
less sensitive. Thus, shells with wrap angles around 45° would appear 

r \ 

to he an attractive configuration in boron-epoxy cylinders. Ehot’s 
imperfection study (ref. 14) indicated that the 45° wraps in 0, -a., +a. 
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and 90 , -a,, +a, configurations gave higher buckling loads for small 

amplitude imperfections. He also concluded that boron-epoxy cylinders 
were less sensitive to initial imperfections than glass-epoxy cylinders. 

Boron-Aluminum Cylinders 

In the studies of glass-epoxy and boron-epoxy cylinders, it was 
concluded that cylinders with wrap angles of V?° had reduced imperfection 
sensitivity. A shortcoming of this wrap angle in the simple configuration 
investigated {t a) is its relatively low buckling load compared to other 
angles, especially in the case of boron-epoxy cylinders (see fig. 5) • It 
seemed desirable to investigate more complex configurations, perhaps 
with other materials, which might possess both reduced imperfection 
sensitivity while still retaining relatively large buckling strengths. 

An interesting concept in fiber-reinforced structures now being 
pursued by the MSA is that of employing a metal '’substrate and overlaying 
on it boron-epoxy tapes. , This concept of selective reinforcement is 

i - , * t ‘ * 

attractive in that it reduces practical attachment problems and .costs- in 
using expensive filamentary material, yet at the same time provides good 

t 

possibilities for structural weight savings. A recent study of tubular 
columns (ref. ^5) in which metal tubes were overlayed with axially f 

oriented boron filaments has produced some spectacular results for . .* 

/ ’ * *' 

compressive buckling. 

Because of current interest in the selective reinforcement concept, 
boron-aluminum cylinders were investigated. The cylinders were composed 
of two layers of equal thickness, one of aluminum, the other of boron- 
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epoxy with filaments at Based on total thickness, the cylinders 

had the same radius-thickness ratios as the glass and boron cylinders 

already mentioned. The boron-epoxy layer was located on the outer 

surface of the aluminum. Preliminary calculations using classical 

buckling theory indicated that substantial differences in buckling load 

occurred depending on whether the boron layer was located on the inner 

or outer surface of the aluminum. The highest buckling loads were 

developed for the -configuration chosen. This effect of layer sequence 

is the counterpart of the stiffener eccentricity effects discussed in 

• * 

references 9 and 38 . * 

Buckling predictions for the boron-aluminum cylinders are presented 

“t 

in figure 7 . Although classical theory predictions suggest axisymmetric • 
collapse for 30 5s a S consistent theory predictions yield only 
asymmetric modes so that the imperfection sensitivity index 'exists , 
everywhere. Note from- figure 7 that substantial differences exist between 
consistent and classical theory predictions.; /The largest differences 
occur for -a equal to zero where the consistent theory prediction is 
about 73 percent of the classical prediction. The maximum load carried 

-» r 

by the boron-aluminum cylinders occurs for a wrap angle of 6o°. The 
reason why cylinders with this wrap angle perform so- well is not obvious, 
since the effects of prebuckling deformations and eccentricity in 
multilayered -shells have not been thoroughly investigated. The only 
unique feature of the 6o° wrap is its large value of |iy (see table i). 

Imperfection sensitivity indices for boron-aluminum cylinders are 
presented in figure 8. Again, reasonable agreement in trends occur for 
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Figure 8.- Imperfection sensitivity indices for boron -aluminum cylinders. 
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the two versions of the Koiter theory. An interesting and encouraging 

*■ $ | 

result is the negative values of the normalized imperfection sensitivity 

o ^ 

index which occur for a equal to 90 . Negative values indicate a' 

* i ’ 4 

positive slope of the load- shortening curve for the cylinder in the post- 
buckling region. While this result is somewhat unexpected, it maybe, 
explained by observing that for a equal to 90 °, the cylinder behaves 
as a ring -stiffened shell with very closely spaced, stiff (boron-epoxy)' 
rings located on the outer aluminum surface. There is- some evidence 
in the literature that ring- stiffened shells have unusual postbuckling 
behavior. In postbuckling studies, both Thielemann (ref. 10 ) and 
Shang (ref. 12) have found relatively large minimum postbuckling loads 
for ring- stiffened cylinders, thereby suggesting that they possess little 
imperfection sensitivity. 

Relatively small imperfection sensitivity is also indicated for the 
^5° configuration. The magnitude of the normalized index, however, is 
much closer to unity than that for boron and glass cylinders with the 
same wrap angle. Thus imperfect boron-aluminum shells in this config- 
uration would he expected to perform about the same as' imperfect 
isotropic shells. 

A somewhat frustrating law appears to evolve from the present 
investigation. As the compressive-load carrying capability of fiber- 
reinforced cylinders is improved by varying the filament orientation,, 
the structure's sensitivity to imperfections becomes larger. An 
attractive fiber-reinforced cylinder configuration has been identified, 
however. Because of a positive slope in its postbuckling load- shortening 
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curve and its reasonably good load- carrying ability, the 90 ° or circum- 
ferentially wrapped boron-aluminum cylinder appears to be the most 
desirable of the configurations investigated herein. 'Whether its 
performance will be superior to configurations with greater potential 
load-carrying capability but greater imperfection sensitivity can only 
be determined by experiment. 



X. EXPERIMENTAL INVESTIGATION - 


Test Specimens 

To investigate the buckling behavior of fiber-reinforced cylindrical 
shells, compression tests were conducted on 12 filament-wound, glass- 
epoxy cylinders. The test specimens had an inside diameter of J>0 inches, 
were 30 inches long, and had a nominal thickness of 0.072 inch. The 
cylinders were wrapped in helical layers composed of two half layers, 
with filaments oriented at angles +cx or -a measured from the cylinder 
axis. The variables in the test program were the helical wrap angle a 
and- the cylinder wall configuration. 

Six of the cylinders had walls composed entirely of five helically 
wrapped layers in the sequence lx, ix, lx, lx, lx. The walls of the other 
six cylinders were composed of four helical layers and four circumferen- 
tially wrapped layers in the sequence ia, 90, lx,' 90, lx, 90, lx, 90, where 
the last circumferential wrap (+90) formed the external surface of .the 
cylinder. Photomicrographs of the two cylinder wall configurations are 
shorn in figure 9- For cylinders containing circumferential wraps, the 
helical layers at lx can be seen to be roughly twice, the thickness of 
the circumferentially wrapped layers. 

i *, , 

Table II contains material data for the ' cylinder, constituents. The 

, ■- * 

properties shown for glass fibers are typical of those usually employed 
for type E glass filament. Values of modulus and density shovpi'for the 

V * ' ■ 

ERL-2256 epoxy resin were obtained from compression tests of .four v 3- inch- 

« t 

diameter, solid resin cylinders about 6 inches long. The cylinders ' were 
cured with the same curing cycle as the test specimens, that is, 250- F 

1 1 < 

for 2 hours. The solid cylinders were instrumented with’ both Tnckerman 
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optical strain gages and differential transformer gages. i The four cylin- 
ders were tested in compression at a loading rate of 20,000 pounds per 
minute. A typical stress-strain curve generated on an X-Y plotter by the 
transformer gages is reproduced in figure 10. All cylinders appeared to 
be capable of straining indefinitely .at maximum load so that the tests 
were terminated when about 6 percent compressive strain had been experi- 
enced by the cylinder. Wo apparent damage could be observed at the 
conclusion of the test. Based on the test results, EEL- 2256 appears to 
have a proportional limit of about 10 ksi in compression. 

The dimensions and helical wrap angles of the glass-epoxy test 
specimens are presented in table III. The values of total wall thickness 
t shown are the average of several measurements taken at random loca- 
tions. The scatter in individual measurements was ±4 percent of the 
value tabulated and is attributed to the irregular outer surface of the 
specimens. It should also be noted that the outer surface of each 
specimen was reinforced at each end with circumferentially wrapped 
layers about 1-1. 5 inches wide with thickness equal to that of the shell. 
The purpose of this reinforcement was to prevent local crippling due to 
the end loads. 

The fiber volume fraction vp, the resin, volume fraction v r , and 
the void volume fraction v v , expressed as percentages, are also 
presented in table III. The fiber volume fractions were determined by 
elevated temperature exposures of two or three coupons cut from* the wall 
of each specimen. The coupons were subjected to a temperature of 1100° F 
for a period of 3 hours to permit virtually , all- of the resin in the- 
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coupon to decompose. The remaining glass was weighed and the equivalent 
fiber volume fraction was calculated using the values of density ■ 
presented in table II. Prior to heating, the. volume of each coupon was 
determined using Archimedes’ principle. The resin and void ‘contents 

i " ' t 

were then calculated using the glass volume fraction together with the ■ 

total volume and the density of resin quoted in table II. Although- * - 

' \ * 

widely used, this method of determining void^ content is not extremely 

i j . ’ 1 t 

■ ' - . i 

accurate, so that the true void content probably can be determined only 
to the nearest half a percent. It should be noted that the volume 
fractions of table III together with studies of photomicrographs of < * 
specimen walls suggest that, the quality of the test specimens was only 
fair. Well-made glass-epoxy, filament-wound cylinders usually possess 
fiber volume fractions of about 65 percent while those of the present 
specimens are about 51 percent. 

Test Procedure 

Each of the 12 cylinders was instrumented with l6 wire-resistance 
strain gages with a gage length of 6 inches. The gages were bonded in 
back-to-back pairs on the surfaces of the cylinder and oriented in the 
axial and circumferential directions at the middle of the cylinder about 
90° apart. The purpose of the axial gages was to measure initial strains 
from which the axial extensional modulus could be calculated, and to 
detect any bending deformations which might suggest buckling of the' 
cylinder wall. The circumferential gages were present to measure -strains 
from which Poisson's ratio il x , associated with loading in the axial 
direction, could be determined. Additional Instrumentation consisted of 
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two linear differential transformers which were used to track the 
relative motion of the testing machine platens. 

For testing purposes, a 3 / 8 - inch plywood bulkhead was inserted at 
each end of the cylinder to ensure that the ends of the specimen would 
remain circular. . The cylinders were tested flat-ended in compression in 
the Langley Research Center's 1,200 kip-capacity testing machine. The 
loading head of the testing machine was carefully alined to promote 
uniform loading of the cylinders by checking strain distributions in the 
axial strain gages under small loads. The buckling test consisted of a 
■single load cycle to failure at a loading rate of 5,000 pounds per 
minute. During' the test, data were recorded at a virtually continuous 
rate on the Langley Central Digital Data Recording System. 

Test Results 

Elastic constants .- The data obtained during the tests were reduced 
by computer and presented in the form of load- strain or load-displacement 
curves for each specimen. A sample plot of data for cylinder 10 is shown 
in figure 11. Data from the axial and circumferential gages have been 
plotted so that back-to-back pairs of gages have the same origin. For 
convenience, only the magnitudes of strain have been plotted. The average 
slope .of the axial strain data was linear- over a wide range of the loading 
history, and values of slopes determined from individual pairs of gages 
were in good agreement. The average of the linear slopes of the axial 
gages expressed as a Young's Modulus is presented in table III as E x . 1 
Values of E x deduced from ‘the transformer measurements of overall 
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shortening after corrections for the cylinder end reinforcement were in 
reasonable agreement with the axial strain-gage data. 

Data from the circumferential gages were somewhat more nonlinear, 
and in .some cases, substantial differences in individual gage slopes 
existed. The initial slopes of the circumferential gages were used to 
determine Poisson* s ratio of each cylinder by dividing E x by the 
average of the circumferential slopes. The Poisson's ratio associated 
with circumferential strain induced by axial load (T x is also, shown in 
table III. The .underlined values in the table are values in which -large 
variations in individual circumferential gage slopes were present. 

The test data for elastic constants are' compared in figures' 12 to 15 

. i 

with theoretical predictions for total extensional stiffness of the 

. * 

cylinder wall based on references 55 and 55- It can be seen that theory 
and experiment are in good agreement for the extensional modulus of the 
cylinder wall and in fair agreement for Poisson* s ratios- in spite of the 
scatter in data. The results suggest that the techniques employed « 

herein for predicting elastic constants yield good .engineering estimates 
of structural stiffnesses. - 5 

Buckling .— Failure .of most of the test specimens was accompanied by 
a loud report and the appearance of two tiers of large, diamond- shaped 
buckles uniformly distributed around the cylinder. As suggested in the 
load- strain curves of figure 11 which are plotted up to maximum load, no 
strain reversal was apparent prior to failure. The maximum load carried 
by the cylinder A ex p is presented in table III. A photograph of a 
buckled cylinder is presented as figure 16. An- observation' associated 
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Figure Pj.- Poisson’s ratio associated with radial expansion induced by axial 
loading of test cylinders with helical wraps (cylinders 7 to 12). 
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with the two specimens having 30° helical and circumferential wraps 

(cylinders 3 and. if) should be noted. On both these specimens, buckles 

appeared only on roughly one-half' of the, shell surface. As a matter of 

routine, before unloading, the specimens were prepared for photography 

by removing instrumentation wiring. In both cylinders, during this time 

period, a uniform buckle pattern suddenly developed, replacing the 

original, nonuniform pattern with no large changes in load or platen 

shortening. Since all specimens were tested with the same procedure, it 

> 

appears that nonunifoim buckle patterns might be' an - intrinsic charac- 

- f . 

teristic of that wrap configuration. ’ The particular configuration 

' * % 

(±30, 90) happens to possess two.- dimensional isotropy with regard to 
bending and extensional stiffness (but not .complete isotropy due to lack 
of symmetry about the shell midsurface). It seems, experimentally, that 
less isotropic wrap configurations possess- more uniform postbuckling 
patterns. 


Comparison of Theory and Experiment 
The experimentally obtained buckling loads for the test cylinders 
with alternating helical and circumferential wraps are compared with 
theoretical predictions in figure 17 and table III. The theoretical 
loads A cr were based on the simple support boundary conditions 
previously mentioned. The computations were performed using the con- 
sistent theory taking into account nonlinear prebuckling deformations 
(ref. 4l). A few calculations using the clamped boundary conditions 
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Figure 17.- Comparison of theoretical 'and experimental buckling 
loads for test cylinders with alternating helical and 
circumferential wraps. 
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Figure 18. - Imperfection sensitivity indices for test cylinders 
with helical and circumferential wraps. 




suggested that differences "between clamped and simple support calcula- 
tions' were small because of the length of the test specimens. 

The trends shown - in figure 17 indicate that the greatest difference 
between theory and experiment occurs for a wap angle of 30° in 
cylinders 3 and ij-. As mentioned previously, these cylinders possess 
two-dimensional isotropy with respect to bending and "extensional stiff- 
ness, and behaved somewhat differently at buckling than the other, less 
isotropic cylinders in the test program. ’ If the data in figure 17 for 
wrap angles of 30° and ^5° are compared, it is evident that, for 
cylinders with nearly identical radius -thickness ratios, agreement between 
theory and experiment is much better for the less isotropic configura- 
tion (that is, 45°). 'Moreover, from table III it can be seen that 
analytically, the 30°' configurations were expected to. carry more load 
than the k^° configurations, whereas experiment ally the reverse was true. 
Presumably, the reason for this, .unusual phenomenon is associated with the 

4 1 * , 

1 t-5° -configuration being less imperfection- sensitive than the nearly 

, t % t 

isotropic configuration. Thus, the present tests .offer 'experimental 
evidence that certain fil ame nt orientation's can enhance the buckling 

s " 1 „ r 

strength of a cylinder due to- their decreased imperfection sensitivity. 

* * 

Direct comparisons between the 15° and 30° -configurations are complicated 


due to the differences in radius -thickness ratios. 



Imperfection sensitivity indices for cylinders 1 to 6 are presented 
in figure 18. The indices have been normalized by the index bj_ g0 for 
isotropic cylinders of the same geometry. Results are shown for both the 
present theory and the theory of references 23 and 2k. Agreement between 
the two theories is fair except for a wap angle of 15°. The large value 
obtained for the sensitivity index at this angle from the referenced 
theory is not presently understood. 

The trend of data from the present theory, however, is quite 
encouraging in that maximum values of the sensitivity index occur at 30°. 
Thus the Koiter theory predicts that cylinders 3 and 4 should be most 
imperfection sensitive. This result is in agreement with the experimental 
data, if it is hypothesized that the more sensitive buckling of a struc- 
ture is to geometric imperfections, the lower is its value of A eX p/A cr . 

It was expected that the sensitivity index for wap angles of 30° would 
be about the same as an isotropic cylinder so that b/b igo 1. The 
lower value of the normalized index shown in figure l8 is believed to be 
a consequence of the antisymmetry of the wall of the cylinder. 

The experimentally obtained buckling loads for the test cylinders 
with helical waps are compared with theoretical predictions for simply 
supported cylinders in figure 19 and table III. From figure 19, the only 
obvious trend in the data is the lack of reduction in values of l e xpAVr 
as R/t is increased from 1 60 to 223- In isotropic cylinders, of course, 
some reduction would have been expected even in this relatively small 
change of R/t. 
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For the 15° configuration, the data for cylinders 7 and 8 are 
inconsistent. From table III it can be seen that even though 
cylinder 7 was thinner than cylinder 8 , it carried more load. Volume 
fractions and thicknesses were carefully reinvestigated, but no reasons 
were apparent to explain the discrepancies in the data. 

Imperfection sensitivity indices for cylinders 7 to 12 are presented 
in figure 20. Agreement between the two theories presented is fair 
except for large disagreement for a helical wrap angle of h^°. This 
anamoly was unexpected, since the agreement for this wrap angle for 
glass-epoxy cylinders with' R/t = 100 was much better (see fig. k). In 
figure 20, variations in R/t make comparisons between wrap angles 
difficult. However, the trend of the indices for the present theory 
suggests that cylinders with wrap angles of 30 ° and b^>° should perform 
better than isotropic cylinders of the same radius -thickness ratio. The 
high values of A eX p/ 7 \ cr shown in figure 1% especially for the largest 
R/t (cylinders 11 and 12), tend to support this prediction. The 
sensitivity indices for cylinders with the lowest radius -thickness ratio 
suggest that their buckling behavior should be about the same as isotropic 
cylinders. The values of Aexp/Tkr for cylinders 7 and 8 , and the scatter 
in data seem to support this prediction also. Thus the trends in buckling 
data suggested by the sensivitity indices of the present theory are in 
fair agreement with experimental buckling data. 



XI. CONCLUDING REMARKS 


A theory has been developed to assess the, sensitivity of buckling 
of axially compressed, fiber-reinforced cylindrical shells to small 
geometric imperfections in the shell wall. Following the work of Koiter 
(ref. 19 ), an index of imperfection sensitivity is obtained by investi- 
gating the character of the initial postbuckling region for a geometric- 
ally perfect cylindrical shell. The theory presented herein differs from 
existing applications of Koiter theory in that it is developed by using 
techniques of perturbation theory. The theory represents an extension to 
cylindrical shells of the plate postbuckling theory derived by Stein 
(ref. 25 )• Using Donnell-von Karman strain -displacement relations, the 
postbuckling equations are obtained for a multilayered, orthotropic 
cylinder in which the principal axes of orthotropy in each layer are 
alined with the longitudinal and circumferential coordinates of the shell. 
The imperfection sensitivity index is identified by suppressing secular 
terms in the perturbation and employing a suitable definition of the 
perturbation parameter. The governing equations were solved by finite 
difference techniques and by exploiting Gaussian elimination. A computer 
program to solve the large systems of simultaneous equations resulting 
was developed for the CDC 6600 series computer. 

The results of an investigation of the imperfection sensitivity 
index for fiber-reinforced shells have been presented. Glass-epoxy, 
boron-epoxy, and boron-aluminum shells having various helical wrap angles 
were studied in order to identify wrapping configurations of minimum 
imperfection sensitivity. In conducting the investigation, both the 
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present theory and an extension of existing theory to multilayered, 
orthotropic shells were employed. Agreement between predicted sensitivity 
indices based on the two theories was reasonably good. Computed results 
for buckling in this investigation suggest that substantial differences 
(up to 27 percent) can exist between classical and consistent theory 
predictions for buckling, depending on the wrapping configuration. The 
differences are attributed to the neglect of load-induced prebuckling 
deformations in classical theory. 

The studies of helically wrapped cylinders indicated that k^>° was a 
desirable wrap angle in both glass-epoxy and boron-epoxy cylinders. For 
wraps with small imperfection sensitivity, boron-epoxy shells appeared to 

f 

be more desirable than glass-epoxy shells. The most desirable shell con- 
figuration uncovered in the investigation was that of an aluminum shell 
overwrapped with a circumferentially wound boron-epoxy layer. This con- 
figuration was found to possess a positive postbuckling slope and appears 
to behave much like a ring -stiffened cylinder. 

The results of an experimental program to investigate the buckling 
behavior of fiber-reinforced, glass-epoxy cylinders have been presented. 
Compressive buckling tests were conducted on twelve 30-inch-diameter 
filament -wound cylinders with various wrap angles. Elastic constants 
derived from data prior to buckling were compared with analytical pre- 
dictions. The agreement between theory and experiment suggests that the 
approximations employed to estimate structural stiffnesses of the cylin- 
der wall are valid. Cylinder buckling loads varied from 50 to 90 percent 
of consistent theory predictions, depending on wrap angle. Comparisons 
of calculated sensitivity indices with the cylinder buckling data show 
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that the indices are in fair agreement with the experimental data trends. 
Both the theoretical and experimental results reported herein demonstrate 
that, certain fiber orientations can enhance the buckling strength of 
geometrically similar cylinders. These configurations appear to possess 
reduced sensitivity to geometric imperfections present in the cylinder 
wall prior to loading. The present results suggest the desirability of 
further experimental investigations into the buckling behavior of config- 
urations such as the circumferentially wrapped boron-aluminum cylinder. 
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XV. APPENDIX A 


In this section program SICK is described in detail. A flow diagram 
for the program is presented as figure 21. The input data required are as 
follows : 


R 

A 

DELBAR 


BC 


NXB 

•N 

EX(l) 

ey(i) 

NUX(l) 

NUY(l) 

GXY(I) 

H(I) 

LAYER 

ST 

NSTAT 


radius of cylindrical shell to> reference surface 
length of cylindrical shell 

distance from inner surface of shell to reference surface 
(for the present calculations t/2) 

boundary condition code BC = 1 simple support 

= v i = w i = M x x = °) 

BC = 3’ clamped 

= vq = wq = W 1>x = 0) 

applied axial compression load at buckling, A cr 

number of full waves in circumferential buckle pattern, n 

elastic modulus for shell layer in axial direction, Ej 

elastic modulus for shell layer in circumferential direction, 
T?i 

Poisson's ratio of shell layer associated with load applied in 
axial direction, p^ 

Poisson's ratio of shell layer associated with load applied in 
circumferential direction, pi 

shear modulus of shell layer, G^y 

thickness of shell layer 

total number of layers in shell wall 

code to detect stiffening ST = 0 cylinder is unstiffened 

ST = 1 cylinder is stiffened 

number of finite difference stations, N 
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M0RM code to specify normalization H0RM = 0 |zfj is normalizing 

factor 

R0RM = 1 W max /t is 

normalizing factor 

IPRINT code to specify type of output printed 
IPRINT = 0 gives minimum output 
IPRINT = 1 gives maximum output 

If ST = 1 the following additional input is required: 

ER modulus of ring 

AR area of ring 

IR moment of inertia of ring about its centroid 

ZR distance from reference surface to ring centroid (+ only if 

ring is external) 

GRJR torsional stiffness of ring 

RS ring spacing 

ES modulus of stringer 

AS area of stringer 

IS moment of inertia of stringer about its centroid 

ZS distance from reference surface to stiffener centroid (+ only 

if stringer is external) 

GSJS torsional stiffness of stringer 

BS stringer spacing 

Some conversions between the program notation and the notation of 
the paper which may not be self-evident are: 


WOX 

WOXX 

WOXL 



w 


* 

0 


WOL 
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M0DRES modified residual resulting from taking determinant of left- 
hand side of buckling equations (see left-hand side of 
equation (62 ) ) 

A listing of the program together with a sample case for the circum- 
ferentially overwrapped boron- aluminum cylinder follows. 




SET m = 0 


INTEGRATE TO FIND b BY 
PRESENT THEORY AND BY EXTENSION 
OF REFERENCES 23 AND 2k 


Figure 21. - Flow diagram for program SICK, 
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C PROGRAM SICK (SENSITIVITY INDEX FOR CYLINDERS EX KOITER ) 

PROGRAM SICK ( INPUT , OUTPUT ) 

10 FORMAT ( ///# IMPERFECTION SENSITIVITY OF AXIALLY COMPRESSED CYLINDE 

1R*/ * 

1 CARD-SYKES A2338 RDF367 OCTOBER, 1968) 

11 FORMAT!//# LAYER*, 6X*EX*,9X*EY*,9X*NUX*,8X*NUY*, 8X*GXY#, 9X#H# ) 

12 FORMAT ( 15,-3 X, 6 El 1 .3 ) 

13 F0RMAT(//13X#11#,17X#12*,17X#22*,17X#66#/) 

14 F0RMAT(1X*B*,I3,E16.8,3(3XE16.8) ) 

15 FORM AT ( 1X*C#» 4 ( 3XE1 6.8) ) 

16 FORMAT ( 1X*D*,4 (3XE16.8) ) 

17 FORMAT ( 1 X#K#,4 (3XE16. 8 ) > 

18 FORMAT (#1#,4X#STATI ON#, 9X*W0* , 17X#W0X# , 17X#W0XX#, 16 X#WO XL#, 
116X#W0L# ) 

19 FORMAT! I 8,4 E20. 8) 

20 FORMAT!////* M0DRE$=*E16. 8) 

21 FORMAT !V/5X*W1 VALUES#/) 

22 FORMAT! 15, E20. 8) 

DIMENSION WO! 61),VfOX( 61),WCXX( 61),W0XL( 61 ) , EX! 20) , EY! 20') , NU X! 2 
10 ) ,NUY ( 20) ,GXY(20),H(20),B11(20),B12(20),B22(20),S!20),SS(20),SSS( 
220) ,866(20) , NYOL ( 61), WOL! 61) 

DIMENSION AA( 61,4,4),B! 61,4,4),C( 61 ,4, 4) , AB ( 2 ,4 , 4) , BB( 2 , 4, 4) , 
1CB! 2, 4, 4), DUMA (4, 4) , DUMB (4, 4) ,DUMC(4,4) ,DUMAB(4,4) , DUMB 3(4, 4) , 
20UMCB (4,4),B1(4,4),C1(4,4),AN!4,4) ,BN(4,4» ,'P'( 61,4,4),R1( 61,4,4), 
3PIM1(-4,4),RI(4,4),PI(4,4),U1( 61), VI! 61),W1( 61), Ml! 61),Z1( 61,4 

4) » Z ( 4) ,'PDC 4,4.).’,- Z1SQ(1), ZD(4) 

DIMENSION X! 61) ,U1X( 61),V1X( 61),W1X! 61),W1XX( 61),W02L( 61), 
•1W13X ( 61-),-W03L( 61) ,NX1 ( 61), NYU 61),NXY1! 61), D( 61,4), 

2SMAT ( 6T,4,4),DMAT( 61,4),PP! 61,4,4),QQ( 61,4),Z20< 61,4), 

3Z'22( 61,4) , U20 ( 61), V20 ( 61),W20! 61),M20( 61),U22( 61),V22( 61), 
4W22 ( 61) ,M22( 61),W20X( 61),U20X! 61),U22X( 61),V22X( 61),W22X( 61 

5) ,W20XXl: 61 ) , W22XX! 61),NX20( 61),NY2C( 61),NX22( 61),NY22( 61), 
6NXY22 ! 61 ) r DB (4 > , DUMD(4), W20XXX! 61), W22XXX ( 61) 

DIMENSION DP3! 61,4),EP3( 61 ,4) , DUMDE ( 244) , Z1T! 244 ) ,.FOFX( 611, 
1XTE(1) ,XTD! 1) -,Z1V(244), DPB3(2,4), EPB3C2,4) 

REAL Ml, NYOL 

REAL IR, IS, NUX,NUY-, Kll, K12, K22, K66, NXB 

REAL NX1, NYl, NXY1 , M20, M22, NX20, NY20, NX22, NY22, NXY22 

INTEGER BC, ST 

REAL M, lAMll, LAM12, LAM13, LAM17 , LAM21 , LAM22, LAM23, LAM24, 

1 LAM32, LAM41, LAM42 , LAM43, LAM45, LAM14, LAM15 , LAM16, 

2 LAM25, LAM26, L AM31 ,. LAM33 , LAM34, LAM35, LAM36, LAM44, 
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3 MOORES 

INTEGER- SUBTRAC , TRANSPO 

NAME L I ST/I N PUT /R, A, DELBAR, B C, NXB, N, EX, EY , NUX ,NUY , GXY, H , LAYER, ST, NS 
1TAT,N0RM, I PRINT 

NAMELIST/ST IFF/ER,AR, IR,ZR» GRJR,RS, ES,AS, IS, ZS, GSJS,BS 
1000 READ INPUT 
PRINT INPUT 
PRINT 10 
PRINT 11 
00 50 L=l, LAYER 

, 50 PRINT 12, L, 6X1 L ) , EY ( LJ , NUXIL), NUYIL), GXY(L), HIL) 

CALCULATE B VALUES 

TEST3=0. 

DO 60 L=i; LAYER 

Bll(L)=EX(L)/( l.-NUXI L} *NUY ( L) ) 

B12(L>=NUY(L.)*B11IL) 

B22IL > = EY(D/( 1.-NUXIU*NUYIL)) 

60 B66I L)=GXY ( L) 

CALCULATE C, K, AND 0 VALUES 

L=LAYER 

CALL SUMIl, L,B11,H,C11) 

CALL SUMIl, L,B12,H,C12J 
CALL SUMIl, L,B22,H,C22) 

CALL SUMIl ,L,B66,H,C66I 

sm=Hii) 

Sl=Sl I)**2-2.*DELBAR*S(1) 

Kll = . 5*B11 I 1)*S1 
K12=. 5*B12 l 1)*S1 
K22=.5*B22-( 1)*S1 
K66=.5*B66< 1)*S1 

S2=S 1 1 ) X>*3-3.^DELBAR*SX 1 )**2+3 .*DELBAR**2*S 1 1 ) 
bll=Bll<l)*S2/3. 

D12=B12I 1) *S2/3. 

022=B22( 1) *S2/3. 

D66=B66( 1) *S2/3. 

IFILAYER.EQ.il GO TO 80 
DO 70 1=2, LAYER 
SI I ) = S( I— 1 l+H( I ) 

SSI I ) = SII) **2-Sl I— 1 ) **2—2 •♦DELBAR*' I S l I )— S 1 1 — 1 ) > 
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70 SSS<I )=S(.I >**3-St I-1)**3-3.*DELBAR*'( S(I )**2-=S( i-1 ) **2 ) + 3.*D,ELBAR** 
12*(S( 11-5(1-1)] 

CALL SUM(2, L,B11,SS,ANS) 

K11=K11+.5*ANS 

CALL SUM(2,L,B12,SS,AN$) 

K12=K12+.5*ANS 

CALL SUM(2t L, 822 , SS , ANS) 

K22=K22+. 5*AN$ 

CALL SUM(2,L',B66,SS,ANS). 

K66=K66+.5*ANS 

CALL SUM(2, L,B11,SS.S,ANS) 

D 11=011+ ANS /3. 

CALL SUM (2 , L , B12, SSS » ANS )' 

012=D12+ANS/3. 

CALL SUM(2, L,B22,SSS,ANS) 

D22=D2-2+ANS/3. 

CALL SUM( 2 , L» B66, SSS, ANS) 

D66=D66+ANS/3. 

80. IF (ST-1) 51 ,51, 52 
52 READ STIFF 
PRINT STIFF 

D11=D11+ES*IS/8S+ES*AS*ZS**2/BS 
022=0 22+ ER* IR'/RS+ER*AR*ZR**2'/RS 
K11=K11+ES*AS*ZS/BS 
K22=K22+ER*AR*ZR/RS- 
C11=C1T+ES* AS/ BS 
C22=C22+ER* AR/RS 


51 

PRINT 

13 






DO 53 

L=li 

.LAYER 



5-3 

PRINT 

14., 

L, Bill L) 

, B12I 

CL) , B22 ( L) , B66tL) 


PRINT 

15, 

Cll, 

Cl 2, 

C22.t 

C66 


PRINT 

16, 

Dll, 

D1'2, 

D22, 

D66 


PRINT 

17, 

Kll , 

K12, 

K22-, 

K66 


C PREBUCKLING QUANTITIES 

RG1= (-C11*C22-C 12**2)/ (01 1*C11-K11**2I 
RG2=(K11*C12-K12*C11)./{D11*C11-K11**2) 
RG3=(R*C11) /(2.*D11*C11-2.*K11**2) 

RG8= ( SQRT( RG1 ) =RG2-NXB*RG3 ) / ( 2 . *R ) 

RGG= ( SQRT( RG1 ) +RG2+NXB*RG3) /( 2«*R> 
BETA=SQR,T( RGB) 

GAMA=SQRT( RGG) 

BET.AL=-C11/ (8.*BETA*,( D1 1*C1‘1-K11**2-) ) 
GAMAL=C11/ ('8.*GAMA*( D1 1*C 11-K11**2 ) ). 



113 


Pl= ( R*C12) / ( Cl 1*C22-C12**2) 

SS=SI NH (BET A*A/2 » )*S IN! GAMA*A/2 . > 

SC=SINH(BETA*A/2. )*C0S!GAMA*A/2. 1 
CC=COSH{ BET A*A/2. )*C0S!-GAMA*A/2. 1 
CS=COSH (BET A* A/2 , ) *S I N! GAMA*A/2. ) 

SHCH=SINH(BET A*A/2. ) *CQSH ( BETA* A/ 2. ) 

SNCN=SIN!GAMA*A/2.)*C0S (GAMA*A/2. ) 

GO TO!l,2,2),BC 

1 G1=NXB*P1*( BETA**2-GAMA**2)-(K11*NXB)/!DU*C11~K11**2) 
G2=2.*NXB*P1*BETA*GAMA 

DBAR=D11*CI1-K11**2 

G1L=P1*!BETA**2-GAMA**2+2.*NXB=MBETA*BETAL-GAMA*GAMAL> J-Kll/DBAR 

G2L=P1*(2. *NXB*( BET A*GAMAL+GAMA*BET AL)+2«*BETA*GAMA) 

G12=G1*CC-G2*$S 

G21=G1*SS+G2*CC 

G=2.*BETA*GAMA*!SS**2+CC**2) 

F1=G12/G 

F2=-G21/G 

PRINT 999, BETA, GAMA, FI, F2 

999 FORMAT!'//* BETA=*E16.8, 5X*G AMA=*E16. 8 , 5X*F1=*E16. 8, 5X*P2=*E16. 8 ) 
FOl=2.*G12* (8ETA*GAMAL+GAMA*BETAL )* ( SS**2+CC**2 )/G**2+2. *G12*GAMA* 
1BETA*A*(BETAL*SHCH-GAMAL*SNCN)/G**2 
F1L=G1L*CC/G-G2L*SS/G+A*SC*(G1*BETAL-G2*GAMAL) / !2.*GI-A*CS*! G1*GAM 
lAL+G2*B£TALl/( 2.*G)-F0L 

F2L=.-G1U*SS/G-G2L*CC/G+A*CS*( G2*GAMAL^G1*BETAL J /! 2.*G )-A*SC*(Gl*GA 
1MAI+G2*BETAL)/!2.*G)+F0L*G21/G12 
GO TO 3 

2 G=BETA*SNCN+GAMA*SHCH 
F1=NXB#P1* ( BET A*SC— GAMA*CS’I /G 
F2=— NXB#P1 * ( B£TA*C$+GAM A#SC )/G 
G12=F1*G 

G21=— F2*G 

FGL=NXB#P1#G12*( BET AL*S IN! GAM A* A )-+A*BET A*GAMAL*COS ( GAMA*A)+GAMAL?( 
I-EXP! BETA*A > -EXP (—BET A*A ) ) /2.+A*GAMA*BETAL*! EXP ! BETA*A) + EXP (— BET A* A 
2)1/2. )/! 2,*G**2) 

F1L=P1*G12/G+NXB*P1*!BETAL*SC-GAMAL*CS)/G+A*NXB*P1*CC*( BETA*BETAL- 
1G 

1AMA*GAMAL)/ ( 2. *G)-A*NXB*P1*SS* ( BETA*GAMAL+GAMA*BETAL> /( 2.*G )-FOL 
F2L=-P1*G21/G-NXB*P-1*!BETAU*CS+GAMAL*SC)/G-A*NXB*P1*SS*!BETA*BETAL 
1- 

1GAMA*GAMAL) /( 2 .*G )- A*NX B*P1*CC* ( BETA*GAMAL+GAMA*BETAL ) /'< 2.*G) + FOL* 
2G21/G12 

3 PRINT 18 
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DO 4 1=1, NSTAT 
J=I-1 

IF( I— 1) 5,5, 6 

5 XB=— A/2 • 

GO TO 7 

6 X=FLOAT(J)*A/FLOAT(NSTAT-l) 

XB=X- A/2. 

7 SXSX=SI NH( BET A*X8)*SIN<G AMA*XB ) 

CXCX=COSHLBETA*XB)*COS(GAMA*XB) 

'SXCX= SINH( BETA*XB)*CaStGAMA*XBJ 
CXSX=COSH( BETA*XB)*SIN(GAMA*XB) 

WO( I )=F1*SXSX+F2*CXCX+P1*NXB 

WOXd ) = <F1*BETA-F2*GAMA)*CXSX+(‘F1*GAMA+F2*BETA)*SXCX 
WOXXd > = (F1*(BETA**2-GAMA**2)-2.*GAMA*BETA*F21*SXSX+(2.*6AMA*BETA* 
lFl+( BETA**2-GAMA**21*F2 )*CXCX 
F1LL=F1L*(BETA*CXSX+GAMA*SXCX)+F2L*<BETA*SXCX-GAMA*CXSX> 
B1LL=BETAL* (.F1*CXSX+F2*SXCX+XB*{ F1*BETA— F2*GAMA )*SXSX+XB*(F1*GAMA+ 
1F2*BETA)*CXCX) 

G1LL=GAMAL*(F1*SXCX-F2*CXSX+XB*(F1*BETA-F2*GAMA )*CXCX-XB*(F1*GAMA+ 
1F2*BETA)*SX$X) 

WOXL ( 1 1 =F1L L+B 1LL+G1LL 

WOLd )=F1L*SXSX+F2L*CXCX+XB*!F1*BETAL-F2*GAMAL) *CXSX+XB*( F1*GAMAL 
1+F2*BETAL) *SXCX+P1 
IF( IPRINT.EQ.O ) GO TO 4 

PRINT 502,. I, WOd), WOXd), WOXXd), WOXLd), WOLd) 

4 CONTINUE 
T=0 . 

DO 55 L=l, LAYER 
55 T=H(L)+T 
M=l. 

C FORM L.H.S. OF EQUILIBRIUM EQUATIONS A*Z d-1 ) +B*Z d) + C*Z ( 1+1 ) 

250 CONTINUE 

DEL=A/FLOAT (NSTAT-1 ) 

L AM 11=0 11— K 11**2/ D1 1 

LAM12=— C66* (M*N/R)**2 

LAM 13= { C12+066— K12*K11/D11 ) *M*N/R 

LAM17=K11/D11 

LAM21--(C12+C66)*M*N/R 

LAM22=C66 

LAM23=-C22*(M*N/R)**2 
LAM24={ K12+2.*K66 )*M*N/ R 
LAM32=— 2.*K-66*M*N/R 
LAM41=-K11 



LAM42=-K.12*M*N/R 

LAM43=D11 

LAM45=-j-( D12*( M*N/R) **2+K12/R) 

NEND=NSTAT+1 
IF(ST.NE.O) GO TO 99 
GS JS=0* 

GRJR=0« 

BS=1. 

RS=1 • 

99 DO. 100 1=1, NS TAT 

LAM14=<C11-K11**2/D11>*W0X( I ). 

LAM15=(C11~K11**2/D11)*W0XX(I )+ ( C12-K12*K'11/D1 1) /R+ (K12+2 . *K66-K11 
1#012/D11 J*( M*N/R)**2 
LAM16=-O66*W0X{I1*( H*N/R)**2 
LAM25=— ( Cl 2+C'66)*W0X( I ) *M*N/P 

LAM2‘6=-C66*W0XX( I >*M*N/ R-C22*M*N/ R**2-K22*( M*N/ R) **3 
LAM31=C12/R +{ K12+2. *K6$ J *( M*N/R >**2-Cll*W0XX (II 
LAM33=?C22*M*N/R**2+K22*CM*N/R)**3-C12*W0XX( I.]*M*N/R 
. LAM34=NXB-K12/R-(4.*D66+GSJS/BS+GRJRARS+D12.).*(M*N/R>**2+K11*W0XX(-I 
T) 

LAM35=(C12/R+( K12+2,.*K66 J * ( ^*N/R) **2)*W0X(I )-Cll*>J0X( I J*WOXX( I ) 
LAM36= ( ( C12'*K1 l/Cll+-2.*< K66-K12) J *( M+N/Rl **2-C12/R) *W0XX( I) +1 . /R*{ 
1C 

122-C12**2/C11 )*(M*N/R)**2*W0(I ) +D22* (M*N/R)**4+2.*K22/R* (M*N/ R) ## 
22+C22/R**2-C12*{ M* < N/R)#*2*NX8/C11 
LAM44=-K11*W0X(I )’ 

A A (I ,1,1)=LAM11 

AA(T , 1,2)=-DEL*LA M13/2. 

AA(i;i,3)=LAM14-DEL*LAM15/2. 

AA 1 1 , 1 , 4) =- DE L*L A Ml 7/2 . 

AA( I ,2>l)=-'DEL*TAH21/2. 

AA(.I,.2,2)=LAM22 

A‘A{ I,,2,3)=LAM24— DEL*LAM25/2. 

AA( I « 2»4)=0 « 

AA{ I , 3 , l}=— DEL+LAM31/2. 

AA<I.,3,2) = LAM32 

AA{ 1,3 »3) = LA'M34— DEL*LAM35/2. 

AA( 1 ,3,4)=-l. 

AA{ T, 4, 1)=— DEL*LAM41/2. 

AA(I,4.,2) = 0. 

AAl I » 4,3)=L AM43— DEL+LAM44/2 . 

A A (.1,4,41=0. 

B(I,1 ,l-)=r-2.*LAMll+DEL**2*LAM12 
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8(1 

,1 

,21=0. 

B ( X 

, 1 

,3}=-2.*LAMl'4+DEL4*2*LAM16 

B ( I 

,1 

,41=0. 

B( I 

1 2 

,11=0. 

B ( I 

,2 

,21=-2.*LAM22+DEL**2*LAM23 

B( I 

,2 

,3 )=-2 .*LAM24+DEL**2*L AM26 

B(I 

,2 

,41=0. 

B ( I 

,3 

,11=0. 

B( I 

,3 

, 21=5-2. *LAH32+DEL**2*LAM33 

B ( I 

,3 

,3)=-2.*LAM34+DEL**2*LAM36 

B( I 

,3 

,41=2. 

B ( I 

,4 

,11=0. 

B( I 

,4 

,2>=LAM42*DEL**2 

8(1 

,4 

,3)=-2.*LAM43+DEL**2*LAM45 

B( I 

,4 

,4)=DEL**2 

C(I 

,1 

,1 )=LAM11 

C(I 

,1 

» 2) =DE L*LAM13/2. 

C ( I 

,1 

,3 )=LAH14+DEL*LAM15/2. 

C(I 

,1 

,4)=DEL*LAM17/2. 

C(I 

,2 

, 1 )=DEL*LAM21/2. 

C<1 

,2 

, 2 1 =LAM22 

C ( 1 

,2 

,31=LAM24+DEL*LAM25/2. 

C(I 

,2 

,41=0. 

C( I 

,3 

,11 =DE L*LAM31/2. 

C(I 

,3 

, 2)=LAM32 

C(I- 

,3 

,3)=LAM34+DEL*LAM35/2. 

C(I 

,3 

, 4 1 =- 1 . 

Ctl 

,4 

, 1 )=DEL*LAM41/2. 

C ( I 

,4 

,21=0. 

C(I 

,'4 

,3 l=LAM43+DEL*LAM44/2. 

100 C ( I 

,4 

,41=0. 

GO 

TOC 101 ,102 , 103 1 , BC 


101 AB(l,l,l)=-Cll*0EL/2. 

AB( 1,1,21=0. 

AB (1,1, 3 )=-Kl 1-C1 l*DEL*WOX ( 1 1 /2. 

AB ( 1 , 1 , 41=0 • 

BB ( 1-, 1 ,.l ) = 0.« 

BB(1,1,21=C 12*M#N*DEL**2/R • ' 

BBU,1,3) = C12*DEL**2/R+2.*K'11+K12=K>EL**2=KM*N/R1**2 
BB ( 1,1,41=0* 

BB( 1,2,11=0. 

68(1,2,21=1.0 

BBC 1,2,31 = 0. 



BB( 1 , 2 , 4)=0 . 

BB(1»3»1)=0. 

BB(1,3‘,2)=0. 

BB( 1, 3,3)=1 .0 
BB( It 3,4)=0. 

BB( 1,4,1)=0'. 

6B(lt4,2)=0. 

BB(1 » 4,3 )=0 . 

BB(1,4,4>=1 .0 
CB(l,l,l)=Cll*DEL/2. 

CBI.lt 1, 2)=0 . 

CB(1,1,3)=-'K1 1+C11*0EL*W0X ( 1 ) /2 . 

CBI 1, 1,4)=Q.. 

00- 110 J=2 , 4 
DO 110 K=1 ,4 
AB( 1, J,K)=0 * 

110 CB( It J ,K)=0 . 

DO 111 J=l,4 
DO 111 K=1 ", 4 

C STATION 2 CORR ESPCNOS TO THE NSTAT BOUNDARY IN AB BB AND CB 
ABI2, J,K}=ABQ ,J,K) 

BB(2t.J,K) = BB(l ,J,K). 

CB.(2, J,K)=CB(l,JtK) 

AB(.2,l,3}=-Kll-Cll*DEL*W0X(NSTAT}/2. 

111 CB(2,l,3I=-Kll+CU=f'DEL«W0X(NSTAT)/2. 

IFIM.GT.O. ) GO TO 105 

DO. 510 J=1 ,4 
AB{2»1, J)=0 . 

BB( 2 ,'l , J )=0 . 

510 CB(2,1,JI=0. 

BB(2,1,1»=1 . 

GO TO 105 
102 DQ 112 J=1 , 4 
DO 112- K*r, 4 
AB(1,U,K)=0 . 

BB ( 1 , J , K)=Q « 

CBC1 ,'J, K’)=0 • 

AB1.2 , J t K) =0. 

BB ( 2 , J ,.K) = 0» 

112 CB( 2 ,J , K) =0. 

AB<l,4,3)=-DEL/2. 

ABI 2, 4,3)=- DEL/2. 

BB< l.,l,l) = l . 
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BB( 1, 2, 2)=1 « 
BB(1,3»3)=1 . 

BB (2 ,1,11 = 1. 

BB( 2 ,2,2)=1. 

BB ( 2 ,3,3)=1. 

C8(l,4,3) = DEL/2., 
CB(2,4,3>=DEL/2. 
GO TO 1.05 
103 DO 113 J=1 f4 
DO 113 K=1 1 4 
AB ( 1 » J» K )=0 . 
BB(1, J,K)=0, 
CB(1,J,K)=0. 
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511 


AB ( 2 ,J,K)^0. 

BB ( 2 ,J,K)=0. 

CB( 2 ,J,KJ=0. 

ABCl,l,l)=-Cll*DEL/2. 

AB(l,l»3)=-Kll-Cll*DEL*W0X(l)/2. 

BBU,1,2>=C12*M*N*DEL**2/R 

BB( 1» 1»3}=C12*DEL**2/R+2.*K11+K12*D£L**2*(M*N/R)**2 
CB(l,l,l>=Cll*DEL/2. 

CB ( 1 » 1 ,3)=-Kl 1+C11*DEL*W0X( 1 ) /-2. 

AB<l,4,3)=-DEL/2. 

AB(2,1»1)=AB(1»1,1) 

AB<2,l,3J=-KU-Cll*DEL*W0X(NSTA»T,)/2. 
BB(2,1,2)=BB<1,1»2) • 

BB{2,1,3)=BB(1,1,3)' 

CB(2» It U=CB( 1,1,1) 

CS(2»l,3)=-Kll+Cll*DEL*W0)^NSTAT)/2. 

AB( 2,4,3) =-DEL/2. 

BB( 1,2,2)=1 . 

BB(1»3,3)=1. 

BB( 2 ,2 , 2)=1. 

BB{ 2 ,3,3J=1. 

CB(l.,4,3)=DEL/2. 

CB(2,4,3)=DEL/2, 

IFIM.GT.O. IGO TO 105 
DO 511 J=l,4 
AB( 2, 1,J }=0'« 

BB( 2, 1, J 1=0 • 

CB{ 2, 1, J)=0 • 

BB(2,1,1)=1. 

GO TO 1C5 



no no no no no on 
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C FORM Bl, Cl, AN, BN AT BOUNDARIES 

105 DO 120 J=1 , 4 
00 120 K=1 , 4 
DUMA< J,K)=AA(1,J,K) 

DUMB! J,Kl=8(l , J,K) 

DUMCt J ,K)=C ( 1 , J, K) 

DUMABI J,Kt=AB( 1,J,K » 

DUMBB ( J » K) = BB( 1,J,K) 

.120 DUMCB( J,KJ=CB( 1, J,K1 
INVERT=10 
MULTIPL =20 
SUBTRAC =22 

CALL MATRIX ! I N VERT , 4 ,4, 0, DUMA, 4, DETERM1 ). 

-1 

DUMA=A1 

CAUL MATRIX (MULTIPL ,,4,4,4, DUMAB,4, DUMA, 4, DUMA, 4)' 
1 

DUMA=AB1 * A1 

CALL MATRIX (MULTI PL ,4,4,4, DUMA, 4, DUMB, 4, DUMB, 4) 
-1 

DUMB = AB-1 * A1 * B1 

CALL MATRIX (SUBTRAC, 4, 4,0, DUMBB, 4, DUMB, 4,-Bl, 4) 

» '• , 

B1=BB1-*AB1*A1 *B1 ' 

CALL MATRIX (MULTIPL ,4,4., 4, DUMA, 4, DUMC, 4, DUMC, 4 ) 
-1 

'DUMC=AB1*A1 *C1 

CALL MATRIX (SUBTRAC ,4, 4,0,.DUMCB, 4, DUMC,4,C1 ,41 ‘ 
-1 

C1=C'B1-AB1*A1 *C1- • , , 

DO 121 J=1 , 4 
DO 121' K=1 »>4 

DUMA( J »K1=A A( NSTAT, J, K) , , 

DUMB ( J ,K) = 8 (NSTAT,, J ,K) * . , 

DUMC{ J,K)=C (NSTAT,J,K) 

DUMAB( J,KI=AB(2 ,J,K! ‘ 

DUMBB (J,KJ=8B(2 ,J,K) 

121 DUMCB { J,K>=CB( 2 ,,J,K) 

CALL MATRI X ( IN VERT,4,4,0, DUMC, 4,DETERM1 ) 

CALL MATRIX (MULTI PL .,4,4,4, DUMCB, 4, DUMC, 4, DUMC, 4) 
CALL MATRIX (MULTI PL ,4,4,4, DUMC, 4, DUMA, 4, DUMA, 4) 
CALL MATRIX (SUBTRAC ,4,4,0, DUMAB, 4,.DUMA, 4, AN, 4) 

-1 


C 
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C AN=ABN-CBN*CN *AN 

CALL MATRIX (MULTIPL ,4,4,4, DUMC,4,DUMB,4,DUMB^4) 

-CALL MATRI X ( SU8TRAC ,4,4,0*DUMBB,4» DUMB ,4, BN', 4) 

^■1 

8N=BBN-CBN*CN *BN 
DO 122 J=1 » 4 
DO 122 K=l,4 
D(1,J,K1-BUJ,K) 

CU, J,K)=CHJ,K) 

AAt‘NSTAT,J,K)=AN(J,K) 

122 B(NSTAT,J,K)=;BN!J,K) 

TEST M AND ROUTE TO 1ST OR 2ND -PERTURBATION EQUATIONS 
300 MT=IFIX(M) 

IF (MT— 1)2 10, 2X1, 21 2 
211 CONTINUE 

IF (TEST3.GT.0. )GG TO 4C7 

TRI ANGULATICN OF A B AND C MATRICES FOR BUCKLING EQUATIONS 
FIRST PERTURBATION 

CALL MATRI X( INVERT, .4,4, 0,81, 4, MODRES) 

Bl=( B1HNV 
M0DRES=DET(:B1) 

CALL MATRIX! MULTIPL ,4,4,4, B1.,4,C1,4,.PIM1,4) 

P(I-1)=B1*C1 
DO - 130 J=1 ,4 
DO 130 K=l,4 
P( 1 , J ,.K }=P I Ml { J , K ) 

130 Rill, J,K)=B(1, J,K) 

SCALE=ABS( MODRES) 

MOOR ES= MODRES/ SC ALE 
DO 13 5 1-2,, NSTAT 
DO 136 J=1 ,4 
DO 136 K=l.,4 
PIMK J,K)=P(I-1,J,K) 

DUMA! J,K)=AA( I ,J,K] 

DUMB! J,K) = B(I, J,K) 

136 DUMC ( J ,K)=C!T , J,K) 

CALL MATRIX! MULTIPL ,4,4, 4, DUMA ,4,P IM1 ,4 ,DUMA, 4) 

CALL MATRIX ! SUBTRAC ,4,4,0, DUMB, 4, DUMA, 4, RI , 4) 

C R! I ) = 8-A*P (.1-1 ) 

DO 137 J,=l , 4 
DO 137 K=1 , 4 

137 R1II, J,K)=RI! J ,K) 
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CALL MATRIX (INVERT, 4,4,0, RI-, 4, DETRI ) 

CALL MATRIXIMULTIPL ,'4,4,4, RI ,4, DUMC,4,PI ,4) 

DO 138 J=1 , 4 
DO 138 K=1 , 4 
138 P(I, J»K)=PI (J,K) 

SCALE=A8S( DETR I ) 

IFt I.EQ.NSTAT) SCALE=l. 
detri=detri/scale 
135 MODRES=MODRES*DETRI 
PRINT 20, MQDRES 
C EIGENVECTOR CALCULATION 

DO 160 J=l , 4 
DO 160 K=1 , 4 

160 DUMA( J ,K)=R1 ( NST AT, J ,K) 

U1 l NSTAT) = DUMA (1,2) #DUMA( 2,3)*DUMA( 4»4)+DUMA( 1 » 3 )*DUMA( 2,41 *DUN 
1,2) +DUMA( 1 , 4)*DUMA( 2, 2 ) *DUMA( 4, 31— DUMA ( 4,2) #DUMA( 2 , S')* DUMA ( 1,4) 
2MA(4,3)*DUMA(2,4)*DUMA(1,2)-DUMA(4,4)*DUMA( 2, 2 ) *DUMA( 1 , 3) 

VI ( NSTAT )=DUMA'(4,1)*DUMA( 2, 3 )*DUMA( 1 ,4 )+DUMA( 4, 3>*DUMA( 2,4) *D^ 

1.1) + DUMA(4, 4) *DUMA( 2, 1 ) *DUMA ( 1 ,3 ]-DUMA( 1, 1 )*DUM AT 2 , 3 >*DUMA( 4, 4 ) 
2MA( 1,3)*DUMA(2«4>*DUMA(4,1)-DUMA(1,4)*DUMA(2,1)*D'UMA(4,3) 

Ml ( NSTAT) = DUMA( 1 , 1 ) *DUMA( 2,2)#DUMA( 4,4 )+DUMAl 1,2 )#DUMA( 2,4)=i‘DUf' 
1»1)+DUMA(1»4)*DUMA( 2, 1 ) *DUMA( 4, 2 )-DUMA( 4, 1 ) *DUMA( 2 , 2)*DUMA( 1,4 ) 
2MA(4,2)*DUMA(2»4)*DUMA( 1,1 )-DUMA( 4,4 )*DUMA( 2, 1 ) *DUMA( 1,2) 
M1(NSTAT)=DUMA(4,1)*DUMA( 2, 2 )*DUMA( 1 ,3 )+DUMA( 4, 2)*DUMA( 2,31*0^ 

1.1) +0UMAI4, 3 ) *DUMA( 2,1 )*OUMA U, 2 )-DUMA< 1, 1 )*DUMA( 2, 2 >*DUMA{ 4, 3 1 
2MA(1,2)*DUMA( 2 ,3 ) *DUMA( 4,1 )-DUMA( 1 , 3 )*DUMA{ 2, 1)*DUMA< 4, 2) 

Z1 ( NSTAT, 1 ) =U1 (NSTAT ) 

Z1 (NSTAT,2 ) =V1 ( NSTAT) 

ZKNSTAT ,3) =W1 (NSTAT) 

ZHNSTAT, 41 =M1 (NSTAT) 

WMAX=W1( NSTAT) 

NS1=NSTAT-1 

DO 170 11=1, NS1 

I=NSTAT— II 

DO 171 J=1 , 4 

DO 171 K=1 , 4 

PD( J,K)=-P( I, J,K) - 

171 Z(K)=Z1( 1+1 ,K) 

CALL MATRIX (MULTIPL ,4, 4, 1 , PD ,4, Z, 4, ZD ,4) 

DO 172 K=l,4 

172 Z1 ( I ,K)=ZD( K) 

C Z1(I)=-P(I)*ZI(I+1) • 

UK I )=21(L, 1) 
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VI ( I )=Z1( 1*2) 

Wl( I ) =Z1( 1,3) 

IFIABSIWKIH .GT.WMAX) WMAX=ABS( W1 ( I ) ) 

170 Ml( I )=Z1( I » 4) 

PRINT 23 

23 FORMAT! *1*,4X*STATI ON*, 9X*U1*,17X*V1*, 17X*W1*I 
DO 173 1=1 , NST AT 

Wl( I J=W1(I ) /WMAX 

24 FORMAT ( 18, 3 E20 .8 ) 

173 PRINT 24, It U1 ( I) ,V1 ( 1 1 , W1 ( I ) 

C NORMALIZE Z VECTOR BY ITS MAGNITUDE 

PRINT 30 

30 FORMAT <///4X*STAT ION*, 9X*U1*, 17X*V1*, 17X*W1*, 17X*M1* > 
ZMAG2=0. 

DO 180 1=1 , NST AT 
WKI)=WMAX#W1U) 

180 ZMAG2=U1 (I )**2+VHI)**2-f-WK I) **2+Ml (I) **2+ZMAG2 
ZNORM=SQRT { ZMAG2 ) 

IF(NORM.GT.O) ZNORM=WMAX/T 

DO 181 1=1, NST AT 

Ul( I )=U1(I > /ZNCRM 

VI ( I )=V1 II ) /ZNGRM 

wim=wim/zNORM 

Ml { I )=M1 ( I ) /ZNORM 

Zl( I,1J=U1{ I) > 

Z1 ( I , 2) =V1< I > 

Zl< I,3)=Wll I) 

Zl( 1 , 4)=M1' { I) 

IF( IPRINT, EQ.O ) GO TO 181 

print 31, i, ukij, vim, mm. Him 

31 FORMAT ( I 8* 4E2Q,8 ) 

181 CONTINUE 

C CALCULATE W1X, W1XX, W1XXX,W0XXL, WOXXXL 

0EL=A /FLOAT (NSTAT-1 ) 

X(1)=0 

00 200 1=2 , NSTAT 
200 Xm = X(I-ll+DEL 
PRINT 921 

921 FORMAT ( /// /4X*$T AT I ON*, 9X*W1X*,17X*01X#, 17X#V1X*, 16X*W1 XX* ) 
DO 201 I =1'» NST- AT 
WlX(n=DIF( 1,1, NSTAT, X,W1) . 

Ulxm=DIPI 1,1, NSTAT, X,U1) 

V1X(I)=DIF( I, 1, NSTAT, X, VI) 
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W1XXCI }=K11*(U1X{ I)+WOX( I)*WlX{in/Dll+K12*(N*Vim/R+Wl(I) /R)/ 
1+D12*(N/R>**2*W1( I)/D11-M1( I) /Dll 
W02L ( I )=DI F ( I » 1» NST AT,X»WOXL) 

NY0Lm = tC12*Kll/Cll-K12)*W02L ( I )+ ( C22-C12**2/C11 ) *WOL< 1} /R-Cl 
1C11 

IF t IPRINT. EQ.OI GG TO 201 

PRINT 922, I. WIX(I), Uixm, V1X(I), WIXX(I) 

922 FORMAT ( I8,4E2G .8 ) 

201 CONTINUE 
PRINT 923 

923 F0RMAT(////4X*STATI0N*,9X*NX1*,17X*NY1*,17X*NXY1*, 17X+NY0L*) 

00 202 I=1,NSTAT 

W13X ( I ) =01 F ( I ,1,NSTAT,X,W1XX) 

W03L(I)=DIF(I, 1,NSTAT,X,W02L) 

C CALCULATE NX1,NY1,NXY1 

EPsx=uixm+wox( mwixm 

EPSY=N#V1 ( I )/R+WKI)/R 

EPSXY=-N*UU I ) /R+Vl X{ I I-WOX { I )^N*W1 (I) /R 

NX1 (I >=Cll*EPSX+C12*EPSY-Kll*WlXXm + (N/R)**2*Wl( I)*K1'2 

NY1 { I )=C12#EPSX+C22#EPSY*-K12*W1XX( I ) + ( N/R)**2*W1( I ) #K22 

NXY1 ( I J =C66#EPSXY+2.#K66#N*W1X{ I ) /R 

IF( IPRINT. EQ.O) GO TO 202 

PRI NT "924, I, NXKIlt NYKI), NXYl(I) »NY0LU> 

924 FORMAT! 18»4E20. 8) 

202 CONTINUE 

M=0. 

GO TO 250 
210 CONTINUE 

C SECOND PERTURBATION TO OBTAIN U20,V20,W20 

225 DO 226 1=1, NST AT 

D C I ■» 1 )=-Cll*WlX( I )*W1XX< I)*DEL**2/2.-(DEL*N/RJ**2*C12*Wim*Wr)i 

1 / 2 . 

D( I » 2 )=C . 

DL2=DEL«*2 

E1=DL2*K11#(W1XX( I) ##2 +W1X (I)*W13X ( I) ) /2. 

E2=DL2 * (N/ R>**2* (K12/2 . ) #W1( I ) #W1XX( I ) 

E3=0L2=! < ( K12*( N/R) **2/2.+ ( Cll+WOXX ( I )— C12/RL/4. )*W1X(I)**2 
E4=DL2*( N/R )**2/4.*( C12*W0XX( I J-C22/R )*W1< I)**2 
E5=DL2*NX1 ( I ) *W1XX( I ) /2.^N/R*NXY1 ( I )*W1X< I 1 *DL2-< N/R )**2*DL2/2 . 
11(I)*W1(II 
D{ I ,3)=E1+E2+E3+E4+E5 

226 D( I » 4 )=0 • 

GO TO ( 220 * 221, 221), BC 
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C FORM 01 AND DN AT BOUNDARIES 

220 DBtl )=?-DEL**2*Cll*WlX(l)**2/4. 

DB( 2.) =0 • 

0B( 3)=C. 

08(4) =-DEL**2*Kll*WlX (-1 >**2/4. 

00 700 J=1 , 4 
DO 700 K=1 1 4 
DUMA (J,K)=AA(1,J,K) 

OUMC'( J,,K)=C (NS TAT » J »K) 

OUMAB( J,K)=AB( 1,J,K) 

700 DUMCBU,K)=CB( 2, J,K) 

CALL MATRIX! INVERT ,4, 4,0, DUMA, 4, DETERM I > 

CALL MATRIX { INVERT, 4 ,4,0 , DUMC,4, DETERMI ) 

CALL MATRIX! MULTI PL ,4,4,4, DUMAB, 4, DUMA, 4, DUMA, 4 > 

CALL MATRIX (MULTIPL ,4,4,4,DUNCB,4,DUMC,4,DUMC,4) 

DO 701 J=1 » 4 
7.01 DUMD( J 1=D( 1 , J > 

CALL MAT.RI X ( MULT IPL ,4,4, 1 , DUMA, 4, DUMD, 4 , DUMD, 4) 

CALL M A T RI X ,( SU BTRAC , 4 , 1 , 0 , D B , 4 , DU MD , 4 , DUMD , 4 ) 

‘ DO 702. J=1 » 4 

702 D(l,d >=bUMD(J) 

DB(1)=0. 

DB( 2 ) =0. 

DB(3)=0. 

■'0B(4)=^pEL**2*KlT*WTX(NSTAT)**2/4. 

. 'DO 703 J*l,4 

703 DUMD ( J !—D( N STAT, J.) ■ 

■CALL MATRIX (MULTIPL, ‘4,4,1, DUMC, 4, DUMD, 4, DUMD, 4> 

CALL MATRIX (SU BTRAC ,4,1 , 0,DB,4,DUMD,4.,DUMD,4) 

•DO 704 J=l,4 

704 D(NS-TAT , JJ=DUMO,( J) ‘ 

• . GO TO ,227 

221' DO 223' J=l, 4 
DM, J )=0. 

223 D(NSTAT, J>=0. 

C ’CALCULATE Z20 BY GAUSSIAN ELIMINATION 

227 ' CALL POTTER < A A, B ,C , D,NSTAT,4 ,SMAT, DMAT, PP, QQ, Z-2Q , 61,4) 

M=2. 

GO TO. .250 
212 DL2=DEL**2. 

C SECOND PERTURBATION TO OBTAIN U22,V-22,W22 

235 DO 2 3’ 6 I=1,NSTAT 

DM., 1 ) =-DL2*Cll*WlX ( I ) *W1XX! I ) /2. +DL2* ( N/R) **2* (C12/.2 . +C66 )*W1( 
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l wix ( i j 

El=DL2*N/(2.*R)*(C12+C66)*WlXtI>**2 
E2=-DL2/2.*(N/R)**3*C22*W1( I)**2 
£3;=DL2*N/(2.*R)*C66*{W1XX( I )*W1( I J) 

D ( 1 , 2 )=E1+E2+E3 

E4=Kll/2.*( W1XXU )**2+WlX( I )*W13X(I ) ) 

E5=(N/R)**2*{K12/2. + 2.*K66'>*Wl( I )*W1XXU) 

E6= ( { C11*W0XX( D-C12/R) /4.~ ( N/R) *42* ( 3.*K12/2 . +2.*K66 ) )*W1X( 1)4 
E7= { ( N/R )*44*K22+(N /R) **2/ 4. *(C22/R-C 12*W0XX ( I ) ) )*W1( 1 1 **2 
E8=NX1(I )*W1XX( I)/2.+N/R*NXYl( I )*W1X( I)-(N/R)**2/2.*NY1( I)*W1( i 
0{ I ,3 )=DL2*(E4-E5+E6+E7+E8) 

236 D( I ,4 ) =0 . 

GO TO (230,231,231) ,BC 
230 DB(1)=-DEL**2*C11*W1X(1 )**2/4. 

D8( 2 ) =C • 

OB(3)=jO. 

DB(4)=-DEL**2*Kll*WlX(l)**2/4. 

DO 800 J=l,4 
DO 800 K=l,4 
DUMA(J,K)=AA(1,J»K) 

DUMC( J,K)=C(NSTAT ,J,K) 

DUMAB ( J ,K) = AB ( 1,J,K) 

800 DUMCB ( J , K) =CB ( 2, J ,K ) 

CALL MATRI X ( INVERT ,4,4,0, DUMA , 4, DETERM l ) 

CALL MATRIX( INVERT, 4, 4,0, DUMC.4, DETERMI ) 

CALL MATRIX (MULT I PL ,4,4, 4,DUMAB, 4, DUMA, 4, DUMA , 4) 

CALL MATRIX (MULTI PL ,4,4,4, DUMCB, 4 , DLJMC* 4, DUMC, 4 ) 

DO 801 v|=l , 4 

801 DUMD( J) =D{ 1 ,J ) 

CALL MATRI X( MULT I PL, 4,4,1, DUMA, 4, DUMD,4, DUMD,4) 

CALL MATRIX ( SUBTRAC ,4;,l,0,DB»4,DUMD»4»DUMD, 4) 

DO 802 J=1‘, 4 

802 D(1,J)=PUMD(J-) ' 

DB( 1 ) ~-DEL**2*Cll*WlX (NSTAT)**2/4. 

DB( 2 ) =0 • 

DB( 3) =0« 

DB(4J=-DEL**2*Kn*WlX(N$TAT)**2/4. 

DO 803 J=l,4 

803 DUMD( J )=D( NSTAT, J ) 

CALL MATRI X ( MULTI PL ,4,4, 1 , DUMC, 4, DUMD,4,DUMD,4 ) 

CALL' MATRI X ( SUBTRAC ,4,1»0»DB»4» DUMD,4,QUMD»4) 

DO '804 J=1 , 4 

804 D (NSTAT, J) =DUM D( J ) 



126 


GO TO 237 

231 DO 232 J=l,4 
D(l, J)=0. 

232 D(NST AT »J)=0. 

C CALCULATE Z22 BY GAUSSIAN ELIMINATION 

237 CALL POTTER ( AA ,B ,C , D ,NSTAT ,4,SMAT, DMAT , PP, QQ, Z22, 61,4) 

C CALCULATE W20X ,W2CX X,W20XXX,W22X, W22XX,W22XXX 

PRINT 25 

25 FORMAT ( //#1 *,4X*STATION*9X*U20*, 17X*V20*, 17X*W20*, 17X*M20* ) 

DO 249 1=1 , NST AT 

U20( I )=Z20( 1,1) 

V20(I)=Z20U,2) 

W20U)=Z20( 1,3 I 
M20( I )=Z20( 1,4) ' 

IF ( IPRI NT. EQ«0 ) GO TO 249 

•PRINT 31,1, U20( I‘) , V20 ( I ) , W20(I), M20(I) 

249 CONTINUE 
PRINT 26 

26 F0RMATL//n*,4X*STATI0N*,9X*U2-2*, 17X*V22*, 1 7X*W22*, 17X*M22* ) 

. DO 899 1=1, NST AT 

U22 ( I ) =Z22 ( 1,1) 

V22(I )=Z22( 1,2) 

W22(I)=Z22C 1,3) 

M22( I )=Z22( 1,4) 

IF< IPRINT. EQ.O) GO TO 899 

PRINT 31,1, U22II), V22 ( I ) , h22(I), M22(I) 

899 CONTINUE 
PRINT 500 

500 F0RMAT(////4X*STATI0N*,9X*W2CX*,17X*W22X*,17X*W20XX*,17X*W22XX* 
DO 251 1=1 , NST AT 
W20X ( I )=DI F(I,1,NSTAT,X,W20) 

U20X( I )=DI F ( I , 1,NST AT, X, U20 ) 

W22XI I )=DI F ( I , 1,NST AT,X , W22) 

U22X( I ) =DI F ( I , 1, NST AT, X ,U22) 

V22X( I)=DIF(I,1,NSTAT,X,V22) 

W20XX ( I )=K1 1* ( U20X( I )+WOX< I)*W20X(I ) ) /D11+K12*W20< I) /(R*D11)-M2 
D/Dll 

W22XX ( I )=K1 1#( U22X( I )+WOX(I )*W22X( I ) ) /Dll+K12#( 2.*N*V22 ( I ) /R+W2 
1 )/R) /Dll+D12*( 2.*N/R) *#2*W22 ( I ) /D11-M22 ( I ) / Dll 
IF( IPRINT. EQ.O ) GO TO 251 

PRINT 924, I, W20X( I ) , W22XH), W20XX(I), W22XX(I) 

251 CONTINUE 

C CALCULATE N X20,NY2Q , NX22,NY22, NXY22 
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PRINT 501 

501 F0RMAX{.///y.4X*STATI0N*,9X*NX20*,15X*‘NY20*,15X*NX22*,15X*NY2 2*, 
115X*NXY22*> 

DO 302- 1=1 j NSTAT 

W20XXXI I )=DTF ( I * 1 ,NSTAT» X, W20XX ) 

W22XXX(I)=DIF( 1,1, NSTAT, X,W22XX) 

EPSX=U20X( I )+WOX( I }.*W20X ( I ) 

EPSY=W20(T } /R 

NX20 ( I ) =C1 l*EPSX+Cl 2#EPSY-K 11*W20XX ( I > 

NY20X I >=Cr2*EPSX+C22*EPSY-fU2*W20XXm 
EPSX=U22X( IT+WOX(I) *W22XU ) 

EPSY=2.*N*V22( I)/R+W22(I)/R 

EPSXY=-2.*N*U22< I)/R+V22X[ I >-WOX (I) *2.*N*W22 (IT /R 

NX22 ( I ) =C11*EPSX+C12*EPSY-K11*W22XX ( I ) + ( 2. *N/R ) **2*W22 ( I)*K12 

NY22 ( IT =C12*EPSX+C22*EPSY-K12*W22XX( I ) + ( 2.*N/R) **2*W22 ( I )*K22 

NXY22 ( I )=C66*EPSXY+2.*K66*2.*N*W22Xm/R 

IF ( IPRINT. EQ.O ) GO TO 302 

PRINT 502, I, NX20(I>,' NY20 (I) , NX22 ( I > , NY22 (I) , NXY22( I ) 

502 FORMAT ( 18 , 5 E20 . 8 ) 

302 CONTINUE 

C CALCULATION OF SENSITIVITY INDEX BY PRESENT THEORY 

DO 400 1=1, NSTAT 

DP3 1 I,H=-C11*CW1XX{ I)*WOXL( I)+k!lX< I)*W02L( I ) ) + ( N/Rl **2*W0XL ( I ) 
1WK I )*C66 

DP3( I , 2 J=N/R*C66*W02L (I )*W1 (.D+N/R* { C66+C12 )*WOXL ( I )*W1X( I ) 

DP3( I,3}="W1XX(I)+NX1U )*W02L<n-(N/R)**2*NY0L( I 1 *W1 ( IT 

i+cn*woxx( i i*woxu i)*wixu)-ci 2 /r*woxl< i.»*wixm 

2+Kll*{W0XL ( I J*W13X( I >+2.*W02L< IT *W1 XX ( I J+W03LI I )*W1X( I > ) 

3-2. *( N/R)**2*K66*W02LC I )*W1 ( I )-< N/R >**2*< 2.*K66+K12 ) *W0XL ( I ) 
4*W1X{ I ) 

DP3( I »4)=0. 

EP3( I ,l)=-Cll*(WlXXm*(W20X(I )+.5*W22X( I )) +W1X (I ) * <W20XX<I)+. 5 
1W22XXII ) ) )+ IN/R) **2*C12* ( W1X ( I ) *W22 (I )+Wl( I >*W22X(.I ) I 
2+(N/R)**2*C66*I3.*WlX( I ) *W22 (I) +1 . 5 *W1 ( I }*W22X( I)+W1(I)*W20X( I I 
EP3( I ,2)= N/R*<C12+C66)*IW1X(I )*(W20Xm + .5*W22XU> ) ) 
1+N/R*C66*(W1XX ( I)*W22(I)+Wim*(W20XX< I J-.5*W22XX(I) ) ) 
2+<N/R)**3*C22*Wl( I)*W22(I) 

AK1= Kll*( W 13X m*(W20X(Il+.5*W22X( I> ) + WlXX(I)*(2.*W20XXm 
l+W22XXUn + WlXCI)*(W20XXX(I )+.5*W22XXX(I ) I) 
AK2=-{N/R)**2*K12*(W1XX( I)*W22( I )+WlX(I)*(W20X( I ) +13. /2 .*W22X< I 
l+m-d )4W22XX(I) )+9.4(N/R)**3*K22*Wl( I)*W22( I ) 
AK3=-(N/R)**2*K66*(6.*WTXX( I ) *W22 <I)+W1X( I) *{ 2. *W20X( I ) +9.*W22X 
1) 
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2+W1U )*(2.*W20XX( I) +3.*W22XX( I ) ) ) 

AK4= WOXXt I >*(C11*W1XU )*(W2CX( I ) +. 5*W22X( I ) ) 

1~(N/R)**2*C12*W1( n*W22( I>) 

AK5=1. /R*l C22* (N/ R) **2*W1 ( I »*W22 ( I ) -C12*W1X ( I) * ( W20X( I )+. 5*W22X 
2 ) ) 

AK6= W1 XX( I )*(3./8.*Cll*WlX{I)**2+l./8.*(N/R)**2*C12*Wl ( I )**2> 
AK7=-( N/R) **2* toll I)*(3./8.*C12*W1X( I >**2+1. /8.* ( N/R )**2*C22*W1 ( 
1 * 2 ) 

AK8= . 5*(N/R)**2'*C66*W1X( I }**2*W1 ( I ) 

AK9 = NXl(I)*(W20XXm+.5*W22XX(I) |+ ( NX20 { 1 1 +. 5*NX22C I) ) *W1XX 

AK1P=-(N/R)**2*(2.*NY1< I ) *W22( I )+Wl ( I) *< NY20 H) +. 5*NY22 Cl) > ) 
AK11= ( N/R ) *( 2 .*NXY1 ( I ) *K22X( I )+NXY22{ I ) *W1X( I ) ) 

EP3(I»3)=AK 1+AK2+AK3+AK4+AK5+ AK6+AK7+AK8+AK9+AKJ.0+AK11 

400 EP3(I,4)=0. 

00 410 J=1 » 4 
DPB3 ( 1 T J )=0 . 

EPB3 { 1 1 J )=0 . 

DPB3 ( 2 1 J )=0 . 

410 EPB3 ( 2 . d ) = 0 . 

GO TO (401r 402 » 403 } i BC 

401 DPB3( 1,1>=-C11*W1X( 1)*W0XL<1) 

DPB3 < 1 » 4)=— K11*W1 X( 1 ) *WOXL( 1 V 
DPB3(2»1)=’"C11*W1X( NSTAT ) *WG XL { NSTAT ) 

DPB3 ( 2» 4}=— K11*W1X{ NSTAT )*W0XL{ NSTAT) 

EPB3( l,l}=->Cll*WlXm*(W2QX(l ) + .5*W22X(l) ) 

6PB3(1,4)=-K11*W1X( 1)*(W2GX(1)+.5*W22X(1) ) 

EPB31 2,1)=-C11*W1X( NSTAT )*(W2GX (NSTAT )+.5*W22X( NSTAT) ) 
EPB3(2,4)=-K11*W1X{ NSTAT) *( K20X (NSTAT )+. 5*W22X< NSTAT) ) 

GO TO 402 

403 DPB3( 1 » 1 )*--Cll*WlX( 1 )*WOXL( 1) 

EPB3(1,1)=-C11*W1X(1)*(W2QX(1)+.5*W22X(1) ) 

DPB3(2,1)=-C11*W1X( NSTAT i*WDXL (NSTAT) 

EPB3 ( 2 ♦ 1 )=— Cll *W1X( NSTAT)*( h2QX { NSTAT >+. 5*W22X( NSTAT) ) 

00 TO 402 

402 CONTINUE 
TEST3=1 . 

M=l. 

GO TO 250 
407 CONTINUE 

00 480 J=1 » 4 
00 480 K=1 1 4 
DUMA (J»K)=AA(1»J»K) 

OUMC (JjK)s=C(N$TAT*JtK) 



DUMABl J,K)=AB( 1* J,K> 

480 DUMCBI J,K)=CBI2, J,K> 

CALL MATRIX (INVERTS, 4,0, DUMA, 4, DETERMI ) 

CALL MATRIX! IN VERT, 4, 4, O.DUMC, 4, DETERMI ) 

CALL MATRIX ( MULTI PL, 4 ,4, 4, DUMAB, 4, DUMA, 4, DUMA, 4) 
CALL MATRIX (MULTI PL ,4 ,4 , 4, DUPCB , 4 , OUMC , 4 , DUMC, 4 ) 
DO 481 J=1 , 4 
DLIMDI J ) =DP3 < 1 , J I 

481 DB( J ) =DPB3 1 1 , J ) 

CALL MATRIX (MULTIPL ,4 ,4, 1 , DUMA, 4, DUMD,4, DUMD,4) 
CALL MATRIX ( SUBTRAC ,4,1,0,08,4, DUMD, 4,DUMD, 4) 

DO 483 J=1 , 4 

483 DP3 ( 1 , J I =OUMD! J ) 

DO 484 J=1 , 4 
DUMD! J1=EP3(1,J) 

484 DB( J J = EPB3 ( 1,JI 

CALL MATRIX I MULTI PL ,4,4t 1 , 0UMA»4» DUMD,4, DUMD,4 ) 
CALL MATRI X { SUBTRAC ,4,1,0,DB,4,DUMD,4,DUM0,4) 

DO 485 J=1 , 4 

485 EP3(l,Ji=DUMD{ J) 

DO 486 J=1 , 4 

DUMD ( J ) =DP3 (NSTAT,J) 

486 DB(J)=DPB3t2,Jl 

CALL MATRI XI MULT I PL, 4, 4,1, DU MC, 4, DUMD, 4, DUMD, 4) 
CALL MATRI XI SUBTRAC ,4,1 ,0,,DB,4,DUMD,4,DUMD,4) 

DO 487 J=1 , 4 

487 DP3INSTAT, J )=DUMD!J) 

DO 488 J=1 , 4 

DUMD! J ) =EP3 ( NSTAT , J ) 

488 DB ( J )=EPB3 ( 2, J ) 

CALL MAI RI X I MULT IPL,4,4,1,DUMC,4,DUMC,4, DUMD , 4 } 
CALL MATRI XI SUBTRAC, 4,;, 0,08, 4, DUMD, 4, DUMD, 41 
DO 489 J=1 , 4 

489 EP3INSTAT, J )=D.UMD(J) 

DO 420 J=l» NSTAT 

DO 420 K=1 , 4 

I=K+4*(J-1I 

zivin=zi(j,K) 

DUMDEI I )=EP3( J ,K) 

420 CONTINUE 
XTE=0. 

XTD=0. 

NVAR^4*NST AT 
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00 425 I —1 * NVAR 

425 XTE=Z1V( II *DUMDE{ I) +XTE 
DO 430 J=1 1 NST AT 
DO 430 K=l,4 

1 = K.+4*!J-1) 

DUMDE(I)=0P3( J,K) 

430 CONTINUE 

DO 435 1=1, NVAR 
435 XTD=Z1V(I) *DUMDE!I) +XTD 
DO 460 1=1 ? NST AT 

460 FOFXII)=Ull I)**2+V1 ( I!**2+WI( I 1**2 
C 

CALL SIMP! DEL , NSTAT » FOFX, AREA } 

637 FORMAT (2X*AREA( Ul, VI, hl)=*E16,8> 

PRINT 637, AREA 

638 FORMAT! /2X#XTE=*E16« 8 ) 

PRINT 638, 'XTE 

639 FORMAT(/2X*XTD=*E16 • 8) 

PRINT 639, XTD 

BCAR0=-2.*T**2*A*XTE/(NXB*XTD*AREA> 

PRINT 635, BCARD 

635 FORMAT (2X*SENSITIVITY INDEX BY PRESENT THEORY B=*1E2Q.8J 
C CALCULATION OF SENSITIVITY INDEX BY BUDIANSKY THEORY 

DO' 470 1=1, NST AT 

470 FOFX { I )= W1X< I )**2*!NX20! I >+,5*NX22! I ) )+(N/R*Wl 1 1 ) >**2*(NY20( I ) 
1-'.5*NY22(I) )-<N/R)*NXY22<I)*WlX( n*Wl(I) 

2+3./8.*Cll*WlXm**4+l./4,#(N/Ri*42 : « < C12*WlX( I »**2*W1! I J **2 
3+3./8.#(N/R ]**4*C22*W1( I ) **4+. 5*( N/R )**2*C66* ( W IX ( I J*Wll I H**2 
CALL SIMP! DEL, NSTAT, FOFX, AREA1 I 
PRINT 640, AREA1 

640 FORMAT !2X#AREA1=*E16.81 
DO 471 1=1, NSTAT 

471 FOFX( I )= NX1( I )*WlXm*(W20Xm + .54W22X(in 

1- *-(N/R>**2*NYl! I )*W1! I)*W22(n 

2— N/R*NXY1! I )*W1X! I) *W22 < I )— N/R*NXY1 ! I )*W1 ! I )*! W20XI I 5*W22X ( I 
CALL SIMP! DEL , NSTAT , FOFX, ARE A2 ) 

PRINT 641, ARE A2 

641 F0RMAT!2X*AREA2=*E16.8) 

DO 472 1=1, NSTAT 

472 FOFX! I I = -W1X< I )**2+ ( N/R >**2*NY0L ( I ) *W1 (II **2 
l+W0XLII)*(2.*NXim*WlX(I )-2.*N/R*NXYl( I)*W1! I) J 

CALL S1MP(DEL, NSTAT, FOFX, AREA3) 

PRINT 642, AREA2 
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642 F0RMAT(2X*AREA3=*E16.8J 

BBUD=-r( AREA 1+2.* ARE A2) / < NXB#AREA3) 

PRINT 636, BBUD 

636 F0RMAT{/2X*SEN$ITIVITY INDEX BY BUDIANSKY THEORY B=*1E20.8J 
GO TO 1000 
END 
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PROGRAM LENGTH INCLUDING I/O BUFFERS 
045532 

FUNCTION ASSIGNMENTS 


STATEMENT ASSIGNMENTS 


1 

- 

000632 

2 


001027 

3 

- 

001207 

4 

- 

001451 

5 

-* 

001220 

6 

- 

001223 

7 

- 

001232 

10 

- 

007334 

11 

- 

007356 

12 

- 

007365 

13 

— 

G07370 

14 

- 

007375 

15 

- 

007401 

16 

- 

007404 

17 

- 

007407 

18 

- 

007412 

19 

- 

007423 

20 


007426 

21 

- 

007432 

22 

- 

007435 

23 

- 

007475 

24 

- 

007503 

25 

- 

007554 

26 

- 

007563 

30 

- 

007506 

31 

- 

007514 

51 

- 

000341 

52 

- 

000272 

55 

- 

001457 

80 

- 

000267 

99 

- 

0G1554 

101 

- 

002163 

102 

- 

002322 

103 

- 

002364 

105 

- 

002507 

137 

- 

003213 

138 

- 

003253 

160 

- 

003312 

170 

- 

003515 

172 

- 

003474 

181 

- 

003633 

200 

— 

003644 

201 

- 

0C3764 

202 

- 

004103 

210 

- 

004107 

211 

- 

003052 

212 

- 

0C4517 

220 

- 

004237 

221 

- 

004467 

225 

- 

004107 

227 

- 

004500 

230 

- 

004704 

231 

- 

005137 

235 

- 

004521 

237 

- 

005150 

249 

- 

0Q522Q 

250 

- 

001464 

251 

- 

005402 

300 

- 

003045 

302 

- 

005553 

401 


006216 

402 

- 

006276 

403 

- 

006255 

407 

- 

006301 

425 

■** 

006647 

435 

- 

006674 

460 

- 

006702 

483 


006434 

485 

- 

006500 

487 

— 

006545 

485 

- 

006612 

500 

- 

007572 

501 

- 

007601 

502 

- 

007612 

635 


007642 

636 

- 

007662 

637 

- 

007627 

638 

- 

007634 

639 

- 

007637 

640 

- 

007651 

641 

- 

007654 

642 

- 

00 76 57 

701 

- 

0C435O 

702 

- 

004404 

703 


004424 

704 

— 

004460 

801 

- 

0C5015 

802 

- 

005051 

803 

- 

00 5074 

804 

- 

005130 

899 

— 

005255 

921 

- 

007517 

922 

10C-0 


0C7526 

000003 

923 


007531 

924 


007540 

999 


007451 


BLOCK NAMES AND LENGTHS 
VARIABLE ASSIGNMENTS 


A 

- 041304 

AA 

- 011401 

AB 

- 017161 

AK1 

- 041444 

AK10 

- 041455 

AK11 

- 041456 

AK2 

- 041445 

AK3 

- 041446 

AK4 

- 041447 

AK5 

- 041450 

AK6 

- 041451 

AK7 

- 041452 

AK8 

- 041453 

AK9 

- 0414 54 

AN 

- 017521 

ANS 

- 041342 

AR 

- 041314 

AREA 

- 041460 

ARE A1 

- 041462 

AREA2 

- 0 4-1463 

AREA3 

- 041464 

AS 

- 041321 

B 

- 013321 

BB 

- 017221 
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BBUD 

-* 

041465 

BC 

- 

041245 

BOARD 

- 

041461 

BETA 

- 

041350 

BETAL 

- 

041352 

BN 

- 

017541 

BS 

- 

041324 

B1 

-T 

017461 

BILL 

- 

041406 

811 

- . 

010773 

812 


011017 

B22 

- 

011043 

B66 

- 

011163 

C 

- 

015241 

CB 


017261 

CC 

- 

041356 

CS 


041357 

CXCX 

- 

041402 

CXSX 

— 

041404 

Cl 

- 

017501 

€11 

- 

041327 

C12 

- 

C41330 

C22 

- 

041331 

C66 

- 

041332 

0 

- 

025741 

•DB 

- 

036411 

DEAR 

-r 

041364 

DEL 

- 

041411 

DELBAR 

- 

041305 

DETERM I- 

041440 

DETERM1- 

041416 

DETRI 

- 

041421 

DL2 

- 

041432 

DMAT 

- 

030245 

DPB3 

- 

041216 

DP3 

- 

036613 

DUMA 

- 

017321 

DUMAB 

- 

017401 

DUMB 

- 

017341 

DUMBB 

- 

017421 

DUMC 

- 

017361 

DUMCB 

- 

017441 

DUMD 

- 

036415 

DUMDE 

- 

037563 

Dll 

- 

041335 

D12 

T" 

041336 

D22 

- 

041337 

D66 

- 

041340 

EPB3 

- 

041226 

EPSX 

- 

041427 

EPSXY 

- 

041431 

EPSY 

- 

041430 

EP3 

- 

037177 

ER 

- 

041313 

ES 

- 

041320 

EX 

TX 

010603 

EY 

- 

010627 

El 

- 

G41433 

E2 

- 

041434 

E3 


041435 

E4 


041436 

E5 


041437 

E6 


041441 

E7 

- 

041442 

E8 

- 

041443 

FOFX 

- 

Q40533 

FOL 

- 

041374 

FI 


041372 

F1L 

- 

041375 

FILL 

-r 

041405 

F2 

-T 

041373 

F2L 

- 

041376 

G 

- 

041371 

GAMA 

- 

C41351 

GAMAL 

• 

041353 

GR JR 


041316 

GSJS 

- 

041323 

GXY 


010723 

G1 

- 

041362 

GIL 

- 

041365 

GILL 


041407 

G12 

- 

041367 

G2 


041363 

G2L 

- 

041366 

G21 


041370 

H 

- 

010747 

I 

- 

041341 

II 

- 

041424 

INVERT 


0414T4 

IPRINT 

- 

041312 

IR 

- 

041236 

IS 

- 

041237 

J 

- 

041377 

K 

- 

041413 

Kll 

- 

041240 

K12 

- 

041241 

K22' 

- 

041242 

K66 

T- 

041243 

L 


041325 

LAM11 

- 

041250 

LAM12 


041251 

LAM13 


041252 

LAM14 

- 

041265 

LAM15 

- 

041266 

LAM16 

- 

041267 

LAM17 


041253 

LAP21 

- 

041254 

LA M2 2 

- 

041255 

LAM23 

- 

041256 

LAM24 

- 

041257 

LAM25 

- 

041270 

LAM26 

- 

041271 

LAM31 

“ 

041272 

LAM32 

- 

04-1260 

LAM33 

- 

041273 

LAM34 

- 

041274 

LAM35 

- 

041275 

LAM36 


C41276 

LAP41 


041261 

LAM42 

- 

041262 

LAM43 

- 

041263 

LAM44 

- 

041277 

LAM45 

- 

041264 

LAYER 

- 

041307 

M 

~ 

041247 

MQDRES 

IT 

041300 

MT 

- 

041417 

MULTIPL- 

041415 

Ml 

“ 

023770 

M20 

- 

034374 

M22 

- 

0347.60 

N 


041306 

NENO 


041412 

NORM 

- 

041311 

NSTAT 

- 

041310 

NS1 


041423 

NUX 

- 

010653 

NUY 


010677 

NVAR 

- 

041457 

NXB 

- 

041244 

NXY1 

- 

025644 

NXY22 

— 

036314 

NX1 

— 

025452 

NX20 

- 

035730 

NX22 

- 

036122 

NYOL 

— r 

011207 

NY1 


025547 

NY 20 

- 

036025 

NY22 

- 

036217 

P 

- 

017561 

PD 

- 

024455 

PI 

- 

023461 

PIM1 

- 

023421 

PP 

- 

030631 

PI 

-T 

041354 

QQ 

- 

032551 

R 

m 

041303 

RGB 


041346 

RGG 


041347 

RG1 

- 

041343 

RG2 

~ 

041344 

RG3 

- 

' 041345 

RI 

- 

023441 

RS 

- 

041317 

R1 

~ 

021501 

S 

- 

011067 

SC 

- 

041355 

SCALE 

-r 

041420 

SHCH 


041360 

SMAT 

- 

026325 

SNCN 

-T 

041361 

ss 


011113 
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SSS 

- 

011137 

ST 

- 

041246 

SUBTRAC- 

041301 

SXCX 

- 

041403 

SXSX 


041401 

SI 

— 

041333 

S2 

T* 

041334 

T 

- 

041410 

TEST3 

- 

041326 

TRANSPO- 

041302 

U1 

- 

023501 

U1X 

- 

024577 

U2C> 

a - 

03410 5 

U20X 

- 

035152 

U22 

- 

034471 

U22X 

- 

035247 

VI 

- 

023576 

V1X 

- 

024674 

V20 

- 

034202 

V22 

- 

034566 

V22X 


035344 

WMAX 

- 

041422 

WO 

- 

010217 

WOL 

- 

011304 

wox 

— 

010314 

WO XL 

- 

C10506 

WOXX 

- 

010411 

W02L 

- 

025163 

•W03L 

- 

025355 

W1 

- 

023673 

W1X 

* 

024771 

W1XX 

- 

025066 

W13X 

- 

025260 

W20 

- 

034277 

W20X 


035055 

W20XX 

- 

035536 

W20XXX 

■*P- 

036421 

W22 


034663 

W22X 

■y- 

035441 

W22XX 

— 

035633 

W22XXX 


036516 

X 

- 

024502 

XB 

-* 

C414C0 

XTO 

- 

040631 

XTE 

— 

040630 

Z 

- 

024451 

ZO 

— 

024476 

ZMAG2 


041425 

2NORM 

— 

041426 

ZR 

- 

041315 

zs 

— 

041322 

Z1 

- 

024065 

21SQ 

Z22 

■5“ 

024475 

033521 

Z1T 


040147 

Z1V 


040632 

Z20 

— 

033135 


START OF CONSTANTS 
007326 


START OF TEMPORARIES 
007671 

START OF INDIRECTS 
010053 

UNUSED COMPILER SPACE 
003400 
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000006 

SUBROUTINE SUM ( I ,L T X, Y t ANS ) 
DIMENSION X (L » » Y(D 

000006 

ANS=C. 

000007 

00 1 J=I,L 

000010 

1 ANS=ANS+X( J )*Y ( J-) 

000015 

RETURN 

-000016. 

END 
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SUBPROGRAM LENGTH 
000026 

FUNCTION ASSIGNMENTS 

STATEMENT ASSIGNMENTS 
I - 000010 

BLOCK NAMES AND LENGTHS 

VARIA8LE ASSIGNMENTS 
J - 000025 

START OF CONSTANTS 
000020 

START OF TEMPORARIES 
000021 

START OF INDIRECTS 
000023 

UNUSED COMPILER SPACE 
033100 



137 


000006 

000006 

OOGOIO 

000011 

0C0012 

000023 

000030 

000030 


SUBROUTINE SIMP ( DEL, NSTAT ,FGFX, ARE A> 
DIMENSION FQF-X (NSTAT I 
MS= (NSTAT-1 )/2 
T0T=0. 

DO 480 K^l t MS 

480 T0T=T0T+4.*F0FX(2*K )+2. *FOFX( 2*K+1I 

AR EA= DEL/3 .#(FCFX(1) +TOT-FOFX { 2*MS+1.) ) 

RETURN 

END 
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SUBPROGRAM. LENGTH' 

000051 

FUNCTION ASSIGNMENTS 

•STATEMENT ASSIGNMENTS 
480 - 000013. 

BLOCK NAMES AND LENGTHS 

VARIABLE ASSIGNMENTS 

K - 000050 MS - 000046 TOT 

START OF CONSTANTS 
0C0032 

START OF TEMPORARIES 
000036 

START OF INDIRECTS 
C00C42 

UNUSED COMPILER SPACE 
033100 


- OC 0047 
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SINPUT 

R 

A 

DELBAR 

BC 

'NXB 

N 

EX 

EY ' 

NUX 

NUY 

GXY 

H 

LAYER 

ST 

NSTAT 
• NORM 
IPRINT 
$END 


= 0.1E+02, 

= 0.7E+01, 

= 0.5E-*01, 

= X, 

= 0.4374E+04, 

= 7, 

= 0.105E+08, 0. 203E+07 , I, It It -It It It It It It It It It It It It It It It 
= 0.105E+C8, 0.3025E+08, I, I, I, It It It It It It It It It It It It It It It 
= 0.3E+00, 0 . 23 E -0 1-t It It It It It It It It It It It It It It It It It It 
= 0«3E+00t 0 .'346 E+0 0 1 It It It It It It It It It It It It It It It It It It 
= 0.4038E+07, 0.5249E+06, I, I, I, I, I-, I, I, It It It It It It It It It It It 
= 0.5E-01, 0.5E-01, I, I, I, I, I, I, I, I, I, I, It It It It It It It It 
= 2 , 

= 0 , 

= ■ 61 , 

= It 
= It 


IMPERFECTION SENSITIVITY OF..AXIALLY COMPRESSED CYLINDER 
CARD-SYKES A2338 RDF367 OCTOBER , 1968 ) 


LAYER EX ‘ EY NUX NUY GXY H 

1 1.050E+07 1.O5OE+07 3.000E-01 3.000E-01 4.038E+06 5.000E-02 

2 2.030E+06 3.025E+07 2.300E-02 3.46CE-01 5.249E+05 5.000E-02 


11 


12 


22 


66 



□ (J □ ^ 


1 1. 1538461 5E+07 

2 2.04628433E+06 
6. 79 2 372 93 E+0 5 
5.66031078E+02 

-1.18652215E+04 


3 .46 15 3 84 6 E+0 6 
7.08014378E+05 
2. 0 8 47764 2E +05 
1 .7373136 8E+02 
-3 .44 1905 10 E+0 3 


1 . 15384615E+07 
3 . C492 6606E+07 
2. 10155611E+06 
1.75129676E+03 
2.36927488E+04 


4.03800000E+06 
5. 24900 000 E+0 5 
2.28145000E+05 
1 . 901 20833E+02 
— 4.39137500E+03 


BETA= 8.64861327E— 01 


GAHA= 2 . 60541762E+00 Fl= -3 . 737 26721E-0 3 F2= -5.063 16344E-04 



STATION 

WO 

wax 

1 

2. 77555756E-1 7 

-9.3920E292E-02 

2 

-9. 43477662E-03 

-6 . 733074 3 3 E-02 

3 

-1. 5659 100 6E- 02 

— 3 . 940 24296E-02 

4 

-1. 86874509E-02 

-1.29633109 E-02 

5 

—1. 88273106E— 02 

9 ■ 812 24981 E— 03 

6 

-1. 65973320E-02 

2. 74849851 E-02 

7 

-1.26408642E-02 

3 -93 502793 E-02 

8 

-7.64363528E-03 

4. 537C2557E-C2 

9 

-2. 2 62 40 5 3 9E- 03 

4.6C-568413E-02 

10 

2. 930 8 2 941 E-0 3 

4. 232 53 691 E-0 2 

11 

7.48675307E-0 3 

3 • 5337404SE-02 

12 

1.10984861 E-02 

2. 63490635E-02 

13 

1 . 360 4607 IE- 02 

1.65776227E-02 

14 

1.49781 29 8 E-0 2 

7.09526524E-03 

15 

1.53049 163 E-0 2 

-1.24526635E-03 

16 

1 . 47553472E-02 

-7 . 85204416E-03 

17 

1. 35530000E-0 2 

—1 . 24044245E-02 

18 

1 • 19437024E— 0 2 

— 1.48370413E— 02, 

19 

1 . 01677135E-02 

-1.53032735E-02 

20 

8. 43 70 40 93 E-0 3 

— 1 . 41248 55 SE-02 

21 

6. 91 91 1945 E-0 3 

-1.1734 29 89 E-02 

22 

5. 72732 590E-0 3 

-8. 61621869 E-03 

23 

4.91814755E-03 

-5. 25270695 E-03 

24 

4. 49429144E-0 3 

-2 .076 663 8 5 E-03 

25 

4.412638 76 E-03 

5 . 642 57145E-04 

26 

4. 5957121 9E-0 3 

2 . 43167610E-03 

27 

4.94522392E-03 

3 . 4C815449E-03 

28 

5. 35628760E-0 3 

3.49706533E-03 

29 

5. 73098460E-03 

2. 81190581 E-03 

30 

5. 99014918E-03 

1. 55676931E-03 

31 

6.08245709E-03 

0. 

32 

5.99014918 E-03 

—1. 55676931 E-03 

33 

5. 73098460E-03 

-2.8119C581E-03 

34 

5.35628760E-03 

-3.497C6533E-03 

35 

4. 945 22 392 E-0 3 

-3. 40815449E— 03 

36 

4. 59571219E-03 

-2. 43 1676 10 E-03 

37 

4. 41263876 E-03 

-5. 64257145E-04 

38 

4. 49429144E-0 3 

2.07666385E— 03 

39 

4.91814755E-03 

5.25270695E-03 

40 

5. 72732590E-03 

8.61621869E-03 

41 

6. 91911 945 E-03 

1 . 17342989E-02 

42 

8. 437040 93 E-0 3 

1.41248559E-02 

43 

1.0 1677135E-02 

1. 5303 273 5 E— 02 

44 

1. 19437024E-02 

1 . 4837041 3 E-0 2 

45 

1.35530000E-02 

1 . 2 4044245 E-0 2 

46 

1. 47553472E-02 

7. 8 520 44 16 E-0 3 

47 

1. 53049163E-02 

1.24526635E-03 

48 

1.49781298E-02 

-7.09526524E-03 

49 

1. 360 460 71 E-0 2 

-1 . 65776 22 7 E-02 


WHY y 

-2.12972811E-01 
2. 3809496 3E-01 
2.36612173E-01 
2.13512103E-01 
1. 7487 036 3E—0 1 
1 • 2 70 9430 6 E-0 1 
7.62737380E-02 
2. 76784827 E-02 
-1.4578C364E-02 
-4. 77087467E-02 
—7. 6 271 23 93 E-02 
-8. 206036 8 9 E-02 
-8. 39C77159E-02 
-7.74229136E-02 
-6.47106142E-02 
t 4. 80927827E-02 
-2.9860C010E-OZ 
— 1. 20 6844 96 E— 02 
3. 60792435E-03 
1.5967i612E-02 
2.43136233E-02 
2* 84437569E-Q2 
2. 85928363E-02 
2.53535307E-02 
1.95772032E-02 
1.22679684E-02 
4 .*4780421 2E-03 
— 2.78888593E-03 
-6. 66133500 E-03 
-1 . 2467471 9 E-02 
-1.37842948E-02 
-1. 2467471 9E-0 2 
-8. 66133500 E-03 
-2.78888593E-03 
4. 4780421 2E-0 3 
1 . 226796 84E-Q2 
1. 9577203 2E-0 2 
2.535353 07E— 02 
2.85928363E-02 
2.84437569E-02 
2.43136233E-02 
1. 59671612E-02 
3. 6079243 5 E-03 
-1.20 6844 96 E-02 
-2. 98600010E-02 
-4.80927827E-02 
-6.47106142E-02 
-7 . 74229136 E-02 
-8* 39077159E-r02 


WQXL 

-6. 38340345 E-05 
-5.45861168E-05 
—3. 9669474 8 E-05 
-2. 15 3872 67 E-05 
— 2 . 70 9 9 94 09 E— 06 
1.45444992E-05 
2. 842 73 546 E-05 
3 .77 5 8 14 80 E-05 
4. 20 2 4636 8 E-05 
4. 13511954E-05 
3.64003894E-05 
2. 82 2 83 63 5 E-05 
1.81 1650 75 E-05 
7. 40 0997 34 E-06 
-2.68118246E-06 
-1.11155839 E-05 
-1.71 9 40 708 E-05 
— 2.05554152E— 05 
-2.11839876E-05 
-1.93 718532 E-05 
-1^56525128E-05 
-1.07 1625 48 E-05 
-5.31763699E-06 
-1.851 31 595E-07 
4.05834697E-06 
6. 95 7220 96 E-0 6 
8 . 26277345E-06 
7 .94 8 2 28 9 2 E-0 6 
6.19950036E-06 
3.38379209E-06 
0 .‘ 

-3.38379209E-06 
-6. 1995003 6E-06 
-7.94822892E-06 
-8.26277345E-06 
— 6 .95722096E-06 
-4.05834697E-06 
1.85131595E-07 
5.31 763699E-06 
1.07162548E-05 
1 .565 251 2 8 E-0 5 
1.93718532E-05 
2.11 839876E— 05 
2.05554152E-05 
1 .71940708E-05 
1 . Ill 55839E-05 
2. 681182 46E-06 
—7. 4009 973 4 E-0 6 
-1.81 16 50 7 5 E-0 5 


WOL 

—2. 16840434E-19 
-6.97393851E-06 
-1.25157747 E-0 5 
-1.61050842E-05 
— 1.75 146617E— 0 5 
-1.67992660E-05 
-1.42 5 28 78 5 E-0 5 
-1 . C3441204E-05 
-5.64042483E-06 
-7.31185017E-Q7 
3. 84133164 E-06 
7 . 63667727E-06 

1. G 35 24436 E-0 5 
1.184052Q3E-05 
1.21042672E-05 
1.1279 39 52 E-05 
9. 602 7061 0E- 06 
7.37356574E— 06 
4.91 312030E-06 
2. 52 58 6 750 E-06 
4. 67394790E-07 

-1.07897797E-06 
-2. 01 51 1955 E-06 
-2.33024837E-06 
—2. 093110 60E-06 
-1.43586290E-C6 
-5. 32013051E-07 
4.28834628E-07 
1.26656202E-06 
1. 83371 833E-06 

2. G3393063E-06 
1. 83371 833E-06 
1 . 26656202E-06 
4.28834628E-07 

-5.32013051E-07 
— 1 • 43586290E— 06 
-2. 09311060 E-06 
— 2 . 33024837E— 06 
-2.01 51 1955 E-06 
-1.0 7 897797 E-0 6 
4. 67394790E-07 
2. 52586750E-06 
4.91 312030E-C6 
7, 37356574E-06 
9.60270610E-06 
1. 12793952E-C5 
1. 2104267 2E-05 
1.18405203E-05 
1 .03524436E— 05 


H 




50 1.10984861E-02 -2. 6349C-635E-02 

51 7. 48675 307 E~ 03 -3.53374049E-02 

52 2. 930 82 941 E- 03 -4. 23253691E-C2 

53 -2. 2 624053 9E- 03 -4. 60568413E-02 

54 -7. 64363528E-03 -4. 53702557E-02 

55 — 1. 26408642E-0 2 -3 . 93502793E-02 

56 -1.65973320E-02 -2. 7484985 IE-02 

57 -1.88273106E-02 -9.81224982E-03 

58 -1. 86874509E-02 1 .29633109E-02 

59 — 1. 56591006 E-02 3 . 94024296E-02 

60 -9.43477662E-03* ' 6 . 73307433E-02 

61 2.77555756E-17 ' 9.39208292E-02 


M0DRES= 4 6. 5 790 248 IE ->-03 


-8.20603689E-02 -2. 82283635E-05 - 7. 63667727E406 

—7. 0271 23 93 E-02 -3 . 64003894E-05 3. 84133164E-06 

— 4. 77087467E— 02 -4. 13 51 1 954E-05 -7. 31 18 501 7E-07 

— 1.4578C364E-02 -4.20246368E-05 -5. 64042483E-06 

2. 7 678 482 7 E-02 -3 .77581480E-05 -1 .03441204E-C5 

7 . 627373 80 E-02 -2.84273546E-05 -1.42528785E-05 

1 . 2 70 943 06 E— 01 -1.45444992E-05 -1 .67992660E-05 

1 . 74870363E-01 2.70999409E-06 -1 .75146617E-05 

2* 13512103 E— 01 2 . 1 5387267 E-0 5 -1 . 61050 842E-05 

2. 3661 217 3Et01 3 .96694748E-05 - 1. 25 157747E-05 

2. 3 80 9496 3 E -01 5.45861168E-05 -6»97393851E-06 

2. 12972 81 IE-01 6.38340345E405 -2. 16840434E-19 



143 


STATION 

U1 

VI 

W1 

I 

4.27874577E+02 

0. 

0. 

2 

3. 68884410E+02 

6. 3299 S 795 E+0 1 

— 1 • 79131629E— 01 

3 

3.24059859E+02 

1.23208761E+02 

-3.45844320E-01 

4 

2.92892596E+02 

1.76618281E+02 

-4. 89135773E-01 

5 

2. 71058914E+0 2 

2.2130522CE+02 

-6.00786467E— 01 

6 

2. 530366G3E+02 

2 • 56227664E+02 

—6. 7629610 8E— 01 

7 

2. 3454344GE+02 

2.81514898E+02 

-7.15200574E-01 

8 

2. 13883C40E+02 

2.98232210E+02 

-7.20765051E-01 

9 

1.91937999E+02 

3.08040818E+02 

-6- 991920 72E-01 

10 

1.71158918E+02 

3.12863031E+02 

-6.58553378E-01 

11 

1. 54198476E+02 

3 . 14621605E+02 

-6.07657672E-01 

12 

1 . 42793904E+02 

3.15074002E+02 

-5. 55025237E— 0 1 

13 

1. 37235502E+02 

3. 157 26283 E+0 2 

-5.08079383E-01 

14 

1 * 36449068 E+0 2 

3.17795660E+02 

-4.72603199E-01 

15 

1. 38498385E+02 

3.22193415E+02 

-4.52460365E-01 

16 

1.41231000E+02 

3. 2951271 8E+02 

— 4. 4954695 IE-01 

17 

1 . 42827650E+0 2 

3.4001S966E+02 

—4. 63926709E— 01 

18 

1.421 18 431 E+0 2 

3. 5365747 5 E+02 

—4. 9410091 2 E~0 1 

19 

1.38641740E+02 

3.70067262E+02 

-5.37368404E-01 

20 

1.32506789E+02 

3 • 88641394E+02 

— 5. 902 3 653 8E— 01 

21 

1.24160756E+02 

4.08596884E+02 

— 6. 48 8461 8 8 E—0 1 

22 

1. 14158859E+02 

4. 29065796E+02 

— 7.09374473E-01 

23 

1.030G2921E+02 

4. 4918673 8E+02 

-7. 6837967 8E-01 

24 

9. 106 9 862 1E+ 01 

4. 681 8 345 6 E+0 2 

— 8. 230 5721 8E— 01 

25 

7. 86136552E+01 

4. 85419537E+02 

-8.71384739E-01 

26 

6.5 8042 848 E+0 1 

5.00423956E+C2 

-9.12148347E-01 

27 

5.27678 145 E+0 1 

5. 12888345E+02 

— 9.44858207E-01 

28 

3. 96068892E+01 

5.22642050E+02 

— 9. 69 5771 81 E-0 1 

29 

2. 640G3884E+0 1 

5. 2 96 13 55 6 E+02 

-9. 86698008E-01 

30 

1. 3194551QE+Q 1 

5 . 33790900E+02 

— 9. 967 10530 E—0 1 

31 

5.79721572E-04 

5.35181365E+02 

-l.OOOOOOCOE+OG 

32 

- 1. 31 93405 6E+01 

5. 3379 1041 E+02 

— 9. 96710842E-01 

33 

-2. 63992 83 1E+ 01 

5. 29614238 E+02 

-9.86698591E-01 

34 

-3. 96058450E+01 

5.22642474E+02 

-9.69577954E-01 

35 

—5. 27668447E+01 

5. 1288 891 5 E+0 2 

—9 . 448 5 906 5 E-0 1 

36 

-6.58033942E+01 

5. 0042468 1E+02 

-9. 12149174E-01 

37 

-7.86128400E+01 

4. 8542C436E+02 

- 8. 7 138 542 5 E-0 1 

38 

-9. 10691097E+01 

4.681 84558E+02 

- 8. 2 30 5768 5 E-0 1 

39 

-1.03002 210 E+0 2 

4.49188C92E+02 

-7. 6 83 7990 IE— 01 

40 

—1.1 41 58 159 E+0 2 

4.29067465E+02 

-7. 0937450 2E-01 

41 

-1 . 24160022E+02 

4. 0859 895 2 E+02 

—6. 4884616 9E—0 1 

42 

-1.32505966E+02 

3. 88643962E+02 

-5.9023672 8E-01 

43 

-1. 38640762E+0 2 

3. 70070443E+02 

-5.37369177E-01 

44 

-1.42117232E+02 

3. 53661390E+02 

— 4. 941C2757E-01 

45 

-1. 428261 80E+0 2 

3 . 40024734E+02 

-4. 63930 21 5 E-0 1 

46 

-1. 41229253E+02 

3 . 2951 8447E+G2 

— 4.49552778E— 01 

47 

— 1. 3849642 2E+0 2 

3 . 2220 0193 E+02 

—4. 524 69 197 E—0 1 

48 

-1. 36447C44E+02 

3 . 1780 3 54 2 E+02 

-4.72615684E-01 

49 

-1.37233671 E+ 02 

3 . 1573 5 28 2 E+02 

-5.0 809605 4E-01 



aM 


50 

—1. 42 792 6 01 E+0 2 

. s , - 

3 ." 150 840 7 6 E+0 2 

— 5.’55046419E— 0 1 


51 

-1 .-541-9 8064E+0'2 

3. 1463 2649 E+02 

-6. 07683396 E-01 


52 

—1 . 71159717E+02 

3 .12 874866 E+02 

-6. 58583293 E-01 


53 

— 1.9194021 0E+0 2 

3.08053175E+02 

-6. ,992 2 53 8 8 E-01 


54 

-2.13886'699E+G2 

2 . 98244729E+02 

-7.20800513E-01 


55 

-2.34548436E402 

2.81527127E+02 

-7. 1 523 6500 E— 01 


56 

-2. 53042773 E+02 

* 2. 562 3 9074 E+02 

— 6. 763 3048 4E— 01 


57 

— 2.71066189E+02 

2. 2131S241E+02 

-6. ‘008 171 02 E-01 


58 

-2. 929.01156 E+02 

1.766 26.3 5 2 E+0 2 

-4.891 605 0 6 E— 01 


59 

-3. 24070207E+02 

1 • 232 1440 4E+02 

— 3.45861254E— 01 


60 

—3. 68897304E+02 

6. 3 ?0 2 8 65 3 E+0 1 

- 1'. 7-91 3 9344E— 0 1 


61 

—4. 27890803E+02 

0. 

• 2.26829198E-06 


STATION 

U1 

VI 

4 

W1 

Ml 

1 

1 . 13129929E— 0 2 

0. 

0. 

0. 

2 

9. 75329445E-03 

1. 673649 8 5E-03 

-1. 791 3162 9E— 02 

— 2. 9869731 8E+01 

3 

8.56813445E-03 

3, 2 5763 651 E-03 

-3.45844320E— 02 

-5.62318563E+01 

4 

7.74407280E-03 

4. 66978285E-03 

— 4. 89 135773 E— 02 

-7.57065770E+01 

5 

7.16679C81E-03 

5. 85130441E-03 

-6.C0786467E— 02 

— 8.59047813E+01 

6 

6. 69028T37E-03 

6 . 77465294E--03 

-6.76296108E-02 

— 8.59423329E+01 

7 

6. 20132^60 E-03 

7. 4432467 5E-03 

—7 . 15200574E— 02 

—7 .65 12545 1E+01 

8 

5. 65506215E-03 

7.88525205E-03 

-7. 2076505 IE-02 

-5.95900447E+01 

9 

5.07483583E-03 

8 . 144 5 913 9 E—G 3 

-6.99192072E-02 

-3 .79 2 39207 E+0 1 

10 

4. 52543744E-03 

8.272C9Q61E-03 

-6. 5 855 33 7 8E— 02 

-1.44 8446 99E+01 

11 

4. 07 700378 E-03 

8. 318 5 £7 2 3 E-03 

-6.07657672E-02 

8.01080816E+00 

12 

3. 77546721 E-03 

8. 3305486 1E-C-3 

- 5. 5502523 7E— 02 

2.74201653E+01 

13 

3.62850318E-03 

8 . 34779491E— Q3 

-5. 080793 83E-02 

4.23300436E+01 

. 14 

3. 60770990E-03 

8.402 5C917E-03 

—4. 726031 99E-02 

5.205 943 5 9 E+0 1 

15 

3. 66189379E-03 

8. 5 187 8570 E-03 

— 4.52460365E-02 

5.65638317E+01 

16 

3. 7341 440 6E- 03 

8.7123C789E-03 

-4.49546951E-02 

5 • 629 1 993 5 E+0 1 

17 

3.77635 943E—0 3 

8. 9901 1927E— 03 

—4. 63926709E— 02 

5.20375922E+01 

18 

3 . 7576G770E-03 

9.350 69467E— 03 

— 4.941G0912E-02 

4. 480 7 972 6 E+0 1 

19 

3. 665684C9E-03 

9.78456901E-03 

— 5. 37368404E-02 

3.57137642E+01 

20 

3. 50347613E-03 

1.02756686E-02 

-5.90236538E-02 

2. 5872551 8E+01 

21 

3. 28280724E-0 3 

1 .0803291 2E-02 

— 6. 48 8461 88 E— 02 

1.63190 273E+01 

22 

3. 01 83 5734 E-03 

1 • 13444887E-02 

— 7.09374473E-02 

7.920065S3E+00 

23 

2.7233 9461 E-03 

1 . 18764860E-02 

— 7. 68379678E-02 

1.. 3009661 OE+OO 

24 

2. 40788484E-03 • 

1 « 237 87587E-02 

-8.2305721 8E-02 

— 3.20585198E+00 

25 

2.07854305E-03 

1-.28344802E— 02 

— 8. 71384739E-02 

— 5 .57 837291 E+OO 

26 

1 « 73986362E-0 3 

1. 3231 1 967 E— 02 

—9.1 21 48347E— 02 

— 6 .08 214277E+00 

27 

1 • 3 95 17967 E-03 

1 . 3560754 8 E— 02 

— 9.44858207E-02 

— 5.20791769E+00 

28 

1.04720515E-03 

1 . 381 8 642 5 £-02 

—9. 69577181E-02 

—3 . 57 828578E+00 

29 

6. 98025599 E-04 

1.40029795E— 02 

-9. 8669800 8E-02 

-1.8 3 750 192 E+OO 

30 

3.48863594E-04 

1.411 3417 8E-02 

-9.9671C530E-02 

—5. 41 60031 2 E-01 

31 

1. 53278237E-08 

1.41501817E— 02 

-l.OOQOOOOOE-Ol 

— 6.55972735E-02 

32 

—3 . 488333 09 E— 04 

1.4}134216E-C2 

—9. 967 1C 84 2 E— 02 

—5 .41 7095 8 3 E-01 

33 

—6. 97996375E— 04 

1.4C029870E— 02 

—9.8669 85 9 IE— 02 

— 1 ♦ 83 7705 81 E+OO 

34 

-1.04717754E-03 

1 . 38186537E— 02 

-9. 69577954E-02 

—3. 57855561 E+OO 

35 

-1 . 39515403E-03 

1 .35607699E— 02 • 

—9 • 448 590 6 5 E— 02 

— 5 . 20821202E+00 
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36 

-1. 73984007 E-03 

1. 3231215 EE-02 

-9. 12149174E-02 

-6.0824Q945E+00 

37 

-2.07852149E-03 

1 . 2834504CE— 02 

-8.71385425E-02 

— 5.578551 96E+00 

38 

-2.40786495E-03 

1.23787878E-02 

-8. 230 57685 £-02 

-3.20587938E+00 

39 

-2.72337582E-03 

1 . 1876521 8E-02 

-7. 68379901E-02 

1.30115385E+00 

40 

-3.01833882E-03 

1.13445329E-Q2 

-7.0937450 2E-02 

7.92052636E+00 

41 

-3. 28278784E-03 

1.08033459E-02 

-6. 4884616 9E-02 - 

1.63198065E+01 

42 

—3. 50345437E-03 

1.027 57 36 5 E— 02 

-5.90236728E-02 

2.58736773E+01 

43 

-3.66565823E-03 

9. 7846 5311 E-03 

— 5. 37369177E— 02 

3.57152387E+01 

44 

—3. 757576 OGE-'O 3 

9.35079817E-03 

-4.941C-2757E-02 

4 • 480 9766 6 E +01 

45 

-3.77632058E-03 

8.99024533E-03 

-4. 6393021 5E-02 

5.20396372E+01 

46 

— 3. 73409786E— 0 3 

8. 71 24593 8 E-03 

— 4.49552778E-02 

5.62941760E+01 

47 

—3. 661 841 69E-03 

8 . 5 1896 493 E-03 

—4. 5 24 691 97 E-0 2 

5.65659903E+01 

48 

-3. 60765639E-0 3 

8.40271758E-03 

-4. 726156 84E-02 

5.20613637E+01 

49 

-3. 62845478E-03 

8. 348C3284E— 03 

-5.08096054E-02 

4.23314989E+01 

50 

-3.77543275E-03 

8 . 3308 149 7 E-03 

—5. 55G4641 9E-02 

2.74208924E+01 

51 

-4.07699289E-03 

8.31887S26E-03 

— 6.07683396E-02 

8.01057054E+00 

52 

-4. 52545856 E-03 

8 . 27240351E— 03 

—6. 58583293E-02 

-1.44858485E+01 

53 

-5. 07489428E-03 

8. 14491 81 IE— 03 

—6.9922538 8E-0 2 

-3.79265118E+01 

54 

-5.65515888E-03 

7. 885 5 8307 E-03 

— 7.2080C513 E-0 2 

— 5.95937764E+01 

55 

-6.20145470E-03' ' 

7.44357C09E-03 

-7, 15236500 E-02 

-7.65171808E+01 

56 

—6 • 69044448E— 0 3 

6.774-9 5462E-03 

-6.76330484E-02 

— 8 .59474759E+01 

57 

-7- 16698315E-03 

5 . 85 1 569 34 E— 03 

—6.0081710 2E— 02 

-8.59099117E+01 

58 

—7 . 744299 12E-0 3 

4. 66999626E-03 

-4.89160506E-02 

-7.57111183E+01 

59 

—8. 56840 804E-0 3 

3. 2 577 8 572 E-03 

-3.45861254E-02 

— 5.62352613E+01 

60 

-9. 75363 534E-03 

1. 67372615E-03 . 

-1.791 3 93 44 E-0 2 

-2.98715665E+01 

61 

—1. 131-34220E-0 2 

0. 

* * ► 
* 

- 2. 268291 58E— 07 

0. 

STATION 

W1X 

U1.X 

. V1X 

W1XX 

1 

-1.58863798E-01 

-1.497.4CC93E-02 

1.47298410E-02 

1.11916155E-03 

2 

— 1.4821 8 994E-01 

-1.17636793E-02 

1.-39612993E-02 

9. 1241 17 61 E-02 

3 

-1.32858919E-01 

-8 . 61094 99 IE- 03 

1. 28405700E-02 

1.72074401 E-0 1 

4 

— 1.G9260920E-01 

—6 .0057 5 845 E-03 ' 

1.11157196E-02 

2.32462715E-01 

5 

-8.02115725E-02 

—4. 5 16248 9 8 E-03 

9..02087179E-03 

2. 65 5 26 106 E-0 1 

6 

-4. 90346 170E-0 2 

-4. -13772091 E-C3" 

6.82261004E-03 

2 *6893 5987 E-0 1 

7 

-1. 90581 181E-02 

-4. 43665381E-C3 

4.7597 10 50 E-0 3 

2.44946853E-01 

8 

6. 86078625E-03 

—4. 82780045E-03 

3.C0576276E-03 

1 .9937722 IE-01 

9 

2. 66621454E-02 

-4.841 24876E-03 

1.657879 51 E-0 3 

1 .400746 51 E-0 1 

10 

3. 92290287E-02 

-4.276423C9E-03 

7.45696464E-04 

7.53576338E-02 

11 

4. 43692034E-02 

-3. 21415814E— 03 

2 . 505342 8 6E— 04 

1 .27596463E— 02 

12 

4. 26764094E-02 

—1.9 22 14544 E-03 

1.25175766E-04 

-4. 17789723 E-02 

13 

3.532373C6E-02 

-7.18959870E-C4 

3.08402421E-04 

-8.42669499E-02 

14 

2.38367223E-02 

1.43102627E-04 

7..32817662E-04 

-1. 12653192 E-01 

15 

9. 88124905E-03 

5.4186066CE-04 

1. 32 770 87 8E— 03 

-1. 26583493 E-01 

16 

-4. 91414732E-03 

4. 905 6703 3 E—C 4 

2.02000102E— 03 

-1 . 27051 874E-01 

17 

-1 . 90945547E-02 

1.0055S471E-04 

2. 73594334E— 03 

-1. 16040 824E-01 

18 

—3 . 147 50125 E— 0 2 

—4. 7432288 9E—C4 

3. 40478460E— 03 

— 9 , 619 55957E— 02 

19 

-4. 12009 82 5 E- 02 

-1.0 8 91 3 52 8 E—C 3 

3. 9 641 73 96 E— 03 

-7 .053531 8 2E-02 

20 

-4. 77761930 E-0 2 

-1.6409C079E-03 

4. 36595239E-03 

-4.21 82 5761 E-0 2 
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21 

-5.10591150E-02 

-2 • 0790 805 5 E— 03 

4. 5 806 57 80 E-C 3 

—1 .409608756-02 

22 

-5.12286385E-02 

-2.39746267E-C-3 

4. 5994C 594E-03 

1 . 1189971 1 E-02 

23 

—4. 87211768E-0 2 

—2. 61631070E— G3 

4. 432 5854 2E- 03 

3. 17950871 E-02 

24 

—4. 41450261 E-0 2 

-2. 763 64 957 E— 03 

4. 10568981 E-03 

4.66532105E-02 

25 

-3. 81819123E-02 

-2.862948C9E-03 

3. 65330574E-03 

5.55715968E-02 

26 

-3. 14886292E-02 

— 2. 92870017E— 03 

3. 1 126C537E— 03 

5. 91703990 E-02 

27 

-2.46123573E-02 

— 2. 96853632E— C3 

2. 51762505E— 03 

5. 870 854 80 E-02 

28 

-1.79313434E-02 

—2. 9878G318E— 03 

1 . 895 2 4876 E-0 3 

5. 582 31 191 E-02 

29 

-1.1 62 85781 E-02 

-2.99289236E-03 

1 . 26332277E— 03 

5. 2 2 2428 50 E-02 

30 

-5.70085373E-03 

— 2. 99147259E-03 

6. 30866394E— 04 

4.9393 847 8 E-0 2 

31 

-1. 34012 557E-07 

— 2.99012959E-03 

1. 59941 27 IE— 08 

4. 833277 60E-02 

32 

5 . 700 60391 E-0 3 

— 2.99147873E-03 

-6. 3083441 5 E— 04 

4.93941597E-02 

33 

1. 162838C6E-02 

— 2.9929G383E-C3 

-1.2 63 2 9070 E-03 

5.22248689E-02 

34 

1.79312252E-02 

— 2.98781853E-C3 

-1. 89521618E— 03 

5. 5 82 3 8967 E-0 2 

35 

2 .46 123346E-02 

— 2.96855373E-03 

-2.51759C97E— 03 

5.87094064E-02 

36 

3. 14887031E-0 2 

—2 . 92871 76 8 E-C 3 

-3. 11256814E-03 

5. 91 71 19 72 E-02 

37 

3.81 820664E— 02 

—2. 8629 6376 E— 03 

-3. 65326291E— 03 

5. 5 5 721736 E-0 2 

38 

4.41452246E— 02 

— 2.76366142E— 03 

—4. 10563819E— 03 

4. 66533959E-02 

39 

4. 87213645E-02 

-2. 61631658 E-C3 

— 4. 43252119E-03 

3. 179471 57E-02 

40 

5. 122 8 741 9 E-0 2 

-2. 39748006 E-03 

— 4. 59932496E— 03 

1 . 11888976E-02 

41 

5. 1059045 9E— 02 

— 2.079C 666 8E— C-3 

-4. 5 805 56 02 E-0 3 

-1.40979714E-02 

42 

4.77758537E-02 

— 1. 64087313E-03 

— 4.36582635E-03 

-4. 21853242E-02 

43 

4.12002733E-02 

-1 .08909268E— 03 

-3.96402135E-03 

-7. 05389113 E-02 

44 

3. 14738412E-02 

-4.74267201E-O4 

— 3.40460476E-03 

-9. 61999244E-02 

45 

1.90928481E-02 

1.0062C599E-04 

— 2. 73573767E-03 

-1.16045672E-01 

46 

4. 91186 450 E-0 3 

• 4.90622972E-04 

-2.01977318E-03 

-1.27056905E-01 

47 

-9.88410270E-03 

5. 41*89200 8E-C4 

—1. 32746484E— 03 

— 1.265 88248 E-01 

48 

-2. 3 840081 2 E-0 2 

1 . 430 8762 0E-C4 

-7. 32 5660 83 E— 04 

—1 . 12657099E— 01 

49 

-3. 532 74 576 E— 02 

—7.1 9041 54 5 E— 04 

— 3.C81 5405 9 E— 04 

-8.42693537E-02 

50 

-4. 26 8C2 896 E-0 2 

-1 . 9223061 8 E-03 

-1. 24943 930E-04 

-4.17791948E-02 

51 

-4.43729464E-Q2 

-3. 21439633E— 03 

-2. 50334830E-04 

1.27622210E-02 

52 

—3. 92322 823 E-0 2 

-4.27672C26E-03 

—7 . 45547749E— 04 

7. 536344 93 E-02 

53 

-2. 66645226E-0 2 

-4. 841572t79E-03 

-1. 6 57 80 190 E-03 . 

1.40083860E-01 

54 

-6.86190525E-03 

-4. 8281 1608E-03 

-3.0C577726E-03 

1 .99389580E-01 

55 

1. 90585836E-02 

-4.43693829E-03 

—4. 759836 19E— 03 

2. 449 61 6 5 8E— 01 

56 

4. 90368851E-02 

-4.13797906E-C3 

-6. 82286034E-03 

2. 68 9 52081 E-01 

57 

8.02157C48E-02 

-4. 51651987E-03 

-9. 0212 50 11 E-03 

2.65541973E-.01 

58 

1. 09266792 E-01 

-6.0061C67CE-03 

-1.11162155E-02 

2. 32 4766 66 E-01 

59 

1.32866212E-01 

-8 • 61 14409 5 E-C 3 

— 1.28411576E-02 

1.72084825E-01 

60 

1. 48227224E-01 

-1 . 17643454E— 02 

-1.396 193 88 E-0 2 

9 .1 2467946E— 02 

61 

1.58872683 E-01 

-1.49748 540 E-0 2 

-1.47305055 E-02 

1.11927182 E-03 


STATION 

NX1 

NY1 

NXY1 

NYOL 

1 

-2. 29851527E+01 

-7.27850388E+00 

2.53052Q34E+03 

—3 . 17916866E— 01 

2 

—2. 2814709GE+0 2 

-1. 5683C261E+03 

2. 34620711E+03 

-1.74861622E+0C 

3 

-4.38673734E+02 

-2.98889820E+03 

2. 16034766E+03 

-2 . 88541 230E+00 

4 

-6.14815772E+02 

—4. 13432232E+03 

1.86971 991 E+0 3 

-3.62010975E+00 

5 

-7.4895C954E+02 

—4. 907 2 5825 E+0 3 

1. 50079939E+03 

-3.90657018E+00 
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6 

— 8. 42140079E+0 2 

—5. 2497040 2 E+0 3 

1.0 864 1051 E+0 3 

-3.75656675E+00 

7 

-9.G0767759E+02 

—5. 1491 746 8 E+0 3 

6 • 621 6490 9E+02 

-3 . 2309353 8E+00' 

8 

-9.32482617E+02 

-4. 63949754E+03 

2. 62692 81 8E+02 

— 2 .426261 75E+00 

9 

-9.43416484E+02 

-3. 7953792 8E+03 

-8.18588385E+01 

-1,459281936+00 

10 

—9. 37464050 E+0 2 

-2.72144103E+C-3 

—3. 4862514 1E+02 

-4.51095177E-01 

11 

-9.17188248E+02 

•*1. 53763145E+G3 

-5.23796256E+02 

4.87009673E-01 

12 

-8.853C2018E+02 

-3. 63633603E+02 

— 6.032C7130E+02 

1.26475826E+00 

13 

—8. 45752 85 6 E+0 2 

6.95081348E+02 

-5.91772022E+02 

1.82029297E+00 

14 

»- 8. 03 92 7 663 E+0 2 

1 . 557 50505E+03 

-5.01962552E+02 

2.12347334E+00 

15 

-7. 66033376E+02 

2. 1726165 8E+03 

—3. 51 647 10 9 E+0 2 

2. 17525954E+00 

16 

-7.38040382E+02 

2. 5202182 5E+03 

-1. 61655900E+02 

2.00375664E+00 

17 

-7. 24638397E+02 

2. 60793678E+03 

4. 65892958E+01 

1 . 65777415E+00 

18 

-7. 2852932 1E+0 2 

2.46565133E+C3 

2. 5311791 3E+02 

1 .198903 13 E+00 

19 

— 7.50191921E+02 

2.13863200E+03 

4.4096131 1E+02 

6.93138923E-01 

20 

-7.88088504E+02 

1 . 6 804496 8 E+0 3 

5. 971 410 2 8 E+0 2 

2 .02 9 953 99 E^O 1 

21 

-8. 3 91 79903 E+0 2 

1.14641334E+C3 

7. 1309973 4 E+0 2 

-2.19108625E-01 

22 

’-8. 99572954E+02 

5. 88Q23364E+02 

7 . 8463353 8E+0 2 

-5.35631119E-01 

23 

-9. 65133732E+02 

4. 87413382E+01 

8.11419440E+02 

-7.26544464E— 01 

24 

—1.031 9473 5 E+0 3 

— 4.38745186E+02 

7, 9 62 536 99 E+0 2 

-7.897646G8E-01 

25 

-1.09657446E+03 

— 8 . 53682478E+02 

7 . 4412 881 2 E+0 2 

-7.395328166-01 

26 

-1. 15612420E+03 

-1 . 18686432E+C3 

6.61278834E+02 

-6 .030 9745 3E— 01 

27 

-1.20821159E+03 

-1.43886435E+03 

5.54313647E+02 

—4 . 1 6179227E-0 1 

28 

-1. 250881 26E+03 

— 1 . 6172C243E+C3 

4. 2954143 5 E+0 2 

-2.17783623E-01 

29 

-1.28256052E+03 

-1.732 89 805E+03 

2.92545875E+02 

— 4. 49 4 84240 E— 02 

30 

-1.302C6656E+03 

— 1 .79695563E+03 

1.48043452E+02 

7.20141581E-02 

31 

—1 . 30865394E+03 

-1. 81734356E+G3 

2 . 02 5002 8 4 E— 03 

1.13295408E401 

32 

-1.302C6682E+03 

-1. 79695720 E+03 

-1.4803 946 5 E+0 2 

7.20141581E-02- 

33 

—1. 28256104E+03 

-1 . 73290027E+03 

-2. 92 54203 8 E+G 2 

-4. 49484240 E-02 

34 

-1. 25088202E+03 

-1 • 61 72035 3 E+03 

~4.'29537727E+02 

-2.17783623E-01 

35 

■’-1. 2O821262E+03 

—1.438 86 189 E+03 

-5. 54309875 E+0 2 

-4.16179227E-01 

36 

—1. 15612553E+03 

-1.1868554CE+03 

-6. 612745 87E+02 

-6.03097453E-01 

37 

—1 . 09657 61 6 E+0 3 

— 8 . 53664050E+02 

-7.44123437E+02 

-7.39532816E-01 

38 

—1 . 03 19495 6 E+0 3 

—4. 38714395E+02 

-7. 962463 04E+02 

-7. 89764608E-01 

39 

—9. 6 51 3 6 661 E+0 2 

4. 87867629E+01 

-8. 1140 8923 E+0 2 

-7 . 26544464E-01 

40 

-8.99576929E+02 

5 • 88084697E+02 

-7. 84618653E+G2 

—5. 35631119 E-01 

41 

-8. 39185370E+02 

1 . 1 464904 5E+0 3 

— 7.13079190E+02 

-2.19108625E-01 

42 

-7. 88096C34E+0 2 

1 .680 5406 5 E+0 3 

-5.97113618E+02 

2. 02 99 539 9 E-01 

43 

—7. 50202195E+02 

2. 13 8732 84 E+03 

-4.40926074E+02 

6. 93138923 E-01 

44 

-7.28543091E+02 

2. 4657 557 6 E+03 

— 2. 5 30743 G9 E+0 2 

1 .19890313 E+00 

45 

-7.24656^146+02 

2. 608C3634E+03 

— 4. 6537390 1E+01 

1 .6577741 5 E+00 

46 

-7.38063310E+02 

2. 5203C250E+03 

1.61 71 5267 E+0 2 

2. 00 375664 E+00 

47 

-7. 66061684E+02 

2. 172 673 79 E+03 

3.51712193E+02 

2.17525954E+00 

48 

-8. 03961 530 E+0 2 

1. 5 5-7 5 23 L5 E+0 3 

5.0203063 8E+02 

2.12347334E+00 

49 

—8 • 45792 107E+0 2 

6. 9 504925 5 E+0 2 

5- 91 83945 5 E+0 2 

1 . 82029297E+00 

50 

-8.85346112E+02 

— 3. 63724462E+C2 

6. 0326 946 2 E+0 2 

1.26475826E+00 

51 

-9. 1723 632 5 E+0 2 

-1.53778559E+03 

5. 2384851 6E+02 

4.87009673E-01 

52 

-9. 37515018E+02 

—2. 72 1657 37 E+0 3 

3. 4 86622 2 5 E+0 2 

— 4. 51095177E-01 

53 

-9.43469102E+02 

— 3.79565C5CE+C3 

8.18759885E+01 

-1.4’5928193E+00 

54 

— 9.32535524E+0 2 

— 4.639 8C 570 E+03 

-2. 626 99493 E+0 2 

-2. 42 6 2 617 5 E+00 

55 

-9.00819418E+02 

— 5.149 50791E+03 

-6.621 97928E+02 

—3. 230 93 53 8 E+00 

56 

-8.42188665E+02 

—5. 25003397E+03 

-1.08647060E+03 

—3 .75656675 E+00 

57 

-7.48994272E+02 

— 4.90755819E+03 

—1. 50088519E+03 

-3.9065701 8E+00 

58 

-6. 14851329E+02 

—4. 13456559E+03 

—1. 86982789E+03 

-3 . 620 10 975 E+00 

59 

-4.38699C44E+02 

-2.989C6C8CE+03 

-2.1 604722 1E+03 

— 2.88541230E+00 

60 

-2.28160185E+02 

-1 . 5 683654 8 E+03 

— 2.34634087E+03 

-1.74861622E+0Q 

61 

-2.29864352E+01 

-7. 22995077E+OG 

—2. 5 30 6 6 23 7 E+0 3 

—3. 17 91 686 6 E-01 
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STATION 

U2C 

V2C 

W20 

M20 

1 

2. 41273601 E-02 

0 . 

0 . 

1.01 89651 3E+00 

2 

2. 13269 831E-02 

0. 

— 1.97657 11 4E— 02 

-3 . 20432184E+01 

3 

1. 93071340E-02 

0. 

-3.56372344E-02 

-6. 35451996E+01 

4 

1.79693584E-02 

0. 

-4. 6886 80 29 E-02 

-8.92656410E+01 

5 

1.70953699E-02 

0. 

—5. 31752511 E-02 

— 1.04189883E+02 

6 

1. 64421 698 E—0 2 

0. 

-5. 45 5 8 69 5 4 E— 02 

-1.049C3861E+02 

? 

1.5823 5981 E-02 

0. 

-5. 14701927E-02 

-9.11309865E+01 

8 

1.51554210E-02 

0. 

-4. 4666845 8E-02 

-6 . 58459794E+01 

9 

1.445 5 2 341 E-02 

0. 

-3.51430195E-02 

-3. 41 5478 30 E+01 

10 

1.3 8060730 E-02 

0. 

-2. 40 19 8473 E— 02 

-1.64196988E+00 

11 

1. 33046449E-02 

0. 

-1. 2427121 6E-02 

2.70780646E+01 

12 

1.30165558E-02 

0. 

— 1 . 3948920 5E— 03 

4.91651691E+01 

13 

1 . 2953 0808 E-02 

0. 

8.23092779E-03 

6 .36445928 E+01 

14 

1. 30723843E-02 

0. 

1. 58405432E-02 

7.09840125E+01 

15 

1.32984847E-02 

0. 

2. 107830 8 9E— 02 

7.24552172E+01 

16 

1. 35469518E-0 2 

0. 

2. 38311879E— 02 

6 .954263 26 E+01 

17 

1 . 37473 104E—0 2 

0. 

2. 41 9 5674 2 E—0 2 

6 . 3 55425 72 E+01 

18 

1 . 385 6241 2E- 02 

0, 

2. 24322239E-02 

5. 5474453 6E+01 

19 

1. 38603356E-02 

0. 

1.891 5131 0E-02 

4. 6009 1938 E+01 

20 

1. 37706268E-G2 

0. 

1. 40835843E— 02 

3. 57 274431 E+01 

21 

1 , 36127682E-02 

0. 

8. 3974145 5E-03 

2 .51 9 7 865 8 E+01 

22 

1.34167499E-Q2 

O'. 

2. 29945202E-03 

1. 504975 56E+01 

23 

1. 3208991 8E-0 2 

0. 

—3. 81436439E— 03 

5.93590994E+00 

24 

1. 300 81 21 5 E-02 

0. 

—9. 6162157GE-03 

-1.57654950E+00 

25 

1.28242773 E-02 

0. 

-1.48581851E-02 

-7 . 131791 43 E+QO 

26 

1. 26607744E-02 

0. 

-1. 93742183E-02 

-1 .06698764E+01 

27 

1.25166408E-02 

0. . 

—2. 30727241 E-02 

— 1.243 5028 1E+01 

28 

1. 23888272E-02 

0. 

—2. 59214261E-02 

-1.29046665E+01 

29 

1. 22735347E-02 

0. 

-2. 7927246 8E-02 

— 1 . 26 6 3 53 90 E+01 

30 

1.2 1667257 E-02 

0. 

—2 , 91148460E— 02 

—1.2 2600221 E+01 

31 

1. 20642136E-02 

0. 

—2. 9 5076971 E-02 

-1.20796955E+01 

32 

1. 1961 69 95 E-02 

0. 

-2.91151449E-02 

—1.22608780 E+01 

33 

1. 18548845E-02 

0. 

-2.79278128E-02 

-1.26651 85 8E+01 

34 

1 . 17395827E-02 

0. 

— 2.59221977E-02 

—1 .290 6976 8E+0 1 

35 

1. 16117573E-02 

0. 

— 2. 30736150E— 02 

-1.24378182E+01 

36 

1.14676102E-02 

0. 

— 1. 93751238E— 02 

— 1 .06729160E+0 1 

37 

1 • 1304P930E-02 

0. 

— 1.48589903E— 02 

-7.13481736E+00 

38 

1. 11202342E-02 

0. 

— 9.61680552E— 03 

— 1.57928304E+00 

39 

1.09193494E-02 

0. 

-3.8 1463400E-03 

5.93374460E+00 

40 

1.G7115781E-02 

0. 

2. 29958606 E-03 

1.50 4841 33 E+01 

41 

1. 051554 90E-02 

0. 

8. 3980038 2E— 03 

2.51975647E+01 

42 

1. 0357683 8E-02 

0. 

1.4084639 2E— 02 

3. 57 2 835 18E +01 

43 

1.0 2679 749 E-02 

0. 

1, 89166 137E-02 

4.601141 89E+01 

44 

1.02720775 E-02 

0 . 

2. 24340442E-02 

5.54780252E+01 

45 

1.03810258E-02 

0 . 

2. 41 97 690 IE— 02 

6. 3 5 591 085 E+01 

46 

1.05814C95E-02 

0 . 

2. 383321 11 E-02 

6.95485724E+01 

47 

1.08299C 56E-02 

0 . 

2. 108011 5 5E-0 2 

7. 246 18991 E+01 

48 

1.10560321E-02 

0 . 

1.58418914E-02 

7.09909102E+01 

49 

1. 11753506E-02 

0 . 

8.23158018E-03 

6.36509993E+01 
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50 

1. 1111 8726E-0 2 

0 . 

•*1 • 39514209E— 03 

4.91702391E+01 

51 

1.08237583E-02 

0 . 

-*1« 24284224E-02 

2 .70 8091 63 E+0 1 

52 

1.03222844E-02 

0 . 

-2.40222635E-02 

-1.64208996E+00 

53 

9. 6730635 5E— 03 

0 . 

-3.51465 125E-0 2 

-3.41582972E+01 

54 

8. 97281262E-03 

0 . 

— 4.467 12635 E— 02 

-6. 58527922 E+01 

55 

8. 3045753 9E-0 3 

0 . 

-5. 14752709 E-G2 

-9.11403898E+01 

56 

7. 685 94 871 E- 03 

0 . 

-5.45640712E-02 

-1.04914601E+02 

57 

7. 03268 93 2 E- 03 

0 . 

-5.31804874E^02 

-1 .04200410 E+0 2 

58 

6. 15861785E^03 

0 . 

-4.68914202E-02 

-8. 92 7 448 3 8 E+01 

59 

4. 82071072E-03 

0 . 

-3.56407468E-02 

-6.35513219E+01 

60 

2. 80065950 E—0 3 

0 . 

— 1. 9 767663 IE— 02 

—3. 204 619 07 E+01 

61 

0 . 

0 . 

0 . 

1.01907911E+00 
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STATION 

U22 

V22 

W22 

M22 

1 

—8. 90431061E-02 

0. 

0. 

1.01896513E+00 

2 

-7.30274709E-02 

—1 .4705648 8E-02 

1 .90 802 93 6 E-0 1 

4.00389061E+02 

3 

-5. 98268948E-02 

-2. 81 9 48 330 E-0 2 

3.68258957E-C1 

7 . 56210448E+02 

4 

—4. 942 5 8779 E-0 2 

-3. 938 2 170 5 E-02 

5.17653284E-01 

1 .01719279E+03 

5 

-4.07019003E-02 

-4. 752 23 11 2E— 02 

6.27833896E-01 

1.14673573E+03 

6 

—3. 2251 843 5 E-0 2 

-5. 2 327 87 8 5 E— 02 

6.92788854E-01 

1.13440683E+03 

7 

-2. 3 146 2 475 E-0 2 

-5 .3943 1811 E-02 

7 . 121 6613 9E— 01 

9.97377371E+02 

a 

-1.32884908E-02 

—5 . 2 8227 55 8E— 02 

6. 9066 16 04 E— 01 

7 .733828818+02 

9 

-3.31278379E-03 

-4. 9577S861E-02 

6. 3 65 6 1711 E-0 1 

5.09255811E+02 

10 

5. 81176682E-03 

—4. 484 36764 E-0 2 

5. 59912293E-01 

2.49383911 E+02 

n 

1.32167244E-02 

—3. 9190 974 2 E-02 

4.70780565E-01 

2.72955045E+01 

12 

1. 84165903E-02 

-3.30881266E-02 

3. 77941 921 E-0 1 

-1.38195068E+02 

13 

2.13749670E-02 

-2 . 68956694E-02 

2. 88135857E-01 

-2.42526636E+02 

14 

2.24033374E-02 

-2. 0879 2 29 2 E-0 2 

2. 058659C8E-01 

-2.91979372E+02 

15 

2. 19890953E-0 2 

— 1 . 5227 1 703 E— 02 

1.336C7926E-01 

-2.99054597E+02 

16 

2. 06377118 E- 02 

-1.00 66 547 IE— 02 

7.22493148E-02 

— 2 .781 30674E+02 

17 

1 . 87729674E-02 

-5. 475 637 8 7 E— 03 

2. 159 60 2 30 E-0 2 

—2 .42279200 E+0 2 

18 

1.66994393 E-02' 

—1 . 49364345E— 03 

— 1.91696766E— 02 

-2.015Q2490E+02 

19 

1.46087670E-02 

1 . 87146689E— 03 

-5. 11428886E-02 

-1.62214313E+02 

20 

1.2 60'643 97 E-0 2 

4 .6363 8 191 E-03 

-7, 54764994E— 02 

-1. 27584730 E+02 

21 

1. 07417548 E-0 2 

6.83695614E-C3 

—9. 3 27 17 79 8 E-0 2 

-9 .83 5 17 507 E+0 1 

22 

9. 03222625E-03 

8. 52369422E-03 

— 1.0 554695 6 E-0 1 

—7 .37778267E+01 

23 

7. 480 22923 E-0 3 

9. 7579 805 6 E-03 

— 1 . 1 32 4866 5 E-0 1 

-5.25283102E+01 

24 

6. 08279775 E-0 3 

1.06086897E-02 

—1. 17276720E-01 

-3.33394079E+01 

25 

4.83618762E-03 

1.1 148 81 5 3 E-02 

— 1.1850 219 IE-01 

-1 . 54188429E+01 

26 

3. 73657465E-03 

1 . 14519095 E-02 

-1.17768973E-01 

1.41 45821 9 E+00 

27 

2. 77808997 E-03 

1.15883464E-02 

- 1. 1 587749 9E-01 

1.67979706E+01 

28 

1. 94954076 E-0 3 

1 . 1621 £106 E— 02 

— 1.13554875 E— 01 

3.00159444E+01 

29 

1-.23136027E-03 

. 1 . 1604897 8 E-02 

—1. 11418357E— 01 

4.02329592E+01 

30 

5.94214910E-04 

1.1 57 82 889 E-02 

- 1. 0993 9291 E-0 1 

4.67022969E+01 

31 

1.19329737E-07 

1 . 1566 6698 E— 02 

- 1. 0941 340 5 E-0 1' 

4.89179500E+01 

32 

-5. 9399097, 7E-04 

1.15784734E-02 

-1.09941 48 5E-01 

4. 670 191 57 E+0 1 

33 

-1. 23118019E-03 

1.1 605261 IE-02 

- 1. 11422694E-01 

4.02320725E+01 

34 ' 

— 1. 94943272E-03, - 

1.1622 140 9E-02 

-1. 13561249E-01 

3.00143279E+01 

35 

-2.77808065E-03 

1 . 15890 25 2E— 02 

—1 • 1 588573 9E— 0 1 

1 . 67 9 53403 E+ 01 

36 

-3.73668885E-03 

1.145 27109E-02 

— 1.17778836E-01 

1.41063903E+00 

37 

-4. 83644 820E-03 

, , 1.1149705 IE-02 

—1. 18513354E-01 

-1.54243809E+01 

38 

-6. 08322599E-03 

1.06096249E-02 

l - 1. 172 8 876 9 E-0 1 

-3. 334679 74E+01 

39' 

-7.48084544E-03 

9.7585G899E-03 

' -1. 132 61091 E-01 

— 5.25378040E+01 

40 

-9.03305C65E-03 

8. 5245 546 9E— 03 

-1.05559152E-01 

-7.37897183E+01 

41 

— 1 .07428C 80E-02 

6. 837 6790 7 E— 03 

-9 .3 28 303 6 3 E-0 2 

— 9.83664228E+01 

42 

-1 • 26077426E-02 

4.63689042E— 03 

-7. 5486000 IE— 02 

-1.276C2676E+02 

43 

—1. 46103387E-02 

1 .871678466-03 

-5.11497029E-02 

-1.62236093E+02 

44 

-1.67012931E-02 

-1 .493814836—03 

-1.9 17 275 13 E-02 

—2 .0152 85 86 E+02 

45 

-1.87751033E-02 

— 5.47627903E-03 

2. 159 7 86 8 3 E-0 2 

-2 .’4230 9751 E+0 2 

46 

—2. 0640 10 51 E-0 2 

—1. 00677423 E-02 

7.22573798E-02 

-2. 781650 92E +02 

47 

— 2. 199TA824E-0 2 

' -1.52289975 E-02 

1.33623 5 9 5E-01 

-2.99091122E+02 

48 

—2. 240 600 07 E-0 2 

-2. 08817 55 IE— 02 

2.0589057 4E— 0 1 

— 2.92014672E+02 

49 

-2. 13775249E—02 

-2.6898S433E-02 

2.88170792E-01 

-2.42555611E+02 
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50 

-1, 84187998E-02 

-3 . 30921733E— 02 

3. 779 8 807 8 E- 01 

— 1.33 21 1053 E+02 

51 

-1.32183037E-02 

-3.91957848E-02 

4. 708 3833 2E— 0 1 

2.72999944E+01 

52 

-5.81243829E-03 

-4. 48491961E— 02 

5. 5998121 8E-01 

2.49415844E+02 

53 

3.31323694E-03 

-4.95841015E~.Q2 

6.36640249E-01 

5.09319862E+02 

54 

1. 32901 7-72E-02 

-5 . 2 829 28 2 IE— 02 

6.90746958E-01 

7.73479591E+02 

55 

2. 31491 552E-02 

-5. 39498540 E-02 

7 . 12254264E-01 

9.97501776E+02 

56 

3. 22558814E-02 

-5. 23343576E-02 

6.92874670E-01 

1.13454815E+03 

57 

4. 07069894E— 02 

—4.752 81994E— 02 

6 .279117348-01 

1.14687S49E+03 

58 

4.94320548E— 02 

-3. 93870525E-02 

5. 17717513E-01 

1 .01 731 937 E+0 3 

59 

5.98343705E-02 

-2,. 81983293E— 02 

3. 68304684E-01 

7 . 56304544E+02 

60 

7. 30365966E-0 2 

-1.47074726E-02 

1. 908 2 6646 E— 01 

4.0043 88 89 E+0 2 

61 

8i 90542355E-02 . 

' -7-.18 922196E-17 

4 

0. 

1 .01 90791 1E+00 


STATION 

W20X 

W22X " 

W20XX 

W22XX 

X 

-1. 86109762E— 01 

1 . 692 65480E+00 

2.05070256E-01 

2.00151799E-01 

2 

— 1. 52731005E-0 1 

1 • 57825267E+00 

2.86103633E-01 

-9.80589652E-01 

3 

-1. 16233249E-01 

1. 4007872 1E+00 

3.39572170E-01 

— 2.06167554E+00 

4 

—7. 51 62 92 8 5 E-0 2 

1 .1124 6402E+00 

3. 6449C476E— 01 

—2. 88100759 E+ 00 

5 

-3. 28795391E-0 2 

7 • 5058101 5 E—0 1 

3. 60367627E-01 

— 3.32270110E+00 

6 

7. 30739294E-03 

3. 6142390 1 E-0 1 

•3.28551207E-01 

-3.34856371E+00 

7 

4. 23936411E— 0 2 

—9. 116 787 81 E-0 3 

2.72927332E-01 

-3.00356238E+00 

8 

6. 99735994E— 0 2 

—3 . 240 1 89 8 1 E-0 1 

1.99871953E-01 

-2.394760936+00 

9 

8.84871362E-02 

-5. 60354187E-01 

1. 17502964E-01 

—1.65 669974E+00 

10 

9. 73538482E-02 

-7.1049C625E-0I 

3. 449781 24E-02 

—9 . 1 7067761 E—0 1 

11 

9. 6 964 0941 E—0 2 

—7 * 798 730 2 3 E— 01 

—4. 1 1 793106 E— 02 

—2.72 3 447 80 E— 01 

12 

8. 85344974E-02 

-7. 8 27630 3 3 E-0 1 

-1.03328061E-01 

2 .22 801 749E-01 

13 

7. 38661510E-02 

— 7.3746E63CE-01 

-1. 481 2 93 06 E-0 1 

5.53673738E-01 

14 

5. 506C2047E-02 

-6. 622 62560 £-01 

- 1. 742 5 834 5 E-0 1 

7 , 35 57 3 166 E-0 1 

15 

3. 42456203 E-02 

-5. 72642541E-01 

- 1.825 631 01 E-0 1 

8.0C770026E-01 

16 

1. 33601369E-02 

-4. 800 51015E-01 

-1.7547375 8E-01 

7.86513279E-01 

17 

—5. S9556017E— 03 

-3. 91 79567 8 E-0 1 

-1.56338191E-01 

7.26435355E-01 

18 

-2.26308991E-02 

—3 . 11738192 E— 01 

— 1.28 839047 E-0 1 

6.45 97 8678 E— 01 

19 

-3. 57793838E-02 

-2. 41314955E-01 

-9. 65721 197 E-0 2 

5.61276822E-01 

20 

-4.50759278E-02 

— 1 . 8055239 IE— 01 

-6.27886350EMD2 

4.80367134E-01 

21 

-5.05034242E-02 

-1.28873387E-01 

-3. 02541 605 E-02 

4.05558644E-01 

22 

-5, 23361554E-02 

-8.56152245E— 02 

-1. 16477463E-03 

3.36009854E-01 

23 

-5.10671474E-02 

-5 .027041 72 E-02 

2. 29198844E-02 

2.69901129E-01 

24 

—4. 73306 601 E-0 2 

-2.251 51C95E-02 

4. 11 341 82 9 E-02 

2. 059041 45E-01 

25 

—4. 1820011 IE— 0 2 

-2. 10965406 E-03 

5.33340 850 E— 02 

1 .43903662E— 01 

26 

-3. 5205167QE-02 

1.1248 6807E—02 

6.0C632429E-02 

8 . 50963624E-02 

27 

-2. 80594619E-02 

1 . 80604187E-02 

6. 24345594E-02 

3. 16762890 E-02 

28 

—2. 08050975 E-0 2 

1 • 91106084E— 02 

6.19259728E-02 

-1 .367303 7 8 E-0 2 

29 

-1 . 36860855E-02 

1 . 5495361CE-02 

6.01142322E-02 

-4.83026303E-02 

30 

-6.77335827E-03 

8 . 5 92 649 14 E-0 3 

5. 83896638E-02 

-7.00 2 9573 9 E-02 

31 

-1.28099958E-06 

-9. 4047853 8E-G6 

5. 77030893E-02 

—7.743 42077 E-0 2 

32 

6. 77093 272 E— 03 

—8 . 611 2386CE-G3- 

5.83920030E-02 

-7.00258005E-02 

33 

1.3 6840594E— 0 2 

— 1.5513272CE— 02 

6. Oil 8740 8E- 02 

-4.82947724E-02 

34 

2. 08037045 E-0 2 

—1 • 91 27332 CE-02 

6.19323180E-02 

-1.36605409E-02 



35 

2. 80588861E-02 

-1 . 80 7 537 3 9 E-0 2 

6. 2442 25 7 6 E-0 2 

3 . 16941071E-02 

36 

3. 52055345E-02 

-1.12612G72E-02 

6.00716800E-02 

8. 51201798 E-02 

37 

4.18213642E-02 

2 . 100 28S27E-03 

5. 33425443E-02 

1 . 439 34045 E-01 

38 

4. 7332 95 57 E-0 2 

2 . 2 50 9696 2 E-0 2 

4.U418809E-02 

2. 05941 502 E-01 

39 

5.10702496E-02 

5. 0269 7852 E- 02 

2. 292601 54E-0 2 

2 .69945 740 E-01 

40 • 

5.23398764E-02 

8.56202352E-02 

-1 . 1609855 8E-03 

3.36061974E-01 

41 

5 . 0 50 73 708 E-0 2 

1. 2 88 8 49 3 7 E-0 1 

— 3. 02 53 3963 E— 02 

4.05618630E-01 

42 

4.50797566E-02 

1 . 80571429E-01 

-6.27914184E-02 

4.80435511E-01 

43 

3. 57831641E-02 

2. 41342495 E-01 

-9. 65787385E-02 

5.61354200E-01 

44 

2. 26331847E-02 

3.1 1775305E-01 

-1.28849480E-01 

6.46065401E-01 

45 

5. 99642944E-0 3 

3.91843419E-01 

-1. 56352039E-01 

7 • 26530831E-01 

46 

-1 . 3361 034 3E-0 2 

4. 80110259E-01 

—1 .754901 96E-01 

7.86614998E-01 

47 

-3.42485130E-02 

5. 72713691E-Q1 

-1.82 5 8086 8 E-01 

8.00872407E-01 

46 ' 

-5. 50651512E-02 

6.62345129E-01 

-1.742 75786E-03 

7 . 35666528E-01 

49 

-7. 38730005E-02 

7.3756C732E-01 

-1.481 4448 8 E-01 

5 .53743 81 3 E-01 

50 

-8. 8542 86 82E-0 2 

7 . 82 86C886E— 01 

-1.03338959E-01 

2.22830264E-O1 

51 

-9.69733776E-02 

7.79970600E-Q1 

-4. 11 84C5 91E-02 

-2. 7237 8036 E-01 

52 

-9.73632433E-02 

7.10579642E-01 

3.45006479E-02 

-9. 17181240 E-01 

53 

-8. 84957142E-0 2 

5 . 6042460CE-C1 

1. 17514136E-01 

-1.65690520E+O0 

54 

-6. 99803 93 3 E-0 2 

3.240 60C67E-01 

1.99891365E-01 

-2.39505821E+00 

55 

-4. 23977475E-02 

9. 11 87 63 90 E-0 3, 

2. 729 53992 E-01 

-3.00393556E+00 

56 

-7. 30807043 E-0 3 

-3.61467986E-01 

3. 2 85 833 29 E-01 

-3. 3489 30 16E +00 

57 

3. 28827903E-02 

-7. 5067 35 2 9 £-01 

3. 604028 54E-01 

-3.32311486E+00 

58 

7. 51703 1655-0*2 

-1. 11260164 E+00 

3. 64526 167E-01 

-2.88136709E+00 

59 

1 . 16244673 E-0 1 

-1.4C096086E+QO 

39605659E-01 

-2.06193380E+00 

60 

1 . 52746C58E-01 

"1. 57844865E+00 

2 . 8 61 3 23 67 E-0 1 

-9.80713967E-01 

61 

1. 86128167E-01 

—1. 69286 52 8 E+00 

2.05091820E-01 

2.00174153E-01 


STATION 

NX20 

NY20 

NX22 

NY22 

NXY22 

1 

—4. 27003039E+03 

-1.3515E486E+03 

-4.16901593E+03 

—1 .31 9 59751 E+0 3 

-2. 23 18 72 64 E+0 4 

2 

-4.06421640E+03 

-5.332CC981E+03 

-3.66613286E+02 

6 .26 7 143 77 E+0 3 

-1.95453261E+04 

3 

-3. 37714006 E+0 3 

-8 . 36574437E+03 

3.71615523E+03 

1.40231381E+04 

-1.76081687E+04 

4’ 

— 2.42938575E+03 

-1.0372008QE+04 

7. 4992942 5 E+0 3 

2.11220980E+04 

- 1.4646 281 0E+04 

5 

-1 . 49755248E+03 

-1. 13664915E+04 

1 . 055 5034 4E+04 

2. 6722481 8 E+04 

-1. 08 54 243 8 E+04 

6 

-8.04822677E+02 

-1. 14293 995 E+04 

1. 26172957E+C4 

■3. 00 3 81 5 46 E+0 4 

-6. 5026212 9E+03 

7 

-4. 47337874E+02 

- 1.06792 55 7E+04 

1. 361941 09E+0 4 

3.05571851E+04 

-1.92972325E+03 

8 

-3. 86621 531 E+0 2 

-9. 2597 89 6 9 E+0 3 

1.36579782E+04 

2 . 821 861 19E+04 

2. 4879537 2E+ 03 

9 

-4.98240165E+02 

-7.337C4050E+03 

1. 29205064E+04 

2. 34532354E+04 

6.38553889E+03 

10 

-6. 41 997 742 E+0 2 

-5.09815C89E+03 

1.16178044E+04 

1.70 746541E+04 

9. 4664693 9E+03 

11 

-7. 18598211 E+0 2 

~2. 74443 877 E+03 

9. 942 68943 E+0 3 

1.00716883E+04 

1. 15482566E+04 

12 

-6. 93974291 E+0 2 

•-4. 76566854 E+02 

8.05526466E+03 

3.38241232E+03 

1. 25835847E+04 

13 

-5.91729934E+02 

1 . 525C961QE+03 

6. 08413171 E+03 

-2.27908768E+03 

1 . 26523361E+04 

14 

— 4. 66542021E+0 2 

3. 1192535 5 E+03 

4. 132 8404 9E+03 

-6, 50 572 736 E+03 

1.19298862E+04 

15 

— 3.74207917E+02 

4. 21648400 E+G3 

2. 285 66 176 E+03 

— 9.19084712E+03 

1.06431397E+04 

16 

-3.49923770E+02 

4.78343842E+03 

6. 10979102E+02 

-1.04594777 E+04 

9. 0267872 1E+03 

17 

— 3.99698 834 E+0 2 

4. 8386C283E+03 

-8. 375966 80E+02 

-1.05741349E+04 

7.28916146E+03 

18 

-5.03952428E+02 

4.44179293E+03 

-2. 02199172 E+03 

—9 . 8 46 3 2 641 E+03 

5.59212668E+03 

19 

-6.28820562E+02 

3.68038622E+03 

-2. 921 73 154 E+03 

-8. 571 9 82 6 5 E+0 3 

4. 0448264 3E+03 


H 

UI 

ro 



20 

-7.39598614E+02 

2.65517210E+03 

-3. 53615986E+03 

21 

-8.11515472E+02 

1.46799596E+03 

-3.88517453E+03 

22 

—8. 34986350E+02 

2. 1248030 8E+02 

-4.007C5004E+03 

23 

-8. 148967 C1E+02 

-1. 031 89937 E+0 3 

-3.95298490E+03 

24 

-7 . 65558684E+02 

-2. 20 2 56351 E+0 3 

-3.77944156E+03 

25 

-7 . 041003 85 E+0 2 

-3. 2 542 2 56 8 E+0 3 

-3. 54033885E+03 

26 

—6. 4492 8537E+02 

-4.15758172E+03 

-3.28117.130E+03 

27 

-5. 968 14 115 E+02 

-4. 896883 77 E+03 

-3.03621699E+03 

28 

-5. 62 743 410 E+0 2 

-5.46676731E+03 

-2. 82872040E+03 

29 

-5. 41632867E+0 2 

-5. 8 6 86 2 541 E+03 

-2. 6729766 2E+03 

30 

-5. 30 67 2 591 E+0 2 

-6. 10685814E+03 

-2. 576 976 7 8E+ 03 

31 

-5i 27344328E+0 2 

-6. 1857854 5 E+03 

-2.54462006E+03 

32 

-5. 3067 1976 E+02 

-6. 10695932E+03 

-2.577C5998E+03 

33 

—5. 41 63 1839 E+0 2 

-5. 8 6 874 13 2 E+0 3 

-2. 673141 6 7 E+0 3 

34 

-5. 6 2742352 E+02 

—5. 466 92548 E+03 

-2.82896457E+03 

35 

-5.96813516E+02 

-4. 897C6667E+0B 

-3, 036 535 9 5 E+0 3 

36 

-6. 4492 8 8 50 E+0 2 

-4. 15776801E+03 

-3.28155890E+03 

37 

-7.04101819E+02 

-3. 2 5439189 E+03 

-3. 54078664E+03 

38 

—7. 65 5609 75 E+0 2 

-2. 20 268593 E+03 

-3.77993815E+03 

39 

-8.14898970E+02 

-1.0 3 19 5 62 2 E+0 3 

-3.95351536E+03 

40 

-8. 34987 174 E+02 

2.12506610E+02 

—4. 00759535 E+03 

41 

—8. 11 51 3255 E+02 

1.46812C56E+03 

-3. 88571137E+03 

42 

-7, 39592261E+0 2 

2. 6553895 5 E+03 

-3, 53666070E+03 

43 

-6.28810306E+0 2' 

3. 68069279E+03 

-2.92216522E+03 

44 

-5.03940539E+02 

4. 44216S57E+03 

-2.02232442E+03 

45 

-3.99689844E+02 

4.83901912E+03 

-8. 3779326 5 E+02 

46 

-3. 49 923886 E+02 

4. 7838 5390 E+03 

6. 10953220E+02 

47 

—3. 742 23 470 E+0 2 

4.21685088E+03 

2.28583840E+03 

48 

-4. 665 77 147 E+02 

3. 11952C95E+03 

4.13324601E+03 

49 

-5. 91784143E+02 

1. 52521542E+03 

6. 084 7 8409 E+0 3 

50 

-6. 940 40 987 E+0 2 

-4. 766360 93 E+0 2 

8.05617046E+03 

51 

-7. 186 65421 E+0 2 

-2.744723496+03 

9. 943 8407 9 E+03 

52 

-6.42051 8665+02 

^5.09866C39E+03 

1.1 6191760E+04 

53 

-4. 98272441 E+0 2 

— 7 • 33776436 E+03 

1. 2922C 52 3E+04 

54 

-3.866352 5 1E+0 2 

-9. 26069792E+03 

1. 36596286E+04 

55 

-4.47352913E+02 

t1 • 068C3Q03E+04 

1.36210706E+04 

56 

-8. 04873731 E+02 

-l. 143C5170E+04 

1.261 8847 0E+ 04 

57 

-1.49767934E+03 

—1. 13676044E+04 

1.05563471E+04 

58 

-2 .42 96 178 4 E+0 3 

-1.03 7 30 2 72 E+04 

7. 50024509E+03 

59 

— 3. 37748 139E+03 

—8 . 366 571 5 1E+C3 

3.7166517 8E+03 

60 

— 4.06463 862 E+0 3 ■ 

-5.33254283E+03 

■^3. 6661 24 57 E+02 

61 

-4. 27047979E+0 3 

-1.3517271CE+03 

-4.16948156E+03 


ARE A( Ul, VI, Wl)= 3. 329 26 879 E-0 2 
XTE= -1.37353058E+04 
XTD- 1.40413805E+00 

SENSITIVITY INDEX BY PRESENT THEORY 8= 9.40435287E+00 

AREA1= 4. 1 82 55826E+01 
AR£A2= 1.97430955E+02 
AREA3= 1 . 97430 955 E+02 


-6 .995 73546 E+0 3 2.70807466E+03 

-5.29968743E+03 1.60495442E+03 

-3.60821173E+03 7. 3339 5067E+02 

-2.00035746E+03 7.75432299E+01 

-5.23709516E+02 -3. 8390I715E+02 

7 .93686732 E+02 -6.72925490E+02 

1.93455763E+03 -8. 10558067E+02 

2.88677090E+03 -8. 17329540E+02 

3.64032816E+03 ' -7. 14649243E+02 

4. 18634649E+03 -5 . 26118866E+02 

4. 51 736682E+03 -2 .78160077E+02 

4.62829596E+03 2.57939538E-01 

4 « 51 731 32 9 E+0 3 2.78671261E+02 

4.18623641E+03 5. 26615773E+02 

3.64015558E+03 7. 15 12 1731E+C2 

2.88652720E+03 8 . 17 76653 8E+02 

1.93423176E+03 8.1094719BE+02 

7 . 93265502E+02 * 6.73252689E+02 

-5. 242 41 070 E+02 3. 84150822E+02 

-2.00101533E+03 -7.73908428E+01 

— 3.60901151E+03 -7.33360719E+02 

-5.30064159E+03 -1. 60506197E+03 

-6.99684899E+03 "2 . 70 834973E+03 

-8 ..5 7 3 24672 E+0 3 -4.04529453E+03 

-9. 84771039E+03 -5. 59281010E+03 

-1 .05755778E+04 -7. 29007469E+03 

-1.04608810E+04 -9.0279311 5E+03 

-9. 19207 2 95 E+0 3 -1.06444954E+04 

-6.50660414E+03 -1 . 19314093E+04 

-2.27942751E+03 -1.26539537E+04 

3.38278478E+03 -1 . 25851962E+04 

1.00728985E+04 -1.15497401E+04 

1. 707673 97E+04 ^9 . 46769355E+03 

2 . 345 611 86E+04 -6. 38637825E+03 

2-82220931E+04 -2 .48830576E+C3 

3. 05609648 E+04 1. 92992391E+03 

3.00418799E+04 6. 50339429E+03 

2.67258067E+04 1.08555614E+C4 

2. 112473 82 E+04 1.4648C731E+04 

1 .402490426+04 1.76103309E+04 

6.26794750E+03 1. 95477296E+04 

-1. 3197448 9E+03 2. 232 14739E+04 
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SENSITIVITY INDEX BY BUDI ANSKY THEORY B= 


2.39133287E+00 
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TABLE I.- ELASTIC CONSTANTS FOR MATERIALS USED IN 


IMPERFECTION SENSITIVITY STUDY 


Constituent Young's Poisson's 
material modulus, ksi ratio 

Glass fiber 10,500 0.23 
Boron fiber 60,000 0.3 0 
Aluminum 10,500 0.30 
Epoxy 500 0 . 4o 


Properties of Helically Wrapped Layers (see refs. 35 and 55) 


Material 

Vf 

a 

E x , ksi 

Ey, ksi 

^x 


Gxy, ksi 

Glass -epoxy 

0.65 

0 

7,000 

2,562 

O .283 

0.104 

709 


15 

6,102 

2,419 

0.420 

0.167 

1,054 



50 

5,878 

2,117 

0.649 

0.354 

1,744 

- 


45 

2,270 

2,270 

0.600 

0.600 

2,089 



6 o 

2,117 

3,878 

0.354 

0.649 

1,744 



75 

2,419 

6,102 

0.167 

0.420 

1,054 



90 

2,562 

7,000 

0.104 

O .283 

709 

Boron-epoxy 

0.50 

0 

- 30,250 

2,030 

0.346 

0.023 

525 



15 

23,800 

1,900 

1.179 

0.094 

2,339 



50 

7,762 

1,608 

1.655 

0.342 

5,967 



45 

1,977 

1,977 

0.844 

0.884 

7,782' 



6 o 

1,608 

7,762 

0.343 

1.655 

5.967 



75 

1,900 

23,800 

0.094 

1.179 

2,339 



90 

2,030 

30,250 

0.023 

0.346 

525 



TABLE II.- TEST CYLINDER CONSTITUENT PROPERTIES 


Material 

ERL-2256 

epoxy 

TYP? E 
glass fiber 


Young * s 

Poisson's 

Density 

modulus ) ksi 

ratio 

lb/cu. ii 

635 

o.4o 

0.0445 

10,500 

0.25 

0.0930 


*Cured at 250° F for 2 hours. 



TABLE III.- TEST CXLINDER DIMENSIONS AND EXPERIMENTAL RESULTS 
(inside diameter, 30 in.j length 30 in.) 




a, 

t, 






^exp 

^cr 

Cylinder 

Type* 

deg 

in. 

v f 

v r 

v v 

E , ksi 

Mv , ksi 

Ib/in. 

Ib/in. 








exp 

exp 



1 ■ 

H + C 

15 

0.0796 

53-8 

42.1 

4.1 

4,780 

0.165 

567 

711 

2 

H + C 

15 

0.0948 

55-1 

40.2 

4-7 

4,54-0 

0.233 

804 

1043 

3 

H + C 

30 

0.0722 

47.8 

50.4 

1.8 

3,310 

0.26^ 

360 

511 

4 

H + C : 

; 30 . 

o'. 6717 

47.2 

50.3 

2.5 

3,210 

O.289 

322 

511 

5 ' 

H + C 

4-5 

0.0720 

53.1 

43.6 

3-3 

2,560 

0.288- 

430 

501 

6 

H + C 

45 

O.0697 

55.0 

41.4 

3.6 

2,750 

0.229 

4o4 

501 

7 

H 

15 

0.0943 

45.9 

50.4 

3.7 

5,060 

0/553 

717 

789 

8 

H 

15 

0.0974 

54.5 

4i.6 

3*9 

5,620 

0.484 

711 

1037 

' 9 

H 

30 

0.0765 

43.0 

48.6 

8.4 

3,280 

0.562 

386 

464 

10 

H 

30 

0.0704 

51-9 

44.1 

4.0 

3,500 

0.605 

341 

447 

‘ 11 

H 

45 

0.0670 

57-2 

36.3 

6.5 

2,780 

0.427 

341 

386 

12 

H ' 

45 

0.0669 

52.7 

4i.i 

6.2 

2,180 

0.512 

311 

386 


*H + C -'helically and circumferentially wrapped. 
■H - helically wrapped. 



